VOLUME 84 NO. HY5 
OCTOBER 1958 


PART 1 


JOURNAL of the 


Hydraulics 
Division 


PROCEEDINGS OF THE 


AMERICAN SOCIETY 


OF CIVIL ENGINEERS 


Cl 


BASIC REQUIREMENTS FOR MANUSCRIPTS 


This Journal represents an effort by the Society to deliver inf 
reader with the greatest possible speed. To this end the material 
none of the usual editing required in more formal publications. 


a, 
: ion 


ad 
Original papers and discussions of current papers should be sub 
Manager of Technical Publications, ASCE. Authors should indicate the tech 
division to which the paper should be referred. The final date on which | 
cussion should reach the Society is given as a footnote with each pap 
who are planning to submit material will expedite the review and publ 
procedures by complying with the following basic requirements: 


x 
ey 


1.. Titles should have a length not exceeding 50 characters and spac 


_ 
. 
2. A 50-word summary should accompany the paper. 


3. The manuscript (a ribbon copy and two copies) should be douk 
on one side of 8!4-in. by 11-in. paper. Papers that were originally prepare 
oral presentation must be rewritten into the third person before being 


4. The author’s full name, Society membership grade, and footnote 
stating present employment should appear on the first page of the paper. 


5. Mathematics are reproduced directly from the copy that is su 
Because of this, it is necessary that capital letters be drawn, in black ink, ¥ 
high (with all other symbols and characters in the proportions dict ed 
standard drafting practice) and that no line of mathematics be longer than 6¥ 
Ribbon copies of typed equations may be used but they will be proportional 
smaller on the printed version, 


6. Tables should be typed (ribbon copies) on one side of 8%4-in. by T 
paper with a 6%4-in. by 1014-in. invisible frame. Small tables should be TOU 
within this frame. Specific reference and explanation should be made in the’ 
for each table. 


7. Illustrations should be drawn in black ink on one side of 8Y4-in. by 
paper within an invisible frame that measures 6¥2-in. by 10%-in.; the cap 
should also be included within the frame. Because illustrations will be red 
to 69% of the original size, the capital letters should be ¥%-in. high, Photogr 
should be submitted as glossy prints in a size that is less than 6%-in. by 10) 
Explanations and descriptions should be made within the text for each illustrat 


8. Papers should average about 12,000 words in length and she ile b 
longer than 18,000 words. As an approximation, each full page of typed t& 
table, or illustration is the equivalent of 300 words. } : 


Further information concerning the preparation of technical papers is « 
tained in the “Technical Publications Handbook” which can be obtained | 
the Society. = 


Ey 
Reprints from this Journal may be made on condition that the full title 
the paper, name of author, page reference (or paper number), and dal 
publication by the Society are given. The Society is not responsible for 
statement made or opinion expressed in its publications. ak 
This Journal is published monthly by the American Society of 
Engineers. Publication office is at 2500 South State Street, Ann Arbor, 
Editorial and General Offices are at 33 West 39 Street, New York 18,1 
$4.00 of a member's dues are applied as a subscription to this Journal, Se 
class postage paid at Ann Arbor, Michigan. > a 


Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


HYDRAULICS DIVISION 
EXECUTIVE COMMITTEE 
Harold M. Martin, Chairman; Carl E. Kindsvater, Vice-Chairman; 
Arthur T. Ippen; Wallace M. Lansford; Joseph B. Tiffany, Jr., Secretary 


COMMITTEE ON PUBLICATIONS 
James W. Ball, Chairman; Carl E. Kindsvater; Joseph B. Tiffany 


CONTENTS 


October, 1958 


Papers 
Number 

sedimentation Studies in the Western Gulf States 

by Graham W. Renfro and Charlie M. Moore.........+.2.. selene aLeOO 
jome Experiments with Emergency Siphon Spillways 

DyeWwarren Bi MeBirney.ts (oie) s sia ceteete oe ed es oe eee aad -.. 1807 
ynthetic Flood Frequency 

RR IMATE BERGE larg cake sais [s, Ve myn! 5.9. Aytace, eux bo 086) ale oa ya, OOS) 
Mivining Rods Versus Hydrologic Data and Research 

By Webs Lanebeine ty.) esheets ws Sti eteae a ey Siete aera eens .. 1809 
Jumerical Solution of Flow Problems in Rivers 

by E. Isaacson, J.J. Stoker and A, Troesch...........4.- were | L810 
Jueuing Theory and Water Storage 

IVR ePE oe ARIE DCI. vp ches 0 9.6, 2 oie v0 1d 6-6 overs eb ideie Lote encir, oe ene OLe 
MSCUSSION... <2. 006 GC Oke, CMO TOTO gee Cibiakavakerskovld, hc eke eos 


‘opyright 1958 by the American Society of Civil Engineers. 


Mh = - ot eg we as (0 y sont 


". . ‘ “ye i i 
niga hts th BA ae eh 
ead 5 nee er : 
= he i " 
; 4 di oly, OO 
4 hy * rs 
‘ ; - “7 
‘ 4 ‘a 


way 


aes Sie 
a 


Paper 1806 October, 1958 HY 5 


Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


SEDIMENTATION STUDIES IN THE WESTERN GULF STATES# 


Graham W. Renfro, | and Charlie M. Moore,” A.M. ASCE 
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SUMMARY 


A summary of some of the more important reservoir sedimentation sur- 
veys made by the Soil Conservation Service in the area is presented. An 
analysis is made of existing suspended sediment and reservoir sedimentation 
survey records in a seven-state area which includes 26 physiographic areas. 
Detailed methods used by the Soil Conservation Service in determining sedi- 
ment yields are described. The effect of watershed protection measures, 
including floodwater retarding structures in reducing sediment yield from 
watersheds is discussed. 


INTRODUCTION 


Methods of making reservoir sedimentation surveys have been discussed 
in detail in previously published proceedings\ and will therefore not be 
discussed further. In surveys now being made by the Soil Conservation 
Service emphasis is placed on determination of the factors influencing sedi- 
ment yield. It is important to understand the effect of each factor on sediment 
yield in order to obtain improved bases for predicting yields from unmeasur- 
ed watersheds. This paper summarizes some results of both reservoir and 
Suspended sediment measurements to illustrate how such data are used for 
purposes of planning, design and construction. 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1806 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 5, October, 1958. 


a. Presented at the February 1957 ASCE Convention in Jackson, Miss. 

1. Geologist (Watersheds), Eng. and Watershed Planning Unit, Soil Conser- 
vation Service, Fort Worth, Tex. 

2. Head, Design and Constr. Section, Eng. and Watershed Planning Unit, Soil 
Conservation Service, Fort Worth, Tex. 
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Summary of Reservoir Surveys 


Detailed and reconnaissance sedimentation surveys have been made on 
approximately 150 reservoirs in the Western Gulf area, comprising the states 
of Arkansas, Louisiana, Oklahoma and Texas. Results of 16 surveys, repre- 
senting 16 distinct physiographic areas, are shown in Table 1. These surveys 
were made by either the range or contour methods described by Eakin. 

The following factors have a marked influence on rates of capacity loss in — 
storage reservoirs: q 

1. Capacity-watershed ratio. This is theratio of the original storage of 
the reservoir per square mile of drainage area. In general, if this ratio is 
100 acre-feet per square mile or more a reasonably long life for the reser- 
voir is assured. As an example, Wilson Reservoir has a C/W ratio of 227 
acre-feet/square mile. Its annual rate of capacity loss is only 0.11 percent. 

2. The amount of sediment inflow to the reservoir or sediment yield. 
Several factors determine this yield, which is the function of the amount of 
gross erosion in the watershed and the efficiency of the stream system in 
transporting the erosional material and stream bed load material to the res- 
ervoir site.(4) Other factors which directly affect this yield are: size, sha 
and topography of the watershed, channel density, stream gradients, amount, 
intensity, and frequency of both rainfall and runoff, inherent erodibility of the 
soil, length and degree of land slopes, land use, kind and amounts of vegeta- 
tive cover and effect of conservation practices applied in the watershed. 

3. Trap efficiency of the reservoir. This is expressed as the ratio be- 
tween sediment accumulation and sediment outflow and is dependent on such 
factors as the ratio between storage capacity and inflow, age of the reservoir, 
shape of the reservoir basin, type and method of operation of the outlets, 
which affects detention storage time, texture of incoming sediment and be- 
havior of fine sediment fractions under various conditions. 

In general a reservoir with a high capacity-watershed ratio will have a 
high trap efficiency since the volume of water passing through the reservoir 
will be small. Conversely, a low trap efficiency would normally occur in 
reservoirs with low C/W ratios. Exceptions to both rules can be found. For 
example a high trap efficiency might occur in reservoirs of an arid region 
where spillways seldom, if ever, discharge. In contrast, a reservoir with 
the same C/W ratio in a very humid area might have a comparatively low 
trap efficiency due to frequent, sustained flows through its spillway. Recent 
studies indicate that use of the capacity-inflow ratio (C/I) offers better 
criteria to estimate trap efficiency than the capacity-watershed ratio (C/W). 


Analysis of Sediment Records 


During the years 1954-1955, a detailed study was undertaken by the Soil 
Conservation Service, Engineering and Watershed Planning Unit, Fort Worth, 
Texas in an attempt to bring together in one publication all known existing 
records on rates of sediment production in the Western Gulf eee 

These records include 75 measurements of suspended sediment and 179 
reservoir sedimentation surveys. Table 2 shows the distribution of these 
sediment records among the 26 physiographic areas. 

The following agencies, in addition to SCS, furnished data on suspended 
sediment measurements and reservoir sedimentation surveys. 
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Table 2 
Distribution of Sediment Records 
Western Gulf States 


: Suspended : Reservoir 
Physiographic : Sediment ; Sedimen- 
Area : Measure- ;: tation Total 


: ments : Surveys 


Arbuckle Mountains 

Arkansas Valley 

Basin and Range Province 
Blackland Prairies if 
Cherokee Prairies 

Coastal Prairie 

East Cross Timbers 

Edwards Plateau 

Flint Hills 

Forested Coastal Plain 
Grand Prairie 

High Plains 

Lampasas Cut Plain 

Llano Basin 

Mississippi Alluvial Valley 
Ouachita Mountains 

Ozark Mountains 

Palo Pinto Section 

Pecos Valley 

Pine Flats 

Red Bed Plains 1 
Red Hills 

Sandstone Hills 

South Texas Coastal Plain 
Southern Rocky Mountains 
West Cross Timbers 


ae ee ae oe 


nN 
> 


Ww 
DANK OUAHr RP NUR RP PANWOUWODYFWWOOkYNe 


ANWAONUNORPNNNOLURYNOHURKFONULPUO 
> 


ee 


Ne 
Wr 


DOP NAWFRFOOWORF OR KF RF NUFROWNUONH 


179 254 


~ 
wn 


Total 


ASCE SEDIMENTATION STUDIES 1806-5 


1. Agricultural Research Service, Blacklands Experimental Watershed, 
tiesel, Tex. 

2. U.S. Geological Survey, at Fayetteville, Ark.; Oklahoma City, Okla.; 
nd Austin, Tex. 

3. U.S. Corps of Engineers, at Tulsa, Okla.; Vicksburg, Miss.; Little 
tock, Ark.; Galveston, Tex.; and Fort Worth, Tex. 

4. U.S. Bureau of Reclamation, Amarillo, Tex. 

5. State of Texas, Board of Water Engineers, Austin, Tex. 

Suspended-sediment records were in most cases summarized by the flow- 
uration curve and sediment-rating curve method. Usually at least 25 sedi- 
1ent samples covering a wide range of water discharges in different seasons 
f the year were required. Normally, 10 years of water records were used 
) compile the flow-duration curve for each station. Suspended sediment 
ecords from the Blacklands Experimental Watershed, however, were com- 
uted by the sediment-hydrograph method, as very detailed sampling was 
vailable for these records. 

Adjustments were made in the records for (1) type of sampler used and 
2) unmeasured bed load. Using all of the records selected, the adjusted 
nnual rate of sediment production in tons was plotted against sediment- 
ontributing area in square miles. There was considerable spread in the 
oints on the graph, because of the many factors which were not evaluated, 
uch as land use, slope, conservation practices, relief, channel density, and 
unoff. However, there was a distinct trend toward lower rates of sediment 
roduction with increasing sizes of drainage area. 

Further details concerning analysis and use of these records are contained 
1 reference number 2, “Rates of Sediment Production in the Western Gulf 
tates” by Brune, Maner, Renfro and Ogle. 

Because of the scarcity of records in some of the physiographic areas 
Fig. 1) additional records from the states of Mississippi, Kansas and New 
lexico were included in the study. The Mississippi records were particu- 
irly valuable since they are applicable to the Forested Coastal Plain physio- 
raphic area in Louisiana where few sediment records exist. Because of 
sn¢th of the table giving detailed information on these 254 records it is not 
1cluded here. 

Adjustments were made for estimated trap efficiency of reservoirs. For 
uspended sediment records adjustments were made for (1) type of sampler 
sed and (2) estimates for bedload, which represented an average of 10-15 
ercent of the total sediment load. 

In addition to its value in having all sediment records summarized in one 
ublication, the study provides a means for making general estimates of 
ates of sediment production by physiographic areas. Such estimates are 
ecessary in preliminary design and construction phases of watershed plann- 
ie by the Soil Conservation Service, which includes various types of reser- 
dirs. In addition it provides approximate information to State and Federal 
yencies and consulting engineers for planning municipal, power or recre- 
ional reservoirs and to local agencies on the sediment load of a particular 
‘ream or the capacity loss of a specific reservoir. 

Additional studies covering individual physiographic areas have been made 
1d more are contemplated. These studies are in more detail and are for the 
irpose of determining the relationship between sediment production and 
atershed causal factors such as length and percent of slope, land use, 
jannel density, relief and conservation treatment. 
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The Lid: detailed studies completed thus far are in the Black Prairie and 
Red Hills 8) physiographic areas. These studies permit a more accurate 
prediction of rates of sediment production and distribution of sediment in 
proposed reservoirs. 


Methods Used by the Soil Conservation Service in Western Gulf States 
To Determine Sediment Yields for Design Purposes 


In designing a floodwater retarding structure, it is important that an ac- 
curate estimate be made of the sediment storage requirements. To assure 
adequate useful life of floodwater retarding structures, sediment pools are 
designed to have a minimum life of 50 years. Under-estimation of the 
average annual sediment yield will result in a depletion of the sediment pool 
in less than 50 years, and consequent encroachment on the detention pool. 
Over-estimation results in excessive costs. 

As pointed out previously, the factors which govern the processes of 
sedimentation in reservoirs are known. All these factors can be measured 
or estimated closely. 

Detailed sediment yield studies are designed to attain two immediate ob- 
jectives (1) to determine as accurately as possible, by using erosional 
measurements and estimates, the total annual sediment yield to the proposed 
structure and (2) to determine the distribution and volume of the sediment 
yield retained within the pools of the structures. 

In estimating annual gross erosion the quantity of material derived from 
sheet erosion and the quantity derived from channel erosion are computed 
separately. 


Sheet Erosion 


Sheet erosion is computed by use of the Musgrave(9) equation which is 
expressed as: 


E = FRs!-35 y0-35 P39 1.75 where 


= Sheet erosion, inches per year 


E 
F = Soil factor, basic erosion rate, inches per year, for each soil unit 
R = Cover factor 

S 


= Degree of land slope, in percent 


L = Length of land slope, in feet 


Pp 
30 

Basic erosion rates have been established on all soil units occurring in 
the Western Gulf states. These soils were classified into groups with similar 
family characteristics and similar erosion potentialities. This classification 
- constitutes a useful basis for evaluating erodibility of soils since erosion 
_ rates of many soils, representing large groups of families, have been evalu- 
ated at experiment stations within the area. 

The erosion rates are based on 10 percent slope, 72.6 feet slope length, 
100 percent row crop cover and a maximum 30-minute rainfall intensity of 
1.375 inches occurring once in two years. 

The above set of conditions is adjusted for actual field conditions of slope 


= Maximum 30-minute rainfall, in inches, 2-year frequency. 
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percent, slope length, rainfall intensity, crop rotations and cover which occur 
in the particular watershed under study. 

Tables have been developed showing: (1) basic erosion rates, in inches of 
annual soil loss, for all soil units; (2) adjustment factors for slope percent 
and length; (3) adjustments for various types of cropland cover; (4) adjust- 
ments for all degrees of range and woodland cover; and (5) adjustments for 
rainfall occurrence frequency within the various rainfall belts. 

In order to compute the volume of gross sheet erosion, using the method 
developed by Musgrave and associates, it is necessary to determine from a 
watershed survey the following: 


1. Total acres in each soil unit by percent slope, slope length in feet and 
land use. 

2. Rotations used on cultivated land (percent of each crop grown). 

3. Cover condition classes on all pasture, meadow or woodland. 

4. Land use capability of the above separations. 

5. Maximum 30-minute rainfall intensity occurring once in two years. 

Soil surveys have been made on a high percentage of the total land area in 
each of the Western Gulf states. From these, soil units, slope percent and 
land use are tabulated. Additional field work is required to obtain slope 
length data, crop rotations and cover conditions. 

Rainfall data is based on USDA Miscellaneous Publication No. 204, “Rain- 
fall Intensity - Frequency Data” by David L. Yarnell. 

The computation of volume of material derived from sheet erosion is made 
by use of the tables referred to on page 10. 

A sample calculation in estimating the annual rate of gross sheet erosion 
for a given set of conditions is as follows: . 


Given: 
Cultivated land 767 acres 
Soil unit 1 
Slope 3 percent 
Slope length 350 feet 
Average land use or cover 90 percent row crop; 10 percent 
fall-planted small grain 
Rainfall 1.375 inches (30-minute rainfall 


intensity, 2-year frequency) 


Problem: Compute the average annual rate of soil loss, in acre-feet, 
from sheet erosion. 


Solution: 1. Given, 767 acres of cultivated land. 

2. The basic erosion rate for soil unit 1 is 0.65", annual soil 
loss. 

3. The adjustment factor for a three percent slope on a 350- 
foot slope length, is 0.335. 

4. The land use adjustment factor for 90 percent row crop and 
10 percent fall-planted small grain is 0.93. 

.. The rainfall adjustment factor for a 1.375" Pgq 2-year 
frequency rain is 1.0. 


Computing the product of the above values: 767 (acres cultivated) x 0.65" 


pa gn et, 


ee 
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(basic erosion rate of soil unit 1) x 0.335 (adjustment factor for 3 percent 
slope and 350-foot slope length) x 0.93 (adjustment factor for 90 percent row 
crop and 10 percent fall-planted small grain) x 1.0 (adjustment factor for a 
P3o rainfall intensity of 1.375", 2-year frequency) = 155.32 acre-inches. 
153.00" = 12 = 12.94 acre-feet, annual soil loss. 


The same procedure is used to compute annual soil loss rates for pasture 
and woodland. 


Channel Erosion 


Channel erosion includes gully erosion, valley trenching and streambed 
and streambank erosion. Measurements of soil loss by these processes are 
usually made and recorded in the field. Annual rates of lateral or vertical 
cutting for various types of channel erosion are estimated by comparing aerial 
photographs of different dates, from comparison of cross section surveys, or 
from historical data obtained from landowners or operators. 

The formula for computing channel erosion is given below: 


Dx Lx Rx Nx P : 43,560 = acre-feet of sediment where, 
D = Depth of channel, in feet 
L = Length of channel, in feet 
R = Annual rate of lateral erosion, in feet 
N = Number of banks affected 
P = Percent of bank length eroding 


Road erosion is a significant source of sediment in many areas within the 
Western Gulf states. Sediment production from this source can be computed 
by two methods: 


1. In the same manner as gully erosion. 
2. By use of the formula: 


ye W x L = 43,560 = acre-feet of sediment, where 
D = Depth of roadbed removal, in feet 
Y = Age of road in years 
W = Width of road surface, in feet 


L = Length of eroded section of road, in feet. 


Rates of Sediment Delivery 


Sediment delivery rate, defined as the ratio between annual sediment yield 
and annual gross erosion‘*/ is dependent upon a number of watershed factors. 
Among these are size and shape of the net sediment-contributing area, channel 
density, and the relationship between certain watershed characteristics such 
as topography, which may be expressed quantitatively by the relief-length 
(R/L) ratio. 

Sediment delivery rates have been investigated in detail in two important 
physiographic areas in the Western Gulf states. MANER and BARNES(1) in 
a study of the Texas Blackland Prairies found watershed area to be an ex- 
cellent indicator of the sediment delivery rate within a given size range. In 
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this study, results of detailed reservoir sedimentation surveys were corre- ~ 
lated with detailed studies of watershed sediment sources. Based on these . 


data the average annual volume, in acre-feet, of sediment derived from sheet 
erosion reaching each reservoir was obtained by making adjustments for 
sediment contribution from other sources, and trap efficiency. 

To develop the delivery rate curve, a non-linear regression equation was 
computed by the “least squares” method, using the size of the watershed and 
the measured delivery rate as the two variates. 11) The exponential re- 
gression curves with “X”, or size of watershed, in square miles as the 
abcissa and “Y”, or calculated delivery rate in percent of gross erosion as 
the ordinate, are shown. (See Fig. 2.) The non-linear standard error of 
estimate is 5.06 percent. This is a measure of the variation or scatter about 
the line of regression. 

Three standard errors of estimate, measured plus or minus from the 
arithmetic mean, include 99.7 percent of the cases. Therefore, 3 x 5.06 
percent = approximately 15 percent as the safety factor required to insure 
that in 99.7 percent of the cases the delivery rate will not be underestimated. 

MANER, in a study of sediment delivery rates in the Red Hills physiogr- 
phic area of Oklahoma and Texas'\*/ found that relief and maximum length of 
watershed expressed as relief-length (R/L) ratio in a simple curvilinear 
correlation, or as individual variables in a multiple curvilinear correlation : 
offered a much closer correlation with sediment delivery rates than size of 
sediment-contributing area or several other factors tested in the study. 

Data from 25 reservoir sedimentation surveys were used in developing 
the watershed sediment yield values for the study, which included the follow- 
ing six steps: (1) tabulation of existing reservoir sedimentation survey data 
after making adjustments for reservoir trap efficiency and size of net con- 
tributing area to determine total sediment yield per square mile of contribu- 
ting area; (2) study of rates of gross erosion during the period of record; 

(3) computation of sediment delivery rates by use of data obtained in steps 

(1) and (2); (4) detailed studies of the watershed of each reservoir to assemble 
data on physical characteristics which were thought to influence downstream 
delivery of erosional material; (5) correlation of data obtained in step (4) with 
sediment delivery rates to develop a regression equation for estimating sedi- 
ment delivery rates; and (6) determination of the characteristics of the inde- 
pendent variables used in the equation. 

The delivery rate curve (Fig. 3) was developed by the following equation: 


Log DRe = 2.94259 - .82362 Colog (R/L) ratio 


where DRe = Estimated sediment delivery rate in percent of annual gross 
erosion 


R/L = Relief-length ratio 


The standard error of estimate for the above equation is + 0.04041 log 
units and the coefficient of curvilinear correlation is 0.987, indicating a 
highly significant correlation between these two variables. 5 
Although delivery rate curves have been developed only in these two a 
physiographic areas within the Western Gulf states, they serve as a guide to 
estimate rates of sediment delivery in several additional physiographic areas. 
Proper adjustments must be made for significant differences in various en- 
vironmental factors when using the curves in such areas. n 
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Summary of Sedimentation Data for Design 


After a study of the watershed above a proposed floodwater retarding 
structure has been completed, including computation of annual gross erosion 
for both present and future conditions and application of the proper sediment 
delivery rate to this eroded material, the information is summarized for 
design purposes. 

Gross erosion, by sources, is totaled for both present conditions and the 
amount expected to occur after proper land treatment measures have been 
applied. The land treatment measures will be most effective in reducing 
sheet erosion. The effect of additional measures, such as special stabiliz- 
ation treatment for gully erosion is also considered in estimating future 
erosion from gullies. 

The delivery rate is then applied to both present and future conditions. 
The prevailing rate of sediment yield is used for the estimated installation 
period for regular land treatment measures, normally 10 years. The rate of 
sediment yield for the remainder of the 50-year period is a realistic estimate 
based on the expected amount and effectiveness of proposed land treatment 
and other sediment control measures applied on the watershed. For example, 
it might be determined that 80 percent of the needed land treatment measures 
will be applied at the end of the 10-year period and maintained at 75 percent 
effectiveness. 

The total sediment yield to the proposed structure for the 50-year period 
is then adjusted for trap efficiency. Available data indicates that the trap 
efficiency of floodwater retarding structures of the type being designed by 
the Soil Conservation Service in the Western Gulf states will be 85 to 95 per- 
cent. Cooperative studies by the U. S. Geological Survey and the Soil Con- 
servation Service are now under way to determine the trap efficiency of 
selected floodwater retarding structures. 

The next step is to make an adjustment for the difference in the dry 
weight of soil in place and sediment in the reservoir. Deposits in the de- 
tention pool or in dry pools will be subjected to aeration and will compact to 
a lesser volume than deposits in a permanently submerged sediment pool. 

In computing the volume of sediment in dry pools and in the floodwater re- 

' tarding pools, ultimate compaction of aerated deposits is assumed. Compu- 
tation of volume of sediments below the normal pool level will be based on 
ultimate compaction under submerged conditions. 

Preliminary results of volume weight analysis of samples in Texas, 
Oklahoma and Arkansas indicate that the ratio of sediment volume in a 
permanently submerged sediment pool to soil in place will average 1.3:1.0. 

The final step in summarizing the sediment yield data for design purposes 
is to allocate the total sediment yield to various pools of the retarding struc- 
ture. The distribution of sediment in the reservoir is dependent on the nature 
_ of the sediment, the topography and shape of the reservoir, nature of the 
approach channel above it and the time required to empty the retarding pool. 
_ Allocation of capacity for sediment accumulation in the reservoir is made 
- in both the sediment pool and the overlying retarding pool. In the Western 
Gulf States the proportion of the total sediment yield deposited in the retard- 
ing pool area has been between 10 and 30 percent. 

In the following table a comparison between estimated and measured 
sedimentation rates in Soil Conservation Service floodwater retarding 
structures is shown: 
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Date Estimated Date Measured ~ 

Floodwater Retarding Structure Constructed Rate Surveyed Rate 7 

(Ac.ft./ (Ac. ft./ t 

sq.mi./yr) sq.mi./yr) 
Sandstone Cr. No. 3 - Okla. 1949 1516 1956 1.46 
Sandstone Cr. No. 5 - Okla. 1949 1.55 1957 5 a 
Sandstone Cr. No. 6 - Okla. 1949 1.40 1957 1.98 
Sandstone Cr. No. 14 - Okla. _ 1949 2.23 1956 2.06 
Six Mile Cr. No. 6 - Ark. 1954 1.05 1956 2.09 
Green Cr. No. 1 - Tex. 1954 1.07 1957 K By-4. 

Double Cr. No. 5 - Okla. 1954 0.45 1957 0.51 7 


Measured rates of sediment accumulation in the Sandstone Creek structures 
checked closely with estimated rates particularly in the case of structures 
Numbers 5 and 14. The measured rate in Double Creek, site 5 also compared 
favorably with the estimated rate. The records in Sandstone Creek are con- 
sidered more reliable than those in the other watersheds listed, since they 
encompass longer periods of time. In addition, rainfall averaged about normal 
during the period of record with no unusual storms occurring. During the 
short term record of Site No. 1, Green Creek one 25-year frequency and one 
100-year frequency storm occurred. In Six Mile Creek, also a short term 
record, two storms exceeding a 5-year frequency occurred, both at times 
when antecedant rainfall had been significant. 

These unusual storm events occurring during such short periods of 
record undoubtedly account for the high rate of sediment accumulation as 
compared to the estimated rate. It is believed that the average annual rate 
of sediment accumulation in these structures will closely approximate the 
estimated rate over a long period of time. Periodic resurveys will be made 
to determine rates of sedimentation and capacity loss in the above listed and 
other Soil Conservation Service floodwater retarding structures. 


Effects of Watershed Protection Measures on Major Downstream Works 


The watershed protection projects of the Soil Conservation Service include 
the treatment of the land with a good soil and water conservation program. 

In addition, they include such structural works of improvement as floodwater 
retarding structures, sediment control structures, floodways and floodwater 
diversions. 

The land treatment phase of the program includes such measures as 
pasture, range and woodland improvement practices, improvement of soil 
fertility, terracing and contour tillage. It usually requires five or more 
years, depending on the extent of the job, to get the conservation work done. 

The effects of conservation practices and watershed treatment upon the 
sediment yield of watersheds have been determined in numerous localities 
in the United States by the Soil Conservation Service.(5) Much of the data 
has been brought together in a recent publication. (6) A few of the outstanding 
results are as follows: 4 

The rate of sediment yield to Lake Waco in Texas was reduced by 38 
percent as a result of land use adjustments and conservation measures on 
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less than one-half of its 1,666 square-mile drainage area; conservation 
practices on 35 percent of the 62.3 square-mile drainage area of High Point 
Reservoir in North Carolina resulted in a 24 percent reduction in sediment 
yield to the reservoir; silting of Lake Issaqueena in South Carolina was 
reduced 53 percent following intensification of erosion control practices on 
about 20 percent of the 14 square-mile drainage area; the rate of sediment 
production of the 2.26 square-mile Jones Creek watershed in western Iowa 
was reduced by an estimated 98 percent as a result of multiple practices for 
erosion control and a highly effective sediment retention structure at the 
lower end of the watershed. 

Although there is particular need for more analysis and research on the 
effects of conservation measures and structural treatment upon sediment 
production from watersheds, there appears to be ample evidence that such 
measures and treatment are not only of tremendous on-site significance but 
that they also contribute in a large degree to the reduction of major sediment 
damages downstream. 
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SOME EXPERIMENTS WITH EMERGENCY SIPHON SPILLWAYS# 
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SYNOPSIS 


Two model designs of low head siphon spillways were tested. One is 
patterned after a standard design used for many years as an emergency -type 
of structure for the protection of canals. The second is a proposed design 
for an improved structure to perform the same function. 

Operational characteristics and peculiarities of each are discussed and 
supporting data and pictures are included. The results of one prototype test 
of a standard design are also given. In the conclusions, studies of both de- 
signs are reviewed and future tests are outlined. 


INTRODUCTION 


The Bureau of Reclamation is conducting a program of research into the 
design and operation of low head siphon spillways. The primary function of 
the Bureau siphon spillway is to provide automatic emergency protection of 
‘canals, tunnels, and related structures. In performing its service, installa- 
tions consist of single as well as multibarrel units, some with staggered 
crest levels, with fixed or adjustable crests, and with various air venting de- 
vices. The siphons to be discussed here are of the single-barrel, fixed- 
crest type. 

Numerous siphons have been build along irrigation canals to prevent an 
excessive rise of the water level above the normal operating level caused by 
incoming storm water or the accidental closing of a control gate. Storm 
water usually enters a canal at such a rate that the rise in water surface is 

not rapid; hence, the siphon may operate for some time as a weir before 
priming, if, indeed, it ever primes. The rate of rise immediately upstream 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1807 
is part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 5, October, 1958. 

. Presented at the October 1957 Convention in New York, N. Y. 

. Hydraulic Research Engr., Division of Engineering Laboratories Bureau 
of Reclamation, Denver, Colo. 
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from a closing radial gate in a canal is somewhat more accelerated, particu- ~ 
larly if the closure is accomplished in a short time; therefore, the siphon 
hardly operates as a weir before the excess head is sufficient to bring about 

a prime. 

Siphons are also used in the forebay upstream from the entrance to a 
power penstock. An emergency closure of the turbine wicket gates, following 
powerplant failure, could cause the forebay level to rise very rapidly if the 
storage area were small. The siphon would then be required to prime imme- 
diately to avoid overtopping of the feeder canal banks or damage to the tunnel 
outlet structure if the forebay were tunnel fed. 

Many designs of siphons have been used in the United States, each one usu- 
ally being an improvement over its predecessor. Because of the cost of 
model and prototype studies, however, not many comprehensive analyses of 
the unsolved problems involved have been made. Difficulty in scaling up the 
results of model data, and the attendant lack of confidence in doing so, may 
be partially responsible for the significant lag in siphon technology in this 
country. At the time the studies reported herein were begun, a search of 
foreign literature revealed that Italy and France, and to a certain extent India, 
were more advanced than the United States in siphon development. The first 
two countries, for instance, are today using large siphon spillways for the 
close control of reservoir levels. These structures make use of so-called 
partialization boxes to accomplish an automatic proportioning of air and wa- 
ter in the barrel; the siphon is thus able to adjust for a rising or falling wa- 
ter surface while operating at a partial prime. 

It had long been noted that some of the Bureau siphon spillways were not 
operating as desired. Some operated much above, and others below, the de- 
sign capacities, the differences estimated to be as high as 30 to 40 percent. 
Others failed to prime even when the head over the crest exceeded the depth 
at the crown. In an outstanding instance, one siphon in a pair of identical 
structures built by the same contractor failed to prime while the other func- 
tioned satisfactorily. 

The above facts led to the inception of a research program to learn first, 
how the standard performed, then why it failed occasionally to perform proper 
ly, and later, how to develop a new design for an improved siphon spillway. 

The material which follows lays greater stress on tests of a standard low 
head siphon spillway design than on the proposed design of a seemingly better 
structure. This disproportion in the discussion which follows can be under- 
stood when it is realized that more comprehensive tests were made on the 
standard design in efforts to gain a firm understanding of basic siphon prin- 
ciples; hence, a greater preponderance of data is available, and tests of the 
proposed design are not yet complete. 


Standard Design of a Low Head Siphon Spillway 


Design Considerations 


A general dimensionless drawing of the standard low head siphon spillway — 
is shown in Figure 1. The nomenclature of its various parts as given will be 
used throughout this paper. Losses in the crest region will be minimized if 
at a given section the pressure and velocity distributions correspond closely 
to those associated with vortex flow. After determining a crest elevation and 
the head available, and selecting preliminary values for the center line 
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radius R tb and the throat size D, one may solve for the maximum unit-width 
discharge in the vortex flow formula! 


R 
a=R, yore) Log, “a (1) 
c ‘ 


q = discharge in second-feet for a unit width 
Rc = radius of the crest in feet 
Rg = radius of the summit in feet 
h = atmospheric pressure at site in feet of water 


where 


In practice, a minimum D = 2 feet is recommended, and by using the ratio 
R 

ss 2.5, the coefficient of discharge will be near the maximum obtainable. 
From field information consisting of topographical and hydrological data, D 
and H may be approximated and a value of the coefficient of discharge ob- 


tained from Figure 2. Substituting, the unit q may be solved for in 


47 OD / ARH. (2) 


H = head in feet from crest to tail water elevations 


The unit discharge in (2) must not exceed that in (1). If it does, the siphon 
should be reproportioned and (1) and (2) solved again. The width of the siphon 
can then be determined by dividing the total discharge by the unit discharge. 
For the purpose of setting the crest elevation, an allowed priming head must 
be assumed which satisfies the normal water surface and freeboard eleva- 
tions. The designer’s judgment is an important factor here. 

The size of the siphon breaker pipe is arbitrarily assumed to have an area 
1/24 that of the siphon throat area, or that of the nearest standard pipe since 
it is usually convenient to use pipe sections and fittings for the breaker. The 
inlet elevation of the breaker should be set at, or just below, the normal water 
surface elevation. 


where 


Laboratory Test Facilities 


The test setup is shown in Figure 3. Part A of this figure shows the fore- 
bay, the entrances to two siphon spillways, and water stage recorder box and 
float well. The net forebay area is 144 square feet. Two siphons were in- 
stalled; one on the left, shown in elevation in Figure 3B, was fixed in position, 
and the one on the right was adjustable. The adjustable could be raised or 
lowered by a hydraulic jack to select the siphon desired for operation. Flow 
conditions of two designs thus could be compared readily by alternate opera- — 
tion. Each siphon discharges into a separate tailbox equipped with a 3-foot 
suppressed rectangular weir. 


1. Design Data of Typical Low-head Siphon Spillways, Reclamation Instruc- 
tions, Volume X, Part 2, Design Supplement 3. Coefficient of 0.7 under 
the radical chosen as a safety factor. 
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Fig. 3-A. Model Forebay, Showing Entrances to Two Siphon 
Spillways and Water Stage Recorder Box. 


Fig. 3-B. View of Standard Siphon Spillway, with Piezometer 
Tubes and Manometer Board. 
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Water was supplied to the model from a large volume channel by a cen- 
trifugal pump. The rates of flow were measured through 6-, 8-, and 12-inch 
venturi meters which had been calibrated volumetrically. 

The forebay water surface was measured with a hook gage during constant 
siphon discharges. When a siphon was operating at less than full-primed 
capacity, the fluctuating level was measured and recorded with a modified 
Stevens A-35 water stage recorder. The modification consisted of an electric 
motor and speed reducer which drove the recorder paper at a speed of 1/10 
inch per second. The pen travel was on a 1:1 ratio with the variation of water 
level. An event marker with remote manual control enabled the operator to 
record the occurrence of any significant operational characteristics. 


General Types of Siphon Operation 


In general, a siphon spillway can operate in any one of three ways. First, 
the excess flow to be siphoned away may be low enough to prevent continuous 
operation. The forebay level will rise slowly to an elevation at which the 
siphon can prime, the forebay drops as the siphon operates a short time fully 
primed, and then as the siphon breaker pipe is uncovered, the prime is lost 
and flow over the crest stops. This is the type of operation one may expect 
on a small canal where the siphon capacity is a high percentage of the canal 
capacity. A typical stage versus time curve is shown in Figure 4A. 

Second, the excess flow is great enough to prime the siphon and sufficient 
to maintain a flow through it at a partial prime. In this case, the inflow keeps 
the siphon breaker pipe covered, or lets an intermittent slug of air enter in a 
quantity too low to fully break the prime, and the water level within the barrel 
remains below the summit. The siphon will continue to operate in this stage _ 
until there is a change of inflow. An increase will bring the siphon closer to 
a fully primed state; however, a decrease may cause a complete break in 
prime and temporary cessation of flow. The stage versus time curve appears 
as in Figure 4B. 

In the third type of operation, the excess inflow must reach the siphon in a 
comparatively short period of time at a rate which causes it to prime imme- 
diately and to remain primed. The siphon barrel will then be full from inlet 
to outlet. The forebay water surface elevation can be lower than the summit 
elevation for this type to exist. Pressures in the throat section reach maxi- 
mum negative values when the barrel is flowing full and the difference in ele- 
vation between crest and forebay water surface is a steady maximum. Any 
increase in excess flow decreases the negative pressures at the throat. This 
type of operation is most apt to occur during a rejection of flow at a power- 
plant. Figure 4C shows the stage versus time curve. 

Of course, there are infinite variations of the above stages of operation, 
but for purposes of this study, the three general ones defined were adequate. 


Operating Characteristics 


The standard siphon spillway design chosen for study in the laboratory was 
of the type constructed on the Mohawk Canal, Gila Project, in Arizona. Pro- 
totype dimensions were reduced on a 1:4 scale. The resulting siphon is dimen- 
sioned in Figure 5A. The locations of 50 piezometers utilized in obtaining 
pressure distribution patterns are shown in Figure 5B. 

In Figure 5A, when the forebay inflow was sufficient to top the crest, the 
overflowing sheet followed the profile to the deflector lip. When the rate of 
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flow was still a very small amount, enough energy was available to allow an 
unbroken jet to leave the lip in a direction tangent to it, cross the barrel, and 
strike the opposite boundary. Along that line of impact, the jet broke into 
flows in all directions, though the principal one was downward. An effective 
seal was accomplished at the point of impact which divided the barrel into 
two parts, an upper air cavity and a lower one. As the head in the forebay 
rose, air in the upper cavity was forced through the pressure equalizer to 
the underside of the nappe. Air flows rather freely to the lower air cavity 
because of an air demand at this point created by the jet of water. In fact, a 
better equalization would result if it were possible to install a larger diame- 
ter pipe. The larger it is, the quicker the siphon will prime. 

At this stage of operation, any but a very slow rate of rise in the forebay 
produced a positive pressure in the lower air cavity. The pressure continued 
to rise only until it was great enough to break through momentarily the free 
jet of water, thereby restoring the equality of pressures in both cavities. 
When this happened, priming was further delayed. 

Following impact of the jet of water against the opposite boundary most of 
the water fell immediately into the lower bend pool where water had sealed 
at the roof of the bend. The air in both cavities must be evacuated by means 
of turbulent mixing past this seal far enough to rise of its own accord to the 
free water surface at or near the exit of the outlet tube. Our experience with © 
this design showed that considerable turbulence was required in the lower 
bend pool to begin to move air out, and until that degree of turbulence did 
exist the forebay head continued to increase. While the reverse curvature of 
the lower leg aided materially in obtaining an adequate seal and separation of 
air cavities, it also caused water to move toward the lower bend with a direc- 
tion component opposite that conducive to discharge through the siphon. 

The design tested was the outgrowth of studies on vertical drop siphon 
spillways in which similar crest profiles fed a vertical conduit with a pool 
and sharp bend at the lower end. The amount of drop was merely varied to 
fit the particular site of the structure. Consequently, only by chance did the 
required discharge, operating head, crest profile, length of drop, and shape 
of lower bend combine hydraulically to operate correctly. The reverse curve 
and the deflector lip of the present design did much to improve the former 
design. 

Very early in the investigations it became obvious that the meaning of 
“prime” must be positively defined. Many of the technical papers studied as 
a prelude to this work made reference to the time required for priming of a 
siphon, with no criterion established for determining what constituted a prime 
other than visual observation. Therefore, a primed siphon was defined as one 
in which the pressure at the crown or summit had reached a maximum nega- 
tive value for the particular forebay inflow. So defining the term allowed us to 
measure with reasonable consistency the length of time required for a siphon 
to reach a degree of prime in keeping with the inflow. The negative pressure 
does not reach the maximum possible value until the forebay inflow reaches 
the capacity discharge of the siphon at a constant head. 

The maximum discharge for the model was 5.70 cubic feet per second; a 
full prime was reached at 5.35 cubic feet per second. The coefficient of dis- 
charge, based on the elevation difference from forebay water surface to the 
outlet roof, was 0.84. This value is high compared to those measured on a 
some prototype siphons of similar configuration, but the model was built un- 
der laboratory controlled standards and tested in an ideal setup, two factors 
conducive to a higher coefficient. 
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PIEZOMETER PRESSURES 
Siphon Spillway Model--Standard Design 
Scale 1:4 


Q = 5.35 cfs 
: Pressure in: : Pressure in: : Pressure in 


Piezometer: ft of water:Piezometer: ft of water:Piezometer: ft of water 
: 0 F 


. : 9 : -0.05 : LS 

1A -0.98 : 9A : 0.79 : 15a 2 -1.69 
1B 1.49 : 9B : -0.91 : 16 y (e217 
2 0.73 ‘ a k0 2 91.13 > 16A : -1.90 
3 0.69 s) SOA : -0.61 : 163 2 «-3.34 
3A -0,.81 : 1OB : 1.89 Memes Blips 2 -2.14 
3B O.TT : 11 : -1.28 4 17B : -3.32 
\ 0.56 = 11A : -0.65 : 18 $ -1.96 
5 0.35 man: Se eke : 18A 2 <1.45 
5A -1.16 3 12 c -1.53 . 18B G -3.58 
6 0.09 sie OR : -0.94 : 119A : 1.22 
if SOG : 13 : -1.77 : 19B : -2.88 
TA -1.63 : 113A : -1.22 : 20A : -0.62 
7B 0.26 ue pale 2 «-1.94 : al : -2.05 
8 -0.47 : 14a ‘ -1.47 $ 22 3 -2.90 
8a -1.27 $ 14B - -2.89 : 23 : -2.90 
8B -0.23 : : : 


Table 1 gives the piezometer pressures measured for a fully primed con- 
dition at points shown in Figure 5B. Along the roof of the model all pres- 
sures were negative, with reference to their respective elevations reaching a 
maximum value just downstream of the siphon breaker inlet. On the side cen- 
ter line, pressures were negative from the inlet to a point beyond the PC of 
the lower bend, and then increasingly positive through the lower bend. Pres- 
sures over the crest and down to the PC of the lower bend along the floor were 
higher negatively than those in the other piezometer sets. Maximum pres- 
sure occured at Piezometer 18B, about -3.6 feet of water. 

Since the design was based on the assumption of vortex flow at the crest, 
it is of interest to learn whether vortex-type flow does exist. With a flow of 
5.35 second-feet, physical dimensions of the model were used to compute the 
head as given in the equation for vortex flow,2 


2 2 
R -C 
Q= V2gh 2,3026b | log(R + C) - log (R - C)| (3) 
2 Ro 
where Q = discharge in second-feet 


R = radius of barrel centerline at the crest section 
c = height of barrel at the crest section 
b = width of the barrel 


The value of h so obtained was 
h = 1.87 feet 


The corresponding head measured in the model and adjusted for elevation 


2. Addison, Herbert, “A Treatise on Applied Hydraulics,” Third Edition, 
page 168. 
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differences, since the siphon was in a vertical instead of horizontal plane, 
was 5 


h = 1.94 feet . 


These two values are close enough to justify acceptance of the vortex flow 
assumption. 

Figure 6 is a composite of two sets of basic operational curves, one of 
maximum and minimum crest heads for the initial cycle of operation, and the 
other of crest heads for those cycles succeeding the initial cycle. About 
0.125 foot of head at a discharge of 0.1 second-foot was required to prime the 
siphon. For discharges to 1.75 second-feet, operation is characterized in 
Figure 4A. The crest was dry between cycles of operation and each cycle 
was therefore independent of the one which preceded it. 

From 1.75 to 3 second-feet, operation was as shown in Figure 4B. The 
initial cycle resulted in a higher head to prime the siphon, but succeeding 
cycles operated over a reduced range of heads. Some varying depth of water 
was continually flowing over the crest. The closer the discharge was to 3 
second-feet, the narrower became the range of operating heads. 

About 3 second-feet, the maximum crest head on the initial cycle continued 
to rise. At about 4.5 second-feet, the maximum head equaled the barrel / 
height at the crest section. For cycles following the first one there was no 
difference in maximum and minimum crest heads, and operation was as shown 
in Figure 4C. As the capacity of the siphon was approached, the head curve 
for cycles following the initial one rose rapidly toward the head curve for the 
initial cycle. 

The curve of priming time versus discharge for the initial cycle of opera- 
tion only is shown in Figure 7. The shape of the curve, particularly the re- 
versal characteristic, can be explained as follows. Below 2 second-feet, the © 
rate of forebay rise was not sufficiently high to cause’a significant compres- 
sion of air in the upper air cavity; however, over 2 second-feet, the rate of 
rise became high enough to cause the pressure in the cavity to be forced 
above atmospheric before air pumping in the lower leg began. The result was 
to delay priming. Also, up to 2 second-feet the lag between the start of air 
pumping and rise in forebay was not as serious as it apparently becomes 
above 2 second-feet. The two explanations offered here are closely related, 
of course. 


Prototype Test of Standard Design 


The results of one- prototype test were available for use in this investiga- 
tion. A siphon spillway on the Mohawk Canal, Figure 8, closely resembling 
that modeled in the laboratory, was used to obtain pressure and velocity pro- 
file data. The barrel was of constant section, 2 feet by 4 feet. At five station: 
along the roof of the structure, pitometer outlets were utilized to obtain . 
velocity profiles across the barrel. | 

Designed for a capacity of 150 second-feet, the structure passed a maxi- | 
mum of 136 second-feet. The head assumed to be required for a prime was 
0.17 foot, but in the tests 0.46 foot proved to be necessary. 

In order to show whether vortex flow prevailed within the siphon, Table 2 
was prepared. Measured velocities were multiplied by corresponding radii to 
determine how nearly constant the products were. An average VR value was 
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Table II 


VARIATION OF VR ® k WITH RADIUS IN PROTOTYPE 
SIPHON eee aE THE oo SEPTION 


Radius is AO ere 2.3.6 : 3. 
Velocity :21. 3 :21.3: :20. 25: 219. 31: 117. 70: "16. a 15. 65: a W3: 13. "6: alee "31: an ‘6h 
VR 763.9 :68.2:68.8 :69.5 :67.3 :65.7 :65.7 :63.5 :63.3 :63.9 : 58.2 


Deviation from:-3. a +3.0:43.9 345.0 :41.7 :-0.8 :-0.8 :-b.1 :-h.b :-3.5 i-12. 1 
VR = 66.2* 3 3 3 : 3 $ 3 : O 


*Average of all VR values, excluding those at boundaries. 


computed neglecting boundary measurements, and the deviation from that 
value was noted. Boundary measurements were not included in the average 
because wall friction upset the measured velocities. Deviations from the 
average ran from +5.0 percent to -4.4 percent. This small range was con- 
sidered adequate proof that near-vortex type flow exists, and it served to lend 
further credence to the model results. 


Analysis of Model Operation, Standard Design 


Several important conclusions were derived from studies of the standard 
design. Although this design assures positive separations of air volumes and 
seals from the atmosphere, the practice of bending back the lower leg to ac- 
complish these objectives causes a serious delay in the removal of air from 
the barrel. This being true, the priming, though positive in action, is re- 
tarded, and the crest head required for priming is higher than desirable. 
Further, the built-in lag in priming produces greater head fluctuations than 
can be tolerated in most installations. 

Numerous attempts were made to design special air-intake partialization 
devices for the model to replace the siphon breaker pipe. Our efforts were 
unsuccessful mainly because of the delay in the start of priming and the com- 
paratively rapid rate of forebay rise. During the rise in the forebay water 
surface, the devices became submerged before sufficient subatmospheric 
pressure in the crown of the siphon could develop to draw in the needed pro- 
portion of air; and when the water surface dropped, siphon action was broken 
before the level reached the minimum stage desired. Thus, it was considered 
impossible to modify the model design to make use of a partialization device. 

Limited success in prototype structures has been had with air slot par- 
‘tializers substituted for the siphon breaker pipe, but where some degree of 
partialization was attained the rates of forebay rise were very much slower 

- than in our model tests. In one instance where the rate of forebay rise was 
very rapid, the air slots failed to accomplish partialization at all, due to flood- 
ing. 

Although erosion downstream of the outlet of a siphon spillway was not 

studied with our models, in many locations the prime-break-prime type of 
intermittent operation, which is a characteristic in the lower third of the dis- 
charge range, could not be withstood. 

The design tested is generally a satisfactory one when used where the rate 
of rise of the forebay water surface is low and erosion damage due to prime- 
break-prime cycling can be tolerated. Where the rate of rise is high, free- 
board allowances in the forebay must be increased to accommodate the 
greater range of heads which must result. 

The model data yielded no positive clue as to why some prototype structures 
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operated satisfactorily while others failed to do so. The answer to this prob- 
lem must be sought in prototype tests on existing structures. 


ee 


Proposed Design of a Low Head Siphon Spillway 


Design Considerations 


With the results of tests on the standard design of a siphon spillway in 
mind, a simplified siphon with more consistent and predictable characteris- 
tics was sought. Several objectives were established as guides for an im- 
proved design. One, a new siphon spillway should continue to make use of a 
vortex flow principle of earlier designs. Second, the crest head required for 
priming should be reduced. Closely related, third, was the desire to reduce 
the priming times for all heads. Fourth, if the second and third objectives 
could be attained, efforts should be made to fit the structure with some de- 
vice for the partialization of air, thereby allowing the siphon to reach some 
degree of prime at low discharges without the necessity of becoming fully 
primed on the first cycle and then settling back on succeeding ones to a par- 
tialized state. It was hoped, too, that a forthcoming design would be less dif- 
ficult to construct and could perhaps be simpler than the standard design. 

Following library research, it was decided to pattern the new siphon after 
ones built and successfully operated in Europe. The design of the structure 
could not be proportioned much differently than the original, i.e., the ap- 
proach hood and crest section must resemble the former design and if test 
results were to be easily compared, the total siphon head must be the same. 
But the conduit downstream of the crest as well as the lower bend were con- 
siderably changed in adapting the European designs to our use. 

Figure 9 shows the proposed siphon spillway model design which evolved 
from these studies. The slope of the lower leg as well as the positions of the 
deflectors were arbitrarily chosen. Each deflector could be rotated from the 
flush position to that of a 45-degree projection. Divided arcs were installed 
to allow setting the deflection within 1 degree. To obtain a positive seal in 
the bucket, the downstream lip was constructed abently higher than the down- 
stream end of the sloping roof section. 

Since the maximum negative pressures near the crown of the standard de- 
sign occurred some 15 degrees downstream from the crest, the outlet of the 
siphon breaker pipe was in the proposed design so positioned. The reason for 
choosing this location was that subatmospheric pressures probably begin 
developing in that location earlier than elsewhere. 

No extensive sets of piezometers for the measurement of pressures were 
installed in this siphon. The pressure patterns and magnitude of pressures 
were sufficiently developed in tests of the standard design, and a repetition of 
this work was not thought to be necessary, despite the different configuranigs 
of the exit section of the structure. 


~~ 


The Model 


The model of the proposed design was installed adjacent to that of the 4 
standard design, Figure 10. The vertical position of the siphon had to be ad- 
justable so that only one siphon would operate at a given time. For this rea- _ 
son, the mounting plate, held between two greased plates, was allowed full 
freedom of movement from below to well above the position of the standard 
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Fig. 10. Model of Proposed Siphon Spillway Design. 


siphon. A hydraulic jack was used to change elevations of the model. A tail- 

box was built to accommodate the discharge over a range of positions. 
Anticipating the need for measurement of the rate of air demand, the 

siphon breaker pipe was provided with an orifice plate section and accompany- 

ing piezometers. The position of the orifice will allow use of a partialization 

device upstream without further modification. 


Operating Characteristics 


Referring again to Figure 9, when water crossed the crest, the overflowing 
sheet accelerated down the slope toward one of the projecting deflectors. The 
direction of flow was then altered to throw the water in an unbroken sheet 
across the conduit to the roof of the descending leg. The position in which it 
initially struck the roof was dependent, of course, on the deflector angle. ; 
Though the jet was broken into flow components in all directions at the line of 
impact, the principal direction remained downward into the bucket pool. Im- — 
pact with the roof and turbulence in the pool served to mix air and water well. 

With the bucket full of water, an air cavity from forebay water surface un- 
der the hood to the bucket pool had to be evacuated before a full prime could 
be realized. If the elevation as well as angular position of the deflector were 
correctly chosen, a jet of water struck the pool very near the lowest edge of 
the roof. The turbulence thus created and the forward component of velocity 
combined to produce immediate air pumping. 

Two deflectors were built in the model to enable evaluating the relative 
merits of each position. Though test results showed the lower position more 
effective than the upper one, it would be unwise to exclude the possibility that 
a third undetermined position might be better than either of the first two. 

Results discussed here are limited to those obtained from the lower de- 
flector. 

The proposed design had a maximum capacity of 4.35 second-feet; a full 
prime was reached at 3.80 second-feet. The coefficient of discharge, based ~ 
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on the elevation difference between forebay water surface and the downstream 
lip of the bucket, was 0.62. So computing the coefficient does not give recog- 
nition to the height the discharging water rises above the lip of the bucket, but 
even if the total head were reduced by that increment, the rise in Cg would 
not bring the value above about 0.69. 

Thus, the coefficient of discharge for the proposed design is much below 
that of 0.84 for the standard design. It can be assumed that part of the re- 
duction was caused by back pressure exerted by the turbulent head of water 
in and above the bucket. As will be seen in the characteristics that follow the 
reduced Cq can be tolerated, or compensated for by a longer crest length, 
when considered with operational improvements. 

Figure 11 is the companion set of operational curves to Figure 6 on the 
standard siphon. About 0.035 foot of head at a discharge of 0.1 second-foot 
was required to prime the siphon. Up to 1.15 second-feet, operation is shown 
in Figure 4A. Each cycle was independent of preceding ones and had no in- 
fluence on those which followed. The crest was dry therefore in the time in- 
terval between cycles. 

From about 1.15 to 1.60 second-feet, the operation is as shown in Figure 
4B. The crest head rose to a higher value on the initial cycle than on any 
following cycles. Water never quite ceased to flow over the crest. The 
closer the discharge was to 1.6 second-feet, the narrower the range of operat- 
ing heads became after the initial cycle. 

Above 1.6 second-feet, the maximum crest head on the initial cycle con- 
tinued to rise, but for cycles following there was no difference in maximum 
and minimum heads. As the discharge approached 4.25 second-feet, the crest 
heads for all cycles became nearly the same. Between 4.25 and 4.35 second- 
feet, the crest head increased from 0.425 to 0.5 foot and the siphon capacity 
was reached. Operation above 1.6 second-feet is as shown in Figure 4C. 

The curve of priming time versus discharge for the initial cycle of opera- 
tion with the deflector set at 45° is shown in Figure 12. The displacement 
which took place about 1.15 second-feet can be explained as follows. At all 
discharges above this value, the rate of flow over the crest during priming 
was sufficient for a very short period of time to break the water seal at the 
low point of the roof. In other words, the turbulence at that location reached 
such a magnitude as to depress the solid water surface below the sealing edge, 
allowing a sudden slug of air into the air cavity of the barrel. The resulting 
momentary reduction of negative pressure in the cavity caused a delay in 
priming while the invading air increment was pumped out and the priming 
process was resumed. 

In the low range of discharges, up to about 0.5 second-foot, the crest head 
required for priming increased as the deflector angle decreased. If the angle 
were less than 20 degrees, priming became erratic and unreliable, the siphon 
yeing more apt not to prime than it was to prime. 


Analysis of Operation and Comparison of Standard and Proposed Designs 


The proposed design of siphon spillway operated with much greater con- 
sistency than the standard design. In obtaining characteristics curves, it was 
lecessary to average a greater number of points to obtain the plotted ones in 
he case of the standard design; with the proposed design the spread of points 
0 be averaged was much narrower. 

The manner in which the crest heads increased with decreasing deflector 
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angles, raises a question as to whether the heads for 50 or 55 degrees might 
be lower yet. There must be a practical limit, of course, and we were unable 
to increase the angle accurately without rebuilding the section containing the 
deflectors. 

The principal advantage of the proposed design over the standard design is 
the significant reduction in the time to prime. The lag between the time 
water flows over the crest and air pumping begins having been nearly elimi- 
nated, efforts to fit the new siphon with a partialization device may meet with 
more success than was had with the standard siphon. 

The new design also is capable of priming at a lower discharge than re- 
quired for priming the standard design, and it is capable therefore of partial- 
izing earlier in the discharge range, even with a siphon breaker pipe. 


Future Investigations 


The studies conducted to date are by no means complete. Attention will 
be concentrated on providing the proposed siphon design with a partialization 
device which will proportion satisfactorily the intake of air in the crown. — 

Many different avenues for investigation have occurred to us which should 
be explored to make the studies complete. The slope of the downstream leg 
of the proposed design should be varied, as well as the position along the 
slope of the deflector. The shape of the bucket and elevation of the down- 
stream lip undoubtedly influence the total operating head; a separate study 
could be made of these features. 

It is believed that the new design, even with conventional siphon breaker 
pipe, is an improvement over that formerly used. Efforts will be made to 
build a prototype structure with test facilities for the purpose of developing 
characteristics at the larger scale. 
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Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


SYNTHETIC FLOOD FREQUENCY 


Franklin F. Snyder,! M. ASCE 
(Proc. Paper 1808) 


SYNOPSIS 


A procedure is developed for computing the flood discharge probability 
associated with a given rainfall-duration-frequency pattern on natural drain- 
age basins, non-channelized overland flow areas and areas with storm sewer 
drainage utilizing basin runoff-producing characteristics of area, length, 
slope, friction and shape. The approach is patterned after the so-called 
rational method and utilizes the time of concentration concept with a unit hy- 
drograph interpretation, but recognizes and evaluates separately the effect 
of storage existing in all types of channels or conduits and an average rainfall- 
runoff relation. The variable factors in this case have been evaluated for 
application in the vicinity of Washington, D. C. 


Time of Concentration 


The time of concentration, Tc, as used herein, is based primarily on hy- 
drograph analyses(1) and can be defined as the time from cessation of effec- 
tive rainfall to the inflection point of the recession side of the resulting run- 
off hydrograph. It is essentially the same as the time of concentration de- 
fined as the time for runoff from the remote portion of the drainage basin to 
reach the point of interest. As such it constitutes the time base of a 
area-shape diagram as visualized in some unit hydrograph analyses. 2)(1) 
(See Fig. 1) 

Table I gives the basic data for some of the drainage areas utilized in 
evaluating the coefficients and exponents of the basic equations, 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1808 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 5, October, 1958. 
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s 0.6 
o WS 10 } 
L' = o.1ve °F a (2) 


in which T, is the time of concentration, C; is a coefficient dependent on the — 


type of drainage system, L' is the length of an equivalent channel or conduit 
having the same time of concentration but with a standard slope of 1% and a 
standard friction factor of 0.1, Lis the length of the principal channel, s is 
‘the weighted slope of the principal channel, in%, determined as the mean 
height of the channel profile above the point of interest divided by one-half [ 
the length, and n is the friction factor for the channel or conduit.(2) Fig. 1 of ; 
Appendix I provides a graphical solution of Eq. 1. Column 8 of Table 1 ‘ 
presents observed values of C; for the drainage areas listed, which canbe ~ 
compared with the values of C; as given in Appendix 1 for use in the vicinity 

of Washington, D. C. 

The general form of Eqs. 1 and 2 were assumed, and varied and successive 
plottings of the data in Table I were utilized to evaluate the exponents and 
coefficients. Assuming slope to the 0.5 power, values of L' were plotted 
against values of T, on logarithmetic paper. The best exponent was slightly 
greater than 0.6 for natural areas, slightly less than 0.6 for the airfields and 
about 0.5 for the sewered areas. A value of 0.6 is consistent with that de- 
veloped by the author(2) from a study of a large number of unit hydrographs 
and was adapted for use for all three types of drainage. Values of length to 
the 0.6 power divided by T, were also plotted against values of slope. The 
best exponent was about 0.4 for natural areas and not definable for the air- 
fields and sewers, partly because of the limited range of slopes available. 

By including the slope factor in the expression for L', an actual value of 0.3 
results for the slope-effect exponent. Use of this value for slope effect did 
not weaken the overall correlation for the airfield and sewered areas. Aver- 
age, equivalent and weighted values of slope were tested, and it was found 

that the weighted form, as adopted, gave the best correlation. Consideration 
was also given to the use of length to the 0.6 power rather than equivalent 
length, L' , but the latter gave slightly better results. A value of 0.10 was 
used as the standard friction factor, n, in converting actual lengths to equiva- 
lent lengths, L'. An overall value of about 2.0 for C; for natural basins was _ 
evident but a value of 1.7 was selected to be slightly conservative for use in 
the vicinity of Washington, D. C. on the basis of values for the natural basins 
in Table I. 

The value of C; of 1.7 hours per mile for natural drainage basins in the 
vicinity of Washington, D. C., is four times larger than that for completely 
sewered areas. The expression for L takes into account all the factors of 
the Manning equation except the hydraulic radius which is indirectly recog- 
nized in the expression for T, through the exponent, 0.6. The question thus — 
arises as to why the value of C; for natural basins should be so much larger 
than that for sewered areas. 

It is believed that this difference, as well as much of the differences be- 
tween values of T, for various natural basins, results from differences in 
the amount of subsurface storm-flow that exists in representative hydrographs 


2. The letter symbols in this paper are defined where they first appear in _ 
the text and are given in Appendix II for ready reference. 
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of direct runoff for the areas. Direct runoff can be defined as that part of the 
total runoff appearing systematically in the drainage channels in direct re- ‘ 
sponse to rainfall or melting snow in contrast to ground-water or base flow 
which enters the channels from the main zone of saturation. Direct runoff 
thus includes all or most of the true surface runoff as well as an appreciable 
portion of the subsurface storm-flow which is characteristic of the area. 
Subsurface storm-flow can be defined as that porticn of the storm flow which 
penetrates the surface of the ground temporarily but soon reappears or drains 
away from the area through the upper soil layers more rapidly than true 
ground-water seepage. 

As such, subsurface storm-flow is subject to greater storage effect and 
time delay than is true surface runoff. The direct runoff hydrographs for a 
drainage basin thus reflect some subsurface as well as surface storage 
characteristics. 

Except for areas of 100 per cent imperviousness, storm runoff hydrographs 
from urban areas with complete storm drainage systems contain some sub- 
surface storm-flow from yards, banks, stone driveways, etc. For such areas, 
however, the possible range in the relative amount of subsurface storm-flow 
is much smaller than that for natural drainage basins. After a value of C; 
for natural basins has been selected for a particular region, it is considered 
that the percentage of open drainage channels which have been eliminated 
and the percentage of the basin which has been storm sewered are the prin- 
cipal factors to use for interpolating between the selected value of Te for the 
natural basin and the value of 0.42 for the complete storm drainage condition. 
Data are not available to define the relationship exactly, but those which are 
available indicate that reasonable results are obtained by averse the two 
percentage factors and interpolating on a straight line between C;’s (1.7 in 
the case of the vicinity of Washington, D. C.) and 0.42. Thus, if §0 per cent 
of a basin has storm sewers, but only 40 per cent of the natural drainage 
channels have been eliminated, a value of C, of 0.9 is appropriate, that is, 

1.7 minus 0.80 + 0.40 (1.7 - 0.42). 
2 


Rainfall-Duration- Frequency 


Fig. 2, Appendix I, shows the rainfall-duration-frequency curves of point 
rainfall recommended for use in the vicinity of Washington, D. C. These 
curves are based on data published by the U. S. Weather Bureau in Technical 
Paper No. 25. 3 

The rainfall intensity-duration-frequency curves of T. P. No. 25 were 
adjusted, as recommended therein, to give partial-duration values and the 
rates were converted to amounts. 

The curves presented in T. P. No. 25 for Washington, D. C. were derived 
from the annual series (maximum value for each year) for the periods 
1896-1897 and 1899-1953, and are spaced according to the Gumbel procedure. 

Also shown on Fig. 2, Appendix I, are curves giving factors to apply to 
values of point rainfall to obtain average depth of rainfall over drainage basins - 
of appreciable size. : 

A dense network of recording rain gages with a reasonably long period of 
record is required to determine the intensity-duration-frequency relationship — 
of areal rainfall to that of a centrally located point observation. The Weather — 
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4 
Bureau has published the results of an investigation of this problem in : 
Technical Paper No. 29.(4) Seven networks having the most suitable density 
of gages and length of record were selected from all available networks of 
this type. Even on this selective basis, the length of record available averaged 
only about 12 years. The results were quite scattered but average relation- 
ships were drawn for durations of 1/2; 1- and 24-hours with intermediate 
values interpolated. 

A study by the State of Ilinois(5) of one of their dense networks of rain 
gages indicated a considerably smaller reduction than the Weather Bureau 
for areas up to 100 square miles and a duration of 6 hours. Part of the dif- 
ference may be due to differences in the methods of analyses. 

Because of the lack of comprehensive information on this matter, it is 
considered that the reduction to point rainfall values for size of drainage area 
should be estimated conservatively. On this basis, it is proposed that use be 
made of average relationships between point rainfall at storm centers and the 
maximum average depths over increasingly larger areas which can be de- 
termined for a given region from storm study data suchas that published by 
the Corps of Engineers.(6) In general, two relationships between maximum 
point rainfall and average depth over surrounding areas are apparent, one 
for intense thunderstorm rainfall of short duration and limited extent 
(summer type) and one for widespread general storms (winter type). The 
summer type rainfall can be considered to control areal variation for dura- 
tions less than four hours with the winter type rainfall controlling for 
durations greater than six hours. 

If it is assumed that all storms are widespread enough to cover areas 
equal to or larger than an area of interest, a rational adjustment to the re- 
lationships based on storm center rainfall.can be evolved. On the basis of 
a random distribution of storm centers over a circular shaped drainage 
basin, it can be reasoned that the average reduction factor to a fixed-point 
record typical of a particular drainage basin should be about one-third that 
based on storm center rainfall. The following results are obtained: 


Reduction Factor to Point Rainfall, % 


Size of Area, Square Miles Summer Type Winter Type 
a 100 100 
10 96 109 
100 90 98 
1,000 83 % 
10,090 -- 92 


It is recognized that small summer thunderstorms often cover only a few 
square miles and thus affect the intensity-duration-frequency relationship at 
point locations with little or no effect on moderate sized drainage basins. P| 
Disregard of this effect, which would tend to decrease the above reduction ? 
factors, is considered to be conservative, but appropriate in view of the ; 
present state of knowledge of the matter. a 

In a region where topographic or other factors produce a variation in ‘ 
annual rainfall, recognition must be given to these differences in selection 
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or adjustment of the point rainfall record to be used. Average annual rainfall 
serves as a Satisfactory basis for adjustment within a region of similar rain- 
fall characteristics. 


Rainfall-Runoff-Relationship 


The average rainfall-runoff relationship shown on Fig. 3 of Appendix I for 
the vicinity of Washington, D. C., is based on rainfall-runoff data for natural 
drainage basins in that area. The long time average annual rainfall at 
Washington is about 41.5 inches. According to records of the U. S. Geological 
Survey the average annual runoff of gaged drainage basins in the vicinity of 
Washington varies from about 12 to 18 inches per year. In the rare cases 
where measurements of runoff are not available for drainage basins in the 
vicinity of or typical of the point of interest, estimates of average annual 
runoff can be obtained from climatological maps or estimated from climato- 
logical factors. On the basis of a simple relationship, developed by the author, 
giving the annual loss (precipitation minus runoff) for humid areas as 


Loss = 0.6 (Average temperature - 11) (3) 
Loss = 0.6 (56.8 - 11) = 27.5 inches, 


the annual runoff in the vicinity of Washington, D. C. would be 14 inches. 
Eq. 3 is based on records of about 30 drainage basins scattered over the 
humid section of the United States with an average error of 1.8 inches of 
runoff. 

The desired relationship between rainfall and direct runoff is one for 
average conditions on the assumption that storms of any given frequency will 
occur on various conditions of runoff potential, and over a period of years 
the average result would be that for average runoff conditions. One alterna- 
tive would be to develop rainfall-runoff relations in terms of rainfall duration, 
amount, season of the year and antecedent conditions. Introduction of the 
frequency factor would require analysis of a large amount of data and a 
number of simplifying assumptions. Application of the results of such an 
analysis would require a reverse application of the assumptions and selection 

- of specific values of the variables for use in synthesizing a flood frequency 
relationship. A practical alternative for natural basins would be to compute 
the actual runoff for the two or three probable maximum flood producing 
storms each year for a number of years under the actual associated soil con- 
ditions and develop an annual flood series. It is believed that the difficulties 
associated with these alternatives would leave the resultant flood frequency 
relation open to more uncertainties than the procedure proposed herein. 

The rainfall and runoff data were plotted on an average monthly basis and 
a mean curve drawn to give slightly higher runoff than the average, starting 
with an initial loss of 0.25 inch. Although it would at first appear that the 
relationship for isolated storms should give considerably more direct runoff 
than average monthly data, it must be remembered that in some drainage 
basins as much as half of the total annual runoff is ground-water discharge, 
whereas during isolated flood events the percentage of the peak discharge 
contributed by ground water is much smaller. The upper end of the curve 
was controlled by assuming all losses (2.75") satisfied (except for a continu- 
ous contribution of 5 %to ground water) with a total rainfall of 6.25 inches. 
These values were based on experience gained from detailed rainfall-runoff 
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investigations of similar areas. The average rainfall-runoff curve was then 
converted to a per cent runoff curve for more convenient use. 

The curves of Fig. 3, Appendix I, for adjusting the percentage of runoff on 
natural areas to obtain the percentage of runoff on a developed area, taking 
into account the portion of the area that is impervious were computed with 
the equation: 


Adj. RO % =% Impervious Area +(100 - %Imp. Area) %Nat. R.O. (4) 
100 


The equation assumes 100 per cent runoff from the impervious area. Wet- 
ting of the surface and ponding of water in small depressions and behind 
various surface irregularities can easily account for from 0.05 to 0.10 inch 
of water even on a fairly smooth airfield pavement. However, assuming 100 
per cent runoff provides a small factor of safety and is the only such un- 
evaluated one in the procedure described herein. Furthermore, in many in- 
stances the depression storage may be filled by antecedent rainfall. 


Effect of Storage and Rainfall Duration 


Fig. 1 and the following discussion illustrate the conditions and assump- 
tions under which the so-called “rational method”, when used with a proper 
runoff coefficient, results in reliable estimates of discharge, even though the 
storage effects of the particular conduits involved are not evaluated 
separately. 

The isosceles triangle of Fig. 1 with a time base of T, is assumed as the 
area-shape diagram of a particular area. The diagram represents the dis- 
tribution of the aréa with respect to distance in time (or length) away from 
the point of interest. If 1 inch of water were deposited instantaneously or 
within a very short time over the whole area and the water moved at a con- 
stant velocity of L/ T, down the drainage paths to the point of interest, the 
hydrograph observed “would be the equivalent of the area-shape (Gischares) 
diagram (1 acre-inch = 1 cfs-hour) under the “rational method” concept of 
time of concentration. For many typical drainage basins an isosceles triangle 
can be accepted as a reasonably close approximation of the actual area-shape 
diagram. Assumption of a rectangular area-shape would change the value of 
q,, in Eq. 8 by only 15%. 

The use of an instantaneous rainfall simplifies the following analysis of 
the assumptions inherent in the “rational method”. However, if a finite rain- 
fall duration t were used for the hypothetical drainage area of Fig. 1, the 
area-shape diagram would become essentailly a trapezoid with top width t, 
and bottom width t + T,. So long as t is made about one- -seventh or smaller 
part of T, the results of this analysis would be affected very little. 

Water *cannot flow without first building up depth, and with depth comes a 
storage effect (also variable velocity) on discharge. Accordingly, observed 
hydrographs of discharge for unit rainfalls are not area-shape (discharge) 
diagrams but rather are area-shape (discharge) diagrams modified by the 
storage and velocity or wave movement characteristics of the collecting 
media (2). 

The effect of drainage channel storage on discharge from headwater areas 
has been found to be closely approximated by reservoir type storage 1) of an 


1808-12 HY 5 October, 1958 


amount 
rte AS 
Aa (5) 
where T, is an index of storage effect with the dimension of time. T, can be 
evaluated from the recession side of discharge hydrographs after inflow has 
essentially ceased by the equation, T, = -t/logeC;) where C, is a recession 


a 


§ 


coefficient, Qo/Q1; for two values of discharge separated by time, t. Because } 


direct runoff recessions can only be estimated by separation from total flow 
for discharge hydrographs of natural drainage basins, values of T, so de- 
termined are subject to a number of inaccuracies as well as the possibility 
of variation with season of the year and size of flood. 

Referring again to Fig. 1 and assuming a unit area of 1 acre with 1 inch 
of rain deposited thereon and time in hours, 


& 


Vol. 1efs-hr (645 cfs-hrs for 1" on 1 sq. mile (6) 
SS aa A So T é 


and,if.T. = T,/2, 
vd q a; /2 
pry (7) 


in which q; is the peak of the isosceles triangular area-shape (discharge) 


diagram in cfs-hours per acre and q, is the peak of the discharge hydrograph ~ 


per unit of area (per acre) per unit of runoff obtained by routing the area- 
shape diagram through a reservoir type storage of T, = T,/ 2 hours. 
Then from Eq. 6 


q, = Si. = Vole = 2b (8) 
P D iy * 
and 
tn (1.01 for 1 acre or 645 for 1 sq. mi.) (9) 


Under the assumption of one inch of rain in a short period the outflow 
hydrograph with peak q,, is an instantaneous unit hydrograph per unit of area 
for the basin with ordinates in cfs per acre. A unit hydrograph for a given 
area may be defined as the discharge hydrograph of direct runoff resulting 
from 1 inch of runoff from a rain occurring in a specified unit of time. The 
unit hydrograph procedure is now the most commonly used procedure for 
synthesizing discharge hydrographs. 

The peak discharge, Qp; for a drainage basin after the unit hydrograph 
procedure using the value of Gp from Eq. 8 is 


Q, = KgA 1 T. Ty = KIA, (10) 
To 

where A equals the area in acres, I, is the average rate of runoff over the 
time of concentration, T,, in inches per hour or the equivalent of Ci in the 
rational approach, Qp = CiA, and Kg is a reduction factor for adjusting the 
result, in that the rainfall and resulting runoff, I,, is the mean rate over the 
time of concentration, T,, rather than over the very short period of time 
assumed in developing the unit hydrograph peak, Ap: 

Thus it has been shown that, with the assumption of an isosceles area- 
shape diagram and a basin storage factor of Ty, = T, /2, that the rational ap- 
proach gives the same result as use of a unit hydrograph except for the 
factor Kg. With the effective rainfall spread over a duration of T, hours in 


: 
t 
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the rational method, the value of Kq is 0.87 for a typical unit hydrograph 
shape and rainfall distribution of duration, T¢. 

The value of Kg was evaluated as 0.87 by applying a typical unit hydrograph 
shape to a typical rainfall distribution over the time of concentration and 
comparing the resulting peak discharge with that which was obtained by as- 
suming the same total effective rainfall in one period. The rainfall distribu- 
tion in per cent of total for a time of concentration divided into seven equal 
periods was: 4, 8, 10, 20, 40, 12 and 6%. These numerical values are 
typical for durations or times of concentration of from one to six hours and 
the arrangement is critical for unit hydrograph usage. For longer durations 
the shape would be somewhat more concentrated. 

However, area-shape diagrams for actual drainage basins vary somewhat 
from isosceles triangles. The variation can be measured by the ratio Lea/L 
(equals 0.5 for an isosceles triangle) where ‘Leg is the distance along the 
principal channel to a point opposite the center of area of the basin. Also, 
although T, = T, /2 is fairly representative of the storage effect of many 
drainage channels, the value of T, does vary from basin to basin and from 
conduit to conduit. Accordingly, equations (8) and (9) need adjustment in 
cases where the actual values of T, and Lea /z, Vary considerably from 
0.5 Tc and 0.5, respectively. This can be accomplished as follows: 

ao tee) 2 fe (11) 
in which 
K, = 2/3 (1 + Lea/L) 
qT; /e + 0.5 


c 


(12) 


Eq. 12 was developed by routing additional triangular area-shape diagrams 
of other than isosceles character (variable values of Leg L) with various 
values of T./ T, and evaluating the resulting effect of the variables on qpnTo. 
Although more approximate in the case of area-shapes differing from the. 
triangular shape assumed in developing the relationship, Eq. 12 serves asa 
practical adjustment for other typical shapes. 

For the Washington, D. C. area a value of K, = 0.86 is appropriate for the 
natural areas. The resulting value of 0.89 for q T, is also applicable to 
sewered areas and to paved surfaces with apparently little limitation as to 
location. For non-paved overland flow surfaces a value of Apt of 0.50 is 
Satisfactory. Values of QpT in Table 1 would be about 570 cfs-hours per 
Square mile and 53 cfs-min/per acre for a value of K, of 0.86 and about 320 
cfs-hours per square mile and 30 cfs-min/per acre for a value of K, of 0.49. 

The adjustment factor, K,, of 0.86 for storage effects and basin shapes, 
when multiplied by 1.03 and the adjustment factor of 0.87 (Kg) for an effec- 
tive rainfall duration of T,, gives a value of 0.77 (0.44 for overland flow on 
turf) as the total adjustment factor, K, to be used in 


Q, = Kal, (13) 
where 
K = 1.03 Ky K. 


Note: 1.03 (658 when area is in square miles) is the actual value of GpTc 
for an isosceles shape and Bi oi = 0.5. 
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APPLICATION 


Appendix I outlines the step by step procedures to be followed in com- 
puting flood frequencies by the procedures described herein with the variable 
factors for natural drainage basins evaluated for use in the vicinity of 
Washington, D. C. Table 2 gives sample computations for Goose Creek near 
Leesburg, Va., Occoquan Creek near Occoquan, Va., Accotink Creek near 


Annandale, Va., and Little Falls Branch, Bethesda, Md. The observed values ; 


of flood probabilities given for comparison are based on U. S. Geological 
Survey data with 26 years of record for Goose Creek, 25 years for Occoquan 
Creek, 10 years for Accotink Creek and 12 years for Little Falls Branch. It 
should be noted that the synthetic frequencies constitute a partial duration 
series and as such the values for recurrence intervals of 10 years and less 
should be somewhat larger than the annual series values with which they are 
compared. 

In cases where it is desired to compute the discharge at a large number 
of locations for a selected frequency of occurrence such as in the design of 
storm sewer systems, working curves to expedite the computations can be 
developed on the basis of the procedures described herein. Using Eq. 13 
values of KI,, K = 0.77 for the vicinity of Washington, D. C., for a specific 
frequency are plotted against indicated values of T, with a family of curves 
covering the desired range of imperviousness. Since a K value of 0.77 is 
adequate for natural areas, overland flow on pavements, and sewered areas, 
one set of curves such as Fig. 4 of Appendix I is all that is needed for each 
frequency of interest in most cases. If conditions of overland flow on turf 
are involved, values of KI, read from the zero or other appropriate impervi- 
ousness curve, reduced = 0.57 (9.44/0.77), are applicable. If overland flow 
on turf is of primary interest, separate KI, curves for this condition can be 
prepared. 

The initial collection or inlet time at the head of each line can be com- 
puted by use of appropriate factor (C; = 0.30) or estimated on a basis of 
judgment (ordinarily about 10 minutes). Because the value of C; for over- 
land flow is twice that for sewered areas, the computation of inlet time is 
kept separate from the computation of T, for the sewered portion of the 
collection system, and the inlet time is added each time to the computed 
value of Te for downstream sections of the line. After a point is reached 
where the time of concentration is greater than about 2 hours, retention of 
the inlet time as an additive factor can be stopped as a practical measure. 

In cases where storage is available along a drainage system such as at 
an inlet for airfield drainage or at a pumping plant, the peak discharge as 
obtained by the procedures of this paper can be adjusted for the effect of 
the storage available. This is accomplished by adding to the computed time 
of concentration, T,, for the location, 1.8 times the value of T, determined 
by successive trials as the average over the range of outflow and available 
storage. This is accomplished by assuming that the peak discharge for the 
location will be reduced by the ratio of the available storage to the volume 
of the runoff. The available storage is then divided by the trial regulated 
discharge so determined to get a trial value of Ty. The regulated peak dis- 
charge is then computed for the location by the regular procedure using an 
adjusted time of concentration equal to the original T, + 1.8T,. If the 


resultant peak discharge is not sufficiently close to the assumed trial value, 
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TABLE 2A COMPUTATION OF SYNTHETIC FLOOD FREQUENCIES 


OCCOQUAN CREEK NEAR OCCOQUAN, VA. 


Natural basin but with several reservoirs which reduce flood peaks slightly. 


Dr. A. = 546 sq. miles; L = 0 miles; s = 5.6 /miles = 0.11%; n = 0.06 
Annual basin rainfall is = Washington rainfall; Adj. for size = 0.96 


a (10 x hO x 206) 036 


= 22.) hours 


Adj. for effect of reservoirs + 1.6 


Adje Te — 
Qp = 500 AT, = 273,000 Ip 


Freq. (Recur. Int., Yrs.) 
Point Rainfall, inches 
Basin Rainfall, inches 
Percent RO, Natural 
Runoff, inches 

Ip, inches per hr. 

Q%, cfs 

Obs. 25 Yr. Record /1 
Obs. 25 Yr. Record /2 


ACCOTINK CREEK NEAR ANNANDALE, VA. 


2) hours 
2 5 10 
3.42 Ig 3s) 5.05 
3.28 4.16 4.85 
37 43 47.5 
Dede ato) abe ee 
0050 07h 096 
13,800 20,300 26,200 
13,500 20,000 2h,),00 
12,700 20,000 25,700 


30, 000 
3h, ,000 


Dr. A. = 23.6 sq. miles; L = 11.0; s = 0.293 n = 0.05 


ie .7 ES OP 327K L016 8 Hours 


. Assume 5% Average Imperviousness, 197-1956 


Areal Adje less than 1% - 
Qp = 500 AIr = 11,800 I, 


Freq. (Recur. Int., Yrs.) 
Point Rainfall, inches 
Percent RO, Natural 
Percent RO, Adjusted 
Runoff, inches 

Iy, inches per hr. 

Qn, cfs 
obs. 10 Year Record fi 


bs 


Ignore 

2 5 10 

2.91 3653 4.06 

Bhs 38.5 (2.3 

373 41.6 4522 
1.08 1.47 1.83 
0.160 2216 0279 
1880 2550 3180 
1550 2530 3200 


U.S. Geological Survey records, after Gumbel 
Best fit, with recurrence interval = (Ney 


25 
47h 
46.5 
49.2 
2.33 
343 
050 
),000 


50 100 
6.72 7240 
665 helo 

57 60 
3057 4.26 
0153 alg 

41,800 8,500 
34,000 (38,000) 
1,000: (48,500) 

50 100 
5e31 5.87 
50.2 53.3 
52.7 56.0 
2.30 3629 
2412 Bly 

860 


hy 5700 
(4600) (5220) 
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TABLE 2B COMPUTATION OF SYNTHETIC FLOOD FREQUENCIES 
LITTLE FALLS BRANCH NEAR BETHESDA, MD. 


Principal channels mostly open but with numerous storm sewers. 

Dre Ae = hell sq. mio3 L = 2.75 mies s = 53'/mile = 1.0%; n = 0.05 
As of 1950 about 22% of area was impervious. 

10% natural channels eliminated; 30% of area storm sewered. 

Cy = 1e7 = O02 (1.7 = 0042) = Lele 

Te = leah (20 2 2e7o x 20s) ean 1.65 hours 


Qp = 500 AT, = 2055 I, 


Freq. (Recur. Int., Yrs.) 2 5 10 25 50 100 
Rainfall, inches 2.15 2.51 2.82 3.27 3.65 4.02 
Percent RO, Natural 28 31 33.5 37 3925 )2 
Percent RO, Adjusted hy 6 7.5 51 53 55 
Runoff, inches 0.95 1.15 1.34 1.67 1.94 2.20 
ire ae per hr. Os57k is : O70: . 05833. o1-OLL eT e eee oe 
1180 14h0 1670 2080 210 2760 
Ee o Yr. Record /1 10,0 shO 1710 2060 (2300) (2570) 


GOOSE CREEK NEAR LEESBURG, VA. 


Dr. A. = 338 sq. miles; L = 33 miles; s = 11.4'/mile = 0.22%; n = 0.06 


= 10 x 33 x 206 0.6 
nl a top ae ) = 16.2 hours 


Annual Basin rainfall is <= Washington rainfall; Adj. for size = 0.96 
Q = 500 AI, = 169,000 Iy 


Freq. (Recur. Int., Yrs.) 2 5 10 25 50 100 
Point Rainfall, inches 3.27 HOO bitSasaStS0eeG.se 6.95 
Basin Rainfall, inches 312 3.92 h.56 5.38 6.05 6.67 
Percent RO, Natural 3505 hi .y 45.5 51 55 58.5 
Runoff, inches Toe 1.68 2.07 2.75 3.33 3.90 
I,, inches per hr. 0.068 2103 2128 Pi ty [0) 2206 22h 
I, efs 11,500 17,500 21,600 28,700 3,800 0,700 

Obs. 26 Yr. Record /1 10,000 20,400 27,300 36,000 i,23500 (49,000) 


{1 U.S. Geological Survey records, after Gumbel 
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the computation is repeated using the computed regulated peak discharge to 
get a new trial value of T,. A similar adjustment was made in Table 2 for 
Occoquan Creek to adjust for the effect of existing reservoirs. 

Values for the coefficient, C,, and adjustment factor, K,, presented here- 
in for overland flow areas and areas with storm sewer drainage, are con- 
sidered to have fairly wide application. The range of slope, s, for the drain- 
age areas listed in Table I is perhaps not great enough to fully support such 
an assumption. However, it is believed that some extrapolation with respect 
to this characteristic is permissible, particularly in view of the fact that 
slope is effective to the 0.3 power only. 

The time and storage factors for the overland flow areas of Table I are 
based on hydrographs of runoff which were adjusted to eliminate the effect 
of storage due to ponding, channel flow to inlets, or flow in collecting 
pipes.(8) Adjustment for appreciable channelization or ponding should be 
made when synthesizing discharge values with the procedures of this paper. 
The time and storage factors were also found to vary with the rate of rainfall, 
especially on the paved areas. The values of these factors reported in Table 
I are averages for storms with rates of runoff varying from 0.5 to 3.0 inches 
per hour for the turfed areas and were selected for supply rates of about 6 
inches per hour for the smaller paved areas. 

The time and storage factors for the storm sewered drainage basins in 
the Louisville, Ky., area as listed in Table I, are based on unit hydrographs 
developed from observed hydrographs of runoff for isolated summer storms 
of short duration.(9) The full capacity of the lower sections of each of the 
systems was approximated in one or another of the storms analyzed, except 
for the Southwestern Outfall. None of the storms caused surcharge of the 
main trunk sewers. Main trunk shapes included were arch, egg-shaped, 
horseshoe, semi-elliptical and circular, but other factors obscured the re- 
lationship, if any, existing between main trunk shape and storage effect. The 
times of concentration, T,, as used herein and obtained from hydrograph 
analyses, were smaller than those developed by the rational approach for 
the sewer systems, especially for the larger systems. 

The value of 0.87 for Kq is based on an average storm distribution and 
typical unit-hydrograph shape. Rather extreme variations from the assumed 
typical conditions would be necessary to change this factor by as much as 
10%. 

Subject to the comments in the preceding paragraphs of this section and 
to subsequent comments with respect to rainfall-runoff relationships, it is 
believed that the procedures and coefficients presented herein for storm 
sewer systems and overland flow are satisfactory for general application 
and will give reasonable results in most locations. ’ 

A number of adjustments and selection or verification of coefficients 
are necessary in application of these procedures to natural drainage basins 
with unusual characteristics or ones generally different than those of the 
Appalachian Highlands. Suggestions for developing these adjustments are 
presented in the following paragraphs. 

The values of Cy and K, to be used can be determined by using Eqs. 1 
and 11, if a tabulation of unit hydrograph data such as that published by the 
Corps of Engineers is available for streams in the vicinity. If such data 
are not available, they can be developed from streamflow records or a time 
of concentration estimated by conventional procedures and a storage factor, 
Kg, selected from general knowledge of runoff characteristics of the region. 


; 
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As previously stated the value of Kg of 0.87 would need little or no adjust- 
ment. In any event values for length, slope and friction can be measured or 
estimated with a satisfactory degree of accuracy for basins with streamflow 
records as well as for the area for which synthetic flood frequencies are 
desired. 

On some drainage basins there is fairly conclusive evidence that the unit 
hydrograph characteristics vary with the magnitude of the particular flood. 
Where such variation can be defined, a more accurate flood frequency would 
be obtained by varying the values of C; and K, on the basis of the recurrence 
intervals desired. 

The rainfall intensity-duration-frequency data of T.P. No. 25(3) provide 
a basic source of rainfall data for general use. In a number of instances, 
the author has obtained satisfactory results by using data from T.P. No. 25 
for the station nearest or most appropriate for application to a particular 
drainage basin. For durations less than 4 hours little adjustment’ is neces- 
sary while for durations greater than 6 hours the station values can be ad- 
justed by the ratio of the average annual rainfall for the basin to the average 
annual rainfall for the station. There is some indication that this results in 
a slight over-adjustment. 

The rainfall-runoff relationship for natural areas of Fig. 3, Appendix I, 
has been used with fairly good results on drainage basins of the Appalachian 
Highlands without adjustment. It is also considered that the relationship 
would provide a satisfactory basis for a storm sewer design in many urban 
areas. However, the rainfall-runoff relationship is the factor that would 
require the most effort in developing adjustments for use on natural drainage 
basins in areas significantly different than those in the vicinity of 
Washington, D. C. 

One basis for adjustment of the relationship is to compare the average 
annual runoff-rainfall ratio for the area of interest with that for the vicinity 
of Washington, D. C., through the medium of the average monthly rainfall 
for the respective areas. The relationship of Fig. 3, Appendix I, was drawn 
slightly to the right of the average annual runoff-rainfall ratio of 33.7% 
(14.0,41.5) plotted against an average monthly rainfall of 3.46" (41.5/12). 
With a selected initial loss for the area of interest (0.25” for vicinity of 
Washington, D. C.) as one point, a relationship can be drawn on Fig. 3, 
Appendix I, relative to the average monthly data for the area of interest. 

Another procedure for developing a rainfall-runoff relationship for use 
in a specific area would be to utilize available flood frequency data for 
streams typical of the region. The computations as illustrated in Table 2 
could be made in reverse order starting with the observed flood frequencies 
to obtain the natural runoff percentages for plotting against the appropriate 
rainfall values. 
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APPENDIX I - PROCEDURE 


1. Determine ie 


Te = Cy Lt 06 Lut =10In 
an oo | 
Solve equation or use Fig. 1 


L = length of principal drainage channel ; 
s = weighted slope of principal drainage channel. t 

= mean height of channel profile divided by L/2 
n = friction factor 


Average Natural Channels, overall 0.05 
Paved Channels 0.014 to 0.040 
Overland Flow and Improved Channels 
Smooth paving 0.016 
Bare packed soil 0.10 
Poor grass cover 0.30 
Average grass cover 0.40 
Dense grass cover 0.50 
Storm Sewers 0.012 to 0.025 


ee a a 
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(Add time of overland flow (inlet time) for small areas where it is a 
significant factor) 


Values of C; Hours/mile Min./foot 
Natural Drainage Basins uy 0.60 
Overland Flow 0.85 0.30 
Sewered Areas (Add inlet time) 0.42 0.15 
Partially Sewered Area C; Adj. = 1.7 - % Sewered + %Nat. Chan, Elim.(1.70-0.42) 
2 x 100 


2. Determine total point precipitation for time T, and frequencies of interest 
from Fig. 2. Adjust for size of drainage area if necessary. 


3. Determine percent.of runoff: 
a. For natural area. 
b. Adjust for percent of area that is impervious. 
Normal Range 
Residential, 10 to 60% 
Industrial, 40 to 70% 
Commercial, 60 to 100 % 


4. Compute total direct runoff in inches and divide by T, to obtain I... 


5. Compute Qp- 


Natural Areas = 0.77 AI, Area in acres 
Overland Flow, pavement} = 500 AI, Area in sq. miles 
Sewered Areas Gp Te = 970 
Overland Flow, Turf = 0.44 AI, Area in acres 
= 280 AI, Area in sq. miles 
dp Te = 320 


6. For Storm Drainage Design 


Sample Tabulation 


Column 

Number Heading Use 
1 Location Identification of location and section 
2 Line of drain. 
3 Invert Elev., feet Basic data for computation of slope, etc. 
4 Le feet ? w " " ’ Ua " 
5 Slope, s, in Basic data and computation of equiva- 
6 Friction Factor, n lent length, leading to evaluation of T,. 
7 Equivalent Length, L' " 
8 Summation L' " 
9 


Tc plus Inlet Time addition of Inlet Time. 
Ti Incremental Area Tabulation of data to arrive at a value 


12 Incremental Area of weighted Imperviousness. Can be 
Imperviousness eliminated if approximation is per- 
13 Incremental Area x Im- missable. 
perviousness " 
14 Accumulated Area x Im- " 


perviousness 


sf 
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15 Accumulated Area Approximation is permissable 4 
16 Weighted Imperviousness Column 14/Column 15 ; 
17 0.77 I, Obtained from curves such as Fig. 4 


of Appendix I for appropriate frequency 
and imperviousness. 

18 Discharge, Qp» in cfs. Design discharge equals column 15 x 
column 17. 

19 Conduit List shape, size and type of conduit. 


¢ 


APPENDIX II - NOTATION 


A Area in acres or square miles 

Ct Coefficient, hours/mile or minutes per foot, in T, =  L 0.6 

I, Average rate of direct runoff for time Tg, in inches per hour 

Kg Reduction factor to be applied to instantaneous unit hydrograph peak to 
adjust for effect of an effective rainfall duration equal to T, 

Kg Adjustment factor to be applied to standard value of Gp T, = 1.03 for ; 
values of T</T, and Lea/1, other than 0.5 - 

K Combined adjustment factor equal to 1.03 Kg K, or 658 Kg Ks { 

Ls Length of principal drainage channel, feet or miles : 

L' Length of equivalent channel with standard slope of 1% and friction factor : 
of 0.10 

n_ Friction factor, generally that used in Manning equation ’ 

qj Peak of area-shape (discharge) diagram for one inch of runoff per unit 

area (1 acre or 1 square mile) of drainage basin 

Peak of instantaneous unit hydrograph (per acre or Square mile) obtained 

by routing area-shape (discharge) diagram through reservoir type : 

storage, T, ; 

Qp Peak discharge of runoff hydrograph for total basin and total direct run- — 
off volume spread over T, hours 

AQ Small increment of discharge 

AS Corresponding increment of storage 

S Weighted slope of principal channel in %, feet per hundred feet, and equal 
to the mean height of the channel divided by one half the length 

Tg Time of concentration in hours or minutes and equal to time from end of — 
rainfall to inflection point on recession side of hydrograph { 

Ts Index of reservoir type storage and equal to AS/AQ. 

M_ Serial number of flood from arrangement of annual floods in order of J 
decreasing magnitude 

N_ Total number of years of record 
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The history of the basic-data programs are reviewed. The divining rod as 
a mystical hydrologic decision-maker is contrasted with modern basic-data 
programs. The extent and kinds of gaps in basic data are discussed with 
major emphasis on the deficiencies in research and interpretation of data. 


Water problems, such as recurring floods and droughts, water surpluses 
and shortages, are everyday news. To solve these problems and to provide 
information for the sound and intelligent development of water resources, 
basic-data programs that are large in extent and scope have been established. 
How well are these programs doing their job? 

Basic-data programs, like so many other things of the 20th century, are 
relatively new to the American scene. Much of the history of man reveals that 
lack of information has led him into enterprises that proved fruitful and that 
it has also wasted much effort. Although Columbus was careful to marshal 
all facts available, he was encouraged to try his westward trip by two errors, 
an underestimate of the size of the globe and an overestimate of the eastward 
extension of Asia. Ignorance of a different sort led those who followed him. 
Indeed, the North American continent was opened by conquistadores and colo- 
nists lured by fancies, not facts, and dreams, not data, even though increasing- 
ly accurate reports of what lay along the westward trails became available to 
those who would read them. In our West, fur traders, explorers, and later, 
scientists, reported what they saw. Yet ideas about the West remained fanci- 
ful, partly because the area was indeed fabulous, but also because the mind of 
the Nation was receiptive to the romantic idea. The mistaken hope that “the 
rain follows the plough” was a much more powerful influence behind the settle- 
ment of the subhumid Great Plains, than J. W. Powell’s reports on rainfall 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1809 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 5, October, 1958. 
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published as early as the 1870’s. Powell warned of the hazards of agriculture ~ 
in areas that are neither humid nor arid and forecast the boom and bust that — 
have marked agriculture on the subhumid-semiarid Plains. 

The year 1879 is a key date in our progress. It was then that John F. 
Stevens, a railroad surveyor, discovered Marias Pass, in Montana, the last 
unknown mountain pass across the Rockies, thus, in a sense, marking the close 
of an era of geographic exploration in the United States. The year 1879 also 
saw the founding of the Geological Survey and the soon-to-follow hydrologic ~ 
investigations as we know them today. With the waning of the era of explo- 
ration, basic data acquired a new meaning. The country now had to develop 
within fixed boundaries. In an earlier stage, action, whether based on fact or 
fancy, was the important thing; now, in a world of decreasing choices, de- 
cisions must be more sharply drawn. 

In the meantime, the hydrologic frontier has been passed. Cheap water is 
no longer abundant everywhere. Projects for the control and use of water are 
now problems of national and even international concern. Sooner or laterin ~ 
planning a water project, the point is reached where a decision must be made: 
“drill here”, “make it so many acre-feet”, or “install so many kilowatts” ac- ? 
cording to the kind of project. In our planning we can use arbitrary methods 
and thereby short-circuit all needs for basic facts, we can resort to the oppo- i] 
site extreme of massive data gathering so as to provide minute details for ¢ 
every problem, or we can balance a moderate collection of facts with research 
and interpretation. 

Man’s early resort to magic, however, is still with us. Of all the hydro- 
logic artifices, the one most commonly used is a forked stick of willow, ash, 
peach, or witch hazel. Such a forked stick, where used in the search for 
water, is called a divining rod or “witching stick”. It is more widely known 
and used than the hydraulic engineer’s current meter. The divining rod is 
believed by many to be capable of indicating the presence of water under 
ground. The user grips a branch of the fork in each hand and walks about with 
the point out in front. When he passes over an “underground river” or a ; 
“ground-water vein”, the point is suddenly and inexplicably jerked down, point- 
ing at a spot on the ground. “Here”, he tells you, “one will find water.” Why — 
does this legacy of magic lore from the dark corners of the past persist along- 
side and despite modern technology? Or can we be scientific about witchery? 

Witching for water is prevalent throughout the world because, rightly or ? 
wrongly, it gives an answer to many who seek water but who have not yet been 
reached by the hydrologic-data programs. Less than 10 per cent of our own ‘ 
country has yet been mapped or thoroughly explored for ground water; even ’ 
where such work has been done, the reports and maps seldom reach all those 
who need them. Water in the ground is out of sight; therefore, to find water I 
at any particular spot where the farmer may wish to dig a well may appeara 
highly uncertain business to the uninitiated. Like the medicine man, the watems 
dowser gives a direct and unqualified answer depending only on the dip of a : 
stock. This is an old trade, with many devotees. The literature dealing with 4 
it is large and ever growing (Ellis, 1917; Roberts, 1951; Sowder, 1955.) 

The dowser gives a specific answer to allay the farmer’s anxiety about ns 
water supply and to enable him to keep his mind on the hard job of farming. 
Vogt (1953) suggests that water dowsing persists because it offers a short cut 
to a decision about a difficult feature in an uncertain environment. In many 
places random drilling has a good chance for success if only a domestic sup- 
ply is needed. The geologist’s advice to drill on a 50-50 chance might be 
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“ejected whereas on the authority of a forked stick much useful ground water 
s developed. 

Modern hydrology offers no support of water witchery. A lack of principles 
whereby experience can add to knowledge condemns divining as an acceptable 
lydrologic tool. The dowser’s efforts to excuse “dry” holes supply a true 
neasure of the trade. It is here that the dowser demonstrates his profession- 
il achievement. The driller is generally the chief target. Here are some of 
he excuses: “The driller busted the structure of the vein,” “The driller 
noved the stake an inch,” “The hole was drilled crookedly,” or “My knife 
nust have deflected the rod on that day.” This kind of fare must compete 
vith experimental science. A drilled hole, wet or dry, adds to the body of 
cnowledge of the geologic structure and sharpens the ability of the geologist 
o make other recommendations in the region. The important job of the 
rround-water hydrologist a generation or more ago was to find water. The = 
ask today, as the late O. E. Meinzer put it, is much greater; it includes also 
he determination of availability, quantity, and quality of ground water. The 
nost severe indictment of water dowsing is its harmful effect on efforts to ob- 
ain dependable information about the occurrence of ground water in an area. 

Divining rods represent a class of techniques for reaching decisions pain- 
essly but with little regard for factual substance and none for proven scien- 
ific principles. By way of contrast, we shall explore an opposite extreme 
-epresented by the idea that proper answers will become directly evident if 
ynly we can collect enough data. 

The medieval appeal to witchery and dogma led Sir Francis Bacon in the 
ate 1500’s to reason that truth could be determined only through observation 
yf nature. It was his further idea, in part, that if one only makes a sufficient 
umber of observations, the truth will thereby become evident. This philoso- 
Ihy was a useful antidote to the divining-rod method of decision, but we have 
earned that, carried to extreme, data collection can be equally sterile. In the 
livining rod’s use there is lack of both data and principles; data collection can 
ead to a surfeit of information to the neglect of principles. Both situations 
lave historical and present interest for us. 

Ever since Bacon’s writings about the scientific method, scientists have 
soncerned themselves about the relations between observation and conclusions. 
Ine of the first to analyze this problem in the modern sense was Augustus 
le Morgan, who, in 1872, write, “Modern discoveries have not been made by 
arge collections of fact, with subsequent discussion, separation, and resulting 
leduction of a truth thus rendered perceptible.” At the Mid-Century Confer- 
nce on Resources for the Future, Neff (1953) stated that “We are pouring too 
nuch of our resources into the mass collection of data, without devoting suf- 
icient resources to the development of hypothesis essential to the collection 
ind the use of data presently available.” Neff and De Morgan both claimed in 
sssence that a special objective of data collection is the test of scientific 
lypotheses. By way of criticism, De Morgan cited the Greenwich Observatory, 
vyhich was established in 1675 “to observe, observe, observe away at the 
noon, until her motions were known sufficiently well to render her useful in 
uiding the seaman.” The discovery of the general law of gravitation, which 
xplained not only the lunar motions but those of all other planetary bodies as 
yell, was not made at Greenwich. Newton needed only 20 or 30 observations 
O prove his generalized theory. Had it not been for Newton, De Morgan, add- 
id, the whole dynasty of Greenwich astronomers might have worked away at 
nightly observations and daily reduction without any remarkable result.” 
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De Morgan did not mean that only a genius could make significant contributions 
to general knowledge, but he used this example because the issues were clear” 
cut. His point was that principles rank ahead of data. To emphasize further 
the importance of hypotheses, De Morgan added, “Wrong hypotheses, rightly — 
worked from, have produced more useful results than unguided observation.” — 

De Morgan was especially critical of the proposed enlargement of meteoro= 
logical observations. He quoted Airy, who, in 1867, had said “Whether the ef- 
fect of this movement will be that millions of useless observations will be 
added to the millions that already exist, or whether something may be expect~ 
ed to result which will lead to a meteorological theory, I cannot hazard a con= 
jecture.” Yet we know today that many meteorologic principles and forecast — 
tools were derived and tested from these data, just as Kepler deduced his : 
three laws of planetary motion from untiring analysis of the great mass of 
data compiled during a lifetime of observations by Brahe. 

The modern view is that both the collection of data and their interpretation 
are necessary for reaching wise decisions. The problem is the proper di- 
vision of effort. If data are too few to evaluate known principles, the designer 
of engineering works must choose between delaying the works in order to col= 
lect the data or to resort to extrapolation—what is facetiously called in the 
trade “galloping transpositions.” On the other hand, data become too many 
when their collection is the main object and sufficient tools are not available © 
for putting them to use—like the watchers of the “dickey” bird who, as Aldo ~ 
Leopold (1953) put it, expend their efforts “labeling species and amassing 
facts about food habits without interpreting them.” 

Today we have extensive programs for the collection of hydrologic and re- 
lated meteorologic data. A recent report (Subcommittee on Hydrology, 1956) 
showed that every Federal agency concerned with water is a source of hydro- 
logic data, running through 30 items from consumptive use of water through ~ 
precipitation, soil moisture, ground water, and streamflow to water waste. 
Agencies such as the Corps of Engineers or the Bureau of Reclamation collect 
information mainly related to their chief activities, flood control andirri- 
gation. The Tennessee Valley Authority is active in collecting data virtually 
throughout the hydrologic cycle within its territory. The Weather Bureau and 
the Geological Survey are the only agencies that collect countrywide hydrologic 
information for public use as one of their major activities. These agencies — 
collect the bulk of the available hydrologic data—precipitation, streamflow, 
ground water, and chemical and physical properties of water. % 

The governmental arrangement for data collection works well, for coordi-— 
nation is the rule. Rarely is there duplication of effort because there is so 
much to be done that there is little occasion for one man to duplicate the work 
of another. Yet, in any particular field of hydrology there is enough talent in 
two or more different agencies to assure that each man will stay on his toes © 
lest his opposite numbers steal a march on him and make a contribution to _ 
the science in the specialty that he thought was his own. This healthy compe- 
tition, constructively directed by informal coordination of the work of all the 
agencies and by the integrity of hydrologists throughout the Government, has 
made the level of achievement high and the volume of basic data large in spit 
of inadequacies of personnel and funds. 

Nevertheless, a new water project nearly always reveals a deficiency. 
Some of these deficiencies are in the class of information. For example, the 
growth of the highway system and the upstream flood-control programs 
showed up a deficiency in information about the flow of small streams. 
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Heretofore, the emphasis in stream gaging has been on large streams, suita- 
ble mainly for power, navigation, irrigation, or water supply, neglecting the 
myriad small streams that are now of such great national concern in the 
Watershed Protection Program enacted in 1953. 

The deficiency may be in coverage—for example, information on ground 
water and on the sediment carried by streams. Ground-water investigations 
have been carried out mainly where overdevelopment has created a crisis. 
Similarly, data on sediment loads have been collected mainly in areas of con- 
cern to Federal agencies planning large dams. Broad areas, such as New 
England, the Gulf region, the Great Basin, and the Pacific Coast are virtually 
without such records. 

Most deficiencies exist because the basic-data programs are built chiefly 
around specific needs. The programs are directed mainly toward a narrow 
professional audience that is concerned with water development and control— 
chiefly those in governmental bureaus (Federal, State, and local). But hydro- 
logic decisions are made in a wider sphere. Consider the builder or pros- 
pective homeowner who erects a house on the flood plain, without any recog- 
nition of perhaps the largest bulk of hydrologic data—those on floods. Or con- 
sider the farmer who had his well “witched” and who finds it inadequate, or 
who sees his farm pond washed out or filled with silt. The major faults, 
troubles, and failures of today often are not caused by a lack of basic data but 
to their disuse. Maximum use of these data can be fostered by their interpre- 
tation and putting the results into simple guides and digests for public use. 

What proportion of the total investigative effort is devoted to the interpre- 
tation of this fund of data cannot be estimated, but it is undoubtedly less than 
10 per cent. About 20,000 pages of water data are published annually. Only a 
few hundred pages are devoted to interpretive and research results. Not that 
tools are lacking. Rainfall, for example, is recorded as the events occur. On 
the other hand, the designer of a storm sewer may need to know the largest 
rainfall he can expect in a 1-hour period once in 20 years. These facts are 
found in the original record. Tabulation, rearrangement, and fitting the 
ordered data to a statistical “model” represents efficient use and interpre- 
tation of the raw facts. Other interpretations can be made from simple charts 
such as a map of mean annual yield of rivers or from complex prognostic 
weather charts. Interpretations may be aimed at specific projects, or at edu- 
cation. For some uses, as in water management, a flow record gives signifi- 
cant results immediately. For the development of a ground-water reservoir, 
the interpretations of well logs, pumping tests, water levels, and well spacings 
require highly involved and elaborate computations. Interpretations of a 
different kind serve to explain the nature of the occurrence and behavior of 
water to the general public, and to students of all ages who wish to enjoy a 
cultural appreciation of nature. The more detail that is available, and the 
ligher the grade of the data, the wider is the range of interpretation that can 
oe made and the more accurate are the applications. 

There are other ways in which water data must be used. In a democratic’ 
society we all share in making decisions about water developments. Pro- 
j0Ssals for this dam or that irrigation scheme engender wide public dis- 
cussion, much of which is not technically well founded. For example, the big 
lam-little dam controversy generates much heat but little light. Wise and 
Iroper decisions can be assured only if the public has direct access to facts 
in a form that they can understand. Raw tabulations of statistics of rainfall, 
streamflow, and ground water, necessary as they are, will not meed the need. 
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They cannot be used by the public until put into a form that bears on problems” 
of public concern. These are proper goals of the basic-data program. The ~ 
demand for specific data by governmental “action” agencies and by private 
groups causes a lag in interpretation of data. This demand has been so great 
that the data-gathering agencies have had little time or opportunity for a 
thorough examination. 

The job is made more difficult because water problems change in empha~ — 
sis. For example, the development of the deep-well pump and increases in 
agricultural yields and industrial demands have combined to make increasing 
demands upon ground water. Yet without reexamination and interpretation, 
the collection of basic data may become so inflexible and sterile that the 
programs cannot readily be adapted to meet changing needs. ; 

Are we sure that we are not building water projects on inadequate data, or 
building homes on flood plains in disregard of known facts? Water needs and 
problems are as broad as the whole modern economy because water is an es- 
sential ingredient of our industrial might and of our hopes and plans for the 
future. Water facts, scientifically interpreted, can tie our plans to the reali-- 
ties of Nature. All of us can recognize a divining rod in the form of a forked 
stick, but it takes a patient skeptic to detect its more subtle counterparts. ¥ 

We have seen that early appeals for collection of basic data by massive at- 
tack as a sole means for solving problems were a matter of concern as early 
as 1867. It must be clear that the usefulness of extensive programs for the 
collection of basic data depends on correct and timely interpretation of the 
data. The matter appears to be one of emphasis. An equally eloquent plea 
can be made both for more analyses of principles and for more quantitative 
data. A basic-data program can go wrong when there is an imbalance on 
either side. 
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1. Outline and Summary 


The purpose of this paper is to describe recent new methods which have 
been devised for dealing numerically with flow problems in large rivers and 
reservoirs for lengthy periods of time of the order of weeks. The basic idea 
of the method is to make use of the differential equations which formulate the 
well-known physical laws governing such flows as a basis for numerical calcu- 
lation of future river depths and velocities. Such a method has become possi- 
ble in a practical way only because of the relatively recent development of 
modern high speed digital computers, and of adequate methods for numerical 
analysis of the differential equations, since the amount of numerical compu- 
tation necessary to solve the equations in question is very great indeed. 

The idea of using the differential equations expressing the laws of conser- 
vation of mass and momentum as a means for treating flows in open channels 
is not at all new. In fact, it goes back to Massau(3)1 as long ago as 1889. 
Since then the idea has been taken up by many others (for the most part in 
ignorance of the work of Massau) for example, by Preiswerk,(4) Thomas, (8) 
Craya, (1) Stoker. (6) Thomas, in particular, attacked the flood-routing 
problem in his noteworthy and pioneering paper 8) and outlined some numeri- 
cal procedures for its solution based on the idea of using the method of finite 
differences. However, some of his methods would not necessarily furnish a 
good approximation to the desired solutions even if a large number of di- 
visions of the river into sections were to be taken. In general, the amount of 
numerical work to be done looked far too formidable to be practical until 
rather recently. 

During the past three years the feasibility of solving numerically such 
problems concerning flows in rivers has been investigated at the Institute of 
Mathematical Sciences of New York University under a contract with the 
Corps of Engineers of the U. S. Army,2 and in cooperation with its Ohio River 
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Division located at Cincinnati, Ohio. In fact, the entire project was greatly 
facilitated by B. R. Gilcrest of the Ohio River Division. Without his constant 
interest, advice, and willingness to help in providing and analyzing the data— 
a crucial task— we would not have been able to carry out our work. 

Three test cases were selected, with the object of comparing the observed 
flood stages with those computed numerically: 


1. The flood of 1945 (one of the largest on record) in the 375-mile stretch 
of the Ohio River between Wheeling, W. Va., and Cincinnati, Ohio, 

2. The flood of 1947 through the junction of the Ohio and Mississippi 
Rivers at Cairo, Ill., and 

3. The flood of 1950 through Kentucky Reservoir, at the end of the TVA 
system of power and flood control dams, which extends from Kentucky 
Dam near the mouth of the Tennessee River to Pickwick Dam 184 miles 
upstream. 


In each case the state of the river or river system was taken at the begin- — 
ning of the flood period, i.e., the river stages and velocities were assumed 
known everywhere at some arbitrarily selected initial time. For subsequent ~ 
times the inflow from tributaries and the local run-off in the main valley were 
taken from the actual records, and the progress of the flood in the main 
stream was then computed for a future period of the order of weeks. In our 
cases the Univac digital computer, built by Sperry Rand Corporation, was 
used. The flood stages—that is, the heights of the river surface above sea 
level—as determined numerically, were then compared with those actually ob- 
served. This paper summarizes the methods used and the results of such 
calculations; for a detailed description the reader is referred to the reports 
of(2) by the authors. 

The success of such an enterprise from a practical point of view depends 
on the answers to two questions: 


1. Can the flood predictions be made accurate and quickly? 


2. Is the cost of preparing and carrying out the numerical calculations 
reasonable? 


The conclusions to be reported here are the resultof a first effort which ought 
to be supplemented by still further investigations aimed at improving accuracy 
and on simplifying the methods, but there is no doubt that both questions can — 
already be answered in the affirmative. We show, for example, in Fig. 1 the 
results of the calculations of river stages at the ends and in the middle of 
Kentucky Reservoir; as one sees the differences in calculated and observed 
stages over a period of three weeks are mostly to be measured in inches. .The 
amount of computing time on the Univac required for calculating the flow over 
the three-week period was a little less than 4 hours. That such calculations — 
can be made is, of course, a matter of considerable practical importance from 
various points of view. Once the basic data for a river, river system, or } 
reservoir have been prepared and coded for a calculating machine it becomes — 
possible to solve all sorts of problems quickly and inexpensively. For ex- 
ample, the effect on a flood wave of closing off a tributary by operating a da 
the relative merits of various schemes for serial operation of a system of 
dams, the preparation of tables showing the influence of changing the operat- 
ing schedules of power dams at the ends of large reservoirs, or the solution — 
of design problems in large reservoirs where, for example, navigation by 
barges is complicated by the occurrence of transient waves due to peaking 
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operations of power plants, are all problems which can be attacked numerical- 
ly. 

These observations refer to the experience gained from the use of the 
Univac computer, which, fast though it is, nevertheless requires roughly a 
minute to compute the flow for an hour in the 375-mile stretch of the upper 
Ohio River, and about one-half minute for a flow of one hour in each of the 
other two cases under discussion here; thus for a flow calculation of 13-1/2 
days in the Ohio River about 6-3/4 hours of Univac computing time were need- 
ed. However, if the much faster operating IBM-704 machine were used, these 
times (which, it should be repeated, resulted on the very first attempts to use 
numerical integration methods) would be cut down in the ratio of something 
like 1 to 15—thus reducing the calculating time needed for a flow lasting for 
weeks from hours down to minutes. Besides, calculating equipment is 
constantly being improved and made more available. It is thus perhaps not 
too visionary to predict that one day—and not too far distant a day either—the 
‘data for the whole of the Mississippi Valley river system will be coded for 
calculating machines, though not necessarily in just the form to be described 
here, ready to deal quickly and efficiently with problems of the most varied 
kinds relating to the system as a whole or to any part of it. 

In discussing the usefulness and practicality of the numerical methods to be 
described here for dealing with large scale flow problems in rivers, it is of 
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course necessary to make comparisons with other methods. Hydraulic engi- 
neers have methods, termed flood-routing, which are relatively simple and 
reasonably accurate for dealing with some flow problems. Flood-routing is 
not always accurate, it is very tedious at a major junction like that of the Ohio — 
and Mississippi Rivers, and it fails entirely in the case of a large reservoir 
such as Kentucky Reservoir. In addition, the flood-routing method requires a : 
good deal of man power to carry out, which is at a premium in the Corps of i 
Engineers, and probably also in the Weather Bureau, which is charged with 

the task of issuing flood predictions. Another method used to study flood 
problems is the method of hydraulic models, which have been built by the 
Corps of Engineers for some portions of the Mississippi valley. (They are 
located at the Waterways Experiment Station at Jackson, Miss.) These models 
are constructed on a large scale (covering acres in some cases). They are — 
not true scale models since one of the essential factors in the problem, the 
roughness of the river bed, cannot be scaled, but must be dealt with empiri- 
cally; such models are therefore really so-called analogue computers. More 
will be said later on this point. 


ee 


2. Description of the Numerical Methods 


The basic laws for flow in open channels are formulated in the following i 
differential equations, in which subscripts denote differentiation with respect 
to the distance x along the channel and the time t:. 


BBs. © CAV) Gece 
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The first equation expresses the law of conservation of mass, the second the ~ 
law of conservation of momentum. In these equations H, the stage, denotes : 
the elevation of the water surface above sea level; V denotes the average ve- — 
locity over a cross-section of the river; A denotes the cross-section area; 
B denotes the breadth of the river at the water surface; q denotes the volume 
of the inflow into the main channel per unit length of channel; and G repre- 
sents the resistance coefficient.3 The symbol | V| means that the numerical _ 
value of V is to be taken: this insures that the term GV lv] represents a re- 
sistance which is in the direction opposed to the velocity. The basic quanti-. 
ties to be determined, the stage H and velocity V, are both functions of the 
distance x along the channel and the time t. The coefficients A, B, and G are 
all functions of the distance x and stage H, while the quantity q characterizing © 
inflows is a function of x and t. (It might be thought strange that the slope of 
the river does not occur in (2.1). However, it is really there implicitly, since 
Hx = yx - § in terms of the depth y of the river and the slope S of its bed.) 4 
All of these coefficients, i.e. the quantities A(x, H), B(x, H),q (x,t), and 
G(x, H), must be known functions for the river or reservoir to be studied. One 
of the major tasks—and no small one, either—involved in solving the problems ~ 
under discussion here is the determination of these coefficients from the basic 
data. The quantities A and B are purely geometrical and thus can in princi- _ 
ple be fixed from topographic maps. The resistance coefficient G(x,H) was 
determined in our cases by using flow data from past floods in conjunction 
with the second of the equations (2.1); the terms V; and VVx, which are small 
compared with the other two in that equation, are neglected for this purpose ~ 
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and thus G is determined from the relation G = -gH,/V|V|. The quantity 
q(x,t) characterizing inflows was known in our cases from the records, since 
we were testing our methods using observed values for past floods; in flood 
prediction problems it would be necessary to obtain this quantity from 
measurements at the mouths of the main tributaries together with estimates 
for the local run-off in the main valley. The details of the methods used by us 
for fixing coefficients are described and analyzed in the reports, (2) especially 
in Report II. 

In addition to the differential equations, it is necessary to prescribe initial 
data and boundary data to have a completely formulated problem with a defi- 
nitely determined solution. As initial data we must have, as is to be expected, 
the values of the stage H and velocity V at some instant of time chosen as the 
initial time at which calculations are to begin. Boundary data are required at 


the ends of the river stretch under consideration and can take different forms. | 


The most reasonable way to prescribe boundary conditions in such a case as 
the upper Ohio River is to give the inflow in the main stream at the upper end 
(in our case, Wheeling) and to use a so-called rating curve, i.e. a relation be- 
tween stage and velocity, at the lower end (Cincinnati) which is the result of 
observations of previous floods.4 In Kentucky Reservoir, on the other hand, 
the appropriate boundary conditions to be imposed are the result of the special 
physical situation there, i.e. that the flow is controlled through releases of 
water at the dams closing the reservoir at both ends; in this case it is there- 
fore the discharge rate of the water at the ends of the reservoir (or, equiva- 
lently, the velocity) which is prescribed for all times under consideration in 
the problem. 

Once the initial data, inflows, and boundary conditions are given, the differ- 
ential equations have a uniquely determined solution for the unknown quantities 
H and V for all future times. As a purely mathematical problem then, every- 
thing is in order. However, the differential equations are nonlinear and can 
not be solved by analytical means in general terms, quite apart from the fact 
that the coefficients are given empirically, and we are forced to turn to 
numerical methods. Even numerical methods were pretty much out of the 
question until recent years because of the intricacy and bulk of the numerical 

_ work necessary to solve such equations over hundreds of miles of a river and 
for periods of several weeks. Fortunately, the equations are exactly analogous 
to the differential equations of gas dynamics, and the methods and calculating 
equipment which were devised for solving such problems numerically during 
and subsequent to World War II are available for solving the flow problems 
under discussion here. 

Thus the equations (2.1) for the determination of the unknown quantities 
H(x,t) and V(x, t) are sufficient to determine these unknowns uniquely once 
appropriate initial and boundary data have been prescribed. However, instead 
of trying to solve the two differential equations exactly for the unknown quanti- 
ties for all values of x and t under consideration, one seeks instead good ap- 
proximate values for H and V at a discrete net of points in an x, t-plane, as 

_ indicated in Fig. 2. For example, in the case of the upper Ohio River, the 
stretch from Wheeling to Cincinnati was divided up into 10-mile intervals Ax 
between net points, and approximate values of H and V were calculated at 
these points. (To have stages and velocities known at points ten miles apart 
is quite sufficient for practical purposes.) At the initial instant t = 0, i.e. on 
the x-axis itself, (which corresponded to Feb. 27, 1945 in this case), stage and 
velocity must be known at all net-points on the x-axis at this time. Once these 
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_ values are prescribed, and the inflows q are given (these are distributed to é 

the nearest net points), the values of H and V are advanced by a time incre- 
ment At through use of the differential equations. Technically, this is done by ~ 
approximating derivatives in the differential equations by difference quotients 
which make use of values of H and V at net points only; the result isa system — 
of algebraic equations which are solved by the calculating machine to yield 
approximate values of H and V at all net points. The values of H and V thus ‘ 
found will be good approximations to the correct values if Ax is small enough 
and if At is such that the ratio of At to Ax is under a certain bound (why this 
important condition is necessary will be explained later). 

In all three cases presented here the interval between net-points was 10 { 
miles, and the time interval At was taken to be 9 minutes. Thus we require | 
data only at discrete points 10 miles apart, but such data must naturally 
represent a reasonable average over a ten-mile interval. Without going into 
details we indicate briefly, with reference to Fig. 2, how the differential 
equations (2.1) were replaced by finite difference equations and how the calcu- 
lation proceeds, starting at the time t = 0, i.e. on the x-axis, when initial 
values for H and V are supposed known at all net points along the river. The 
method used to advance the solution to the next row of points distant At from 
the x-axis is then as follows. A staggered net is used, as indicated on Fig. 2. 
Consider the points L and R, where H and V are known; we give the procedure _ 
for determining their values at P, atime At later. First of all, values forH 
and V at M, midway between L.and R, are computed as the averages of the ‘ 
values at L and R: 


1 
(2.2) Hy = 3p + Hy) , Vy =3lVp+V,) . 


The time derivatives at M are then given approximately by 
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and derivatives with respect to x at M are calculated similarly: 


H, - H Vio - V 
eee ee V0 a) 


(2.4) : Ax ; 4 Ax 


These approximations to the derivatives at M are inserted for the correspond- 
ing derivatives in the differential equations (2.1), the coefficients being evalu- 
ated at point M. The result is a pair of algebraic equations which can be 
solved to yield values for Hp and Vp, as follows: 


2 
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Here qyj,R is the volume of inflow per unit time between the points L and R. 
At the boundaries, a somewhat different procedure is used (see (2)). In 
analogous fashion the solution is advanced from one line of net-points parallel 
to the x-axis to the next line At minutes later. In principle, the idea is simple, 
but in practice a great deal of calculation is needed, since, for example, 1100 
constants must be available to the calculating machine in the case of the Upper 
Ohio River to define the coefficients A, B, G, etc. In addition, since the so- 
lution is advanced only 9 minutes at a time the number of net points, and hence 
the number of values of H and V to be determined is very large for a stretch 
of river 375 miles long and for times of the order of weeks. 

It is conceivable a priori that the amount of computation necessary might 
be so large as to make this scheme impractical. It is even conceivable that 
it might take longer to predict a flood numerically than the time of duration of 
the flood itself.5 Fortunately, as has already been noted, reasonable calculat- 
ing times resulted in the cases reported here. 

The question of the accuracy of the approximation has been studied for 
differential equations of the type in question here. It is known that approxi- 

_ mations to the exact solutions can be made with as much accuracy as desired 
by taking the mesh-widths in the finite difference scheme referred to in con- 
nection with Fig. 2 sufficiently small. However as was already mentioned, it 
is a striking and interesting fact that it is not possible to dispose of the size 
of the mesh width interval At in the time variable in an arbitrary way; instead, 

’ once an interval Ax has been chosen—presumably small enough for the accu- 
racy desired in the approximate solution—the time interval At must not be 
taken larger than a certain multiple of Ax which depends upon the stage H and 
velocity V. If this rule is violated, the consequence is not simply an inaccu- 
racy in the results; rather, totally unusable values for the unknowns H and V 
are obtained. This fact also has decisive practical importance, since the 
amount of calculating machine time needed to solve a given problem is direct- 

_ly proportional to the number of net points and hence inversely proportional 
to the permissible size At of the time increment. For example, in our three 
problems in which an interval Ax = 10 miles was chosen for the space variable 
(on the basis of preliminary calculations in simple cases) the maximum time 
interval permissible was about 10 minutes; we therefore chose 9 minutes as 
time increment to be on the safe side. Since some of the calculated flows had 
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a total time duration of three weeks, it is obvious that the limitation of 9 
minutes as the time step by which the solution could be advanced at any one 
step means that a large number of net points was needed, with correspondingly 
large amounts of calculating machine time. 

The limitation that must be imposed on the ratio of the time increment to 
the space increment is not simply the consequence of a purely mathematical 
fact having to do with the convergence of a numerical scheme to the correct 
solutions; rather, it stems from a fundamental physical property of the flow. 
We are concerned in these problems, as also in the analogous problems in gas 
dynamics, with a wave motion. A characteristic property of wave motions in 
a river (and in the analogous situations in gases and elastic solids also, for 
that matter) is that a small disturbance in the flow created over a limited 
portion of it—an influx of water from a tributary, say—will create a dis- 
turbance (called a wavelet) the front of which propagates both upstream and 
downstream at a rate, relative to the stream flow, which depends on the local 
value of the stage H; in fact the propagation speed of wavelets relative to the 


stream is given by \eq = yea, with d the average depth of the water at a given 


cross-section. In other words, if we interpret this remark in the x, t-plane of 
Fig. 2, it means that a disturbance created in a certain region at a given time 
spreads out in a definite fashion and affects the flow only in definite regions 
of the x, t-plane. Or, upon reversing this consideration to put it in the way 
which is most important for present purposes, the magnitudes of the stage H 
and velocity V at any particular time and place, that is, at any given point of 
the x, t-plane, are determined from the values of these quantities at an earlier 
time only on a definite and uniquely determined segment along the river, and 
are totally uninfluenced by what happens at any other points along the river. 
This has a most important bearing on the problem of calculating a wave 
motion, as we explain with reference to Fig. 3. We suppose that calculations 


Fig. 3 Characteristics 


Pe 
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for a flood wave had already been advanced to a certain time t, i.e. to a 
certain line in the x, t-plane parallel to the x-axis, and that it is desired to 
advance the solution to the next line of net points corresponding to a time At 
later. The new values at the point P, for example, are to be determined from 
the known values at the two nearest points A, B (at the time At earlier) 
through use of the differential equations expressed approximately in terms of 
difference equations. From our discussion of wavelets it seems clear that 
such a method would not work if point P were to fall outside of the region of 
the x, t-plane which is affected by conditions imposed over the segment AB, 
Since if it did, one would not be using all of the relevant data that is required 
physically to fix the motion completely at point P. The curves AC and BC 
which define the region within which P must be chosen are called character- 
istics; they are curves fixed by the solution, which, in physical terms, are 
nothing but the loci along which the wave fronts of small disturbances propa- 
gate—in fact, the curve AC, for example, is a portion of that wave front caused 
by a small disturbance at A which moves to the right relative to the stream. 
Thus, the more rapid the propagation speed of wavelets the smaller is the per- 
missible time increment At, since high propagation speeds mean that the wave 
fronts move along curves inclined at small angles to the x-axis. Thus numeri- 
cal calculation of wave motions is made easier and more rapid if wave propa- 
gation speeds are small; this is the reason it was stated earlier that gas dy- 
namics problems are in some respects more difficult to deal with than flood 
wave problems. 

The flood wave motions in rivers are the result of the nonlinear super- 
position of elementary wave disturbances which move both upstream and down- 
stream from the sources of disturbance—which are, it should be said, not 
caused solely by flows into the stream from tributaries, but also by variations 
along the river in cross section, in roughness, in operating conditions at 
dams, etc. Most engineering methods of flood routing are based on the fact 
that a flood wave in a river is commonly observed to progress downstream 
only, and with a velocity which can be estimated or fixed from data obtained 
from past floods; this observation, together with the use of the law of conser- 
vation of mass, leads to a relatively simple method for calculating the 
progress of the wave. At first sight, this would seem to be in crass contra- 
diction with the analysis made above, which involves waves propagating up- 
stream as well as downstream, particularly when, in addition, the propagation 
speed of a flood wave (defined, say, as the speed of travel of its crest) is actu- 
ally slower in ordinary rivers by a factor of 1:4 or 1:5 than the propagation 
speed of wavelets. There is no real contradiction, however. What happens is 
that the resistance to the flow caused by the roughness of the river bed has 
the effect of slowing down the propagation speed of the main part of the wave, 
which nevertheless is the result of superposition of wavelets traveling in both 
directions. The engineering method thus could be expected to work best when 
the effect of propagation upstream is relatively small, and this is often the 
case for a flood in a long river. However, if there are many obstructions in 
the river, or if tributary inflows, so large as to cause important backwater 
effects occur (such as, for example, the flows at the junction of the Ohio and 
Mississippi Rivers) a method which ignores the finer details involved in the 
method employing wavelets may not be sufficiently accurate. The reason why 
the problem presented by flows in a large reservoir, such as Kentucky Reser- 
voir, have never yielded to the conventional methods of solution is doubtlessly 
the fact that propagation of waves upstream due to reflection from the face of 


the dam at the downstream end, and similarly at the upstream end, is a factor 
of major importance in such cases. In Kentucky Reservoir, for example, a 
wavelet propagates with such a speed as to traverse the length of the reser- 
voir (a distance of 184 miles) in a time of the order of 1/4 day; for a flood 
calculation covering a period of three weeks there would thus be time for 
something like 80 such traverses of waves to be made back and forth between 
the dams at the two ends; and since the resultant effect at any one point comes 
about through a superposition of these waves it is not surprising that it is es- 
sential to consider wave propagation in both directions in this case. 
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3. Description of Results. Comparison with Observations 5 


We proceed to give an account of the actual results obtained in the three 
cases we have treated using numerical methods. 
We begin with the case of the upper Ohio River, shown diagrammatically 
in Fig. 4. Figs. 5 and 6 show calculated versus observed river stages (given H 
as heights above sea level) at Maysville and Pomeroy for somewhat more than — 
two weeks of the 1945 flood; this period includes the crest of the flood. As : 
one sees, the agreement of the calculated with the observed stages is quite 
good, the differences between the two being mostly less than one foot. This : 
case, which was the first one dealt with by us, is, however, the one for which 
7 
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our results are the least accurate, for reasons which will be discussed later 
on. 

Fig. 7 shows a schematic diagram of the junction of the Ohio and 
Mississippi River. The flow velocities into the river at Thebes and Metropolis 
were given, and the stage-discharge relation was given at Hickman on the 
lower Mississippi. Figs. 8 and 9 show calculated and observed stages at 
Hickman and Cairo, and Fig. 10 the same things at Metropolis for a period of 
three weeks during the 1947 flood in the Ohio River. Again the agreement of 
the two is very good, and even minor variations in the curves for the observed 
stages are mirrored in the curves obtained by calculations. The differences 
between observed and calculated stages in this case are of the order of inches. 

Fig. 11 is a diagrammatic sketch of Kentucky Reservoir at the mouth of the 
Tennessee River. In this case the motion of the water in the reservoir was 
almost entirely controlled by the releases of water at Kentucky Dam and 
Pickwick Dam at its ends, since inflows from the local drainage area were 
quite small. Fig. 12 is a graph showing the releases at these two dams during 
three weeks of the 1950 flood. As one sees, the releases varied quite rapidly 
at the beginning—to empty the dams in order to provide storage for the flood 
water coming down. The results of calculating stages at the two ends of the 
reservoir and at Perryville near its center have already been shown in Fig. 1. 
Again the agreement between observed and calculated stages is very good for 
the entire period of three weeks. It is worth repeating that the successful out- 
come of the calculations for Kentucky Reservoir is a fact of particular im- 
portance, since the conventional flood routing methods fail in such cases. 
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It should be emphasized that the amount of Univac computing time required 
in these cases was one minute for determining the flow for an hour in the case 
of the Ohio, and half that in the other two cases. If the IBM-704 were to be 
used these times would be cut in the ratio of approximately 1:15. 


4. Comparison with Model Experiments 


It has already been mentioned that hydraulic models of rivers have been 
built at the Waterways Experiment Station in Jackson, Miss. for the purpose 
of studying flood waves, and that they are large—covering acres of ground. It — 
is of interest to compare this method with the numerical method, for one thing - 
because the two methods have much in common, and besides, a discussion of 
these facts makes it possible to touch upon a number of things that are of 
general interest. 

The models are built with a much exaggerated vertical scale compared 
with the horizontal scale. The contours of the river valley are reproduced 
accurately in concrete. Specially designed pumps have been devised to repro- 
duce given inflows. The river stages at the main gaging stations are continu- — 
ously measured with delicate gauges and recorded remotely on graph paper 
in a special building which houses the electronic equipment. It has already 
been mentioned that one of the main physical features governing the flow in a. . 
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iver, namely the roughness of the river bed, which causes a resistance to the 
low, does not scale simply for a model. Instead it is necessary to introduce 
lbstacles in the bed of the model in such a way that actually observed past 
loods are accurately reproduced; only after this has been done can the models 
used to study floods. At Jackson, these obstacles take the form of brass 
nobs screwed into the bed of the model, and of wire screen. 

The two methods being compared are closely analogous even in details, 
‘he purely geometrical quantities, such as cross-section areas and breadths 
nust be available for both methods; in the numerical method, however, aver- 
ge values of these quantities over ten mile intervals have been proved to 
uffice, and only the area and breadth, not the exact contours, are needed. 
‘he resistance is finely adjusted in a model simply by experimentation: 
rass knobs, or wire screens, are added or removed until the flood stages 
heck with those from actually observed past floods. In the numerical method 
first estimate for the resistance is obtained from records of past floods. 
levertheless, we had expected that modifications of the resistance coefficient, 
ist as in model studies, and perhaps also of the average cross sections, 
mould be necessary in order to reproduce accurately a given flood over a 
engthy period of time by means of a numerical calculation. However, it 
irned out that in only one of the three problems discussed here was a 
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substantial revision in the first estimates for cross section and resistance 
coefficients needed. For the junction problem and the problem of Kentucky — 
Reservoir practically no changes from the values first obtained from the basic 
data for both cross-sections and resistance coefficients were needed to check 
the observed flows. For the problem of the upper Ohio River, however, con- 
siderable revisions, particularly of the average breadths, were necessary in — 
order to reproduce the 1945 flood with the accuracy shown in Fig. 5. Asa 
consequence, when calculations were made for the flood of 1948 in the Ohio, 
it was found that the results, though not bad at many places nevertheless wers 
considerably less accurate—in the vicinity of Huntington, W. Va., particularly 
—for this flood than they were for the 1945 flood. The natural inference is 
that the basic original data available for the Ohio River were less accurate 
than they were for the junction problem and for Kentucky Reservoir. This wa 
indeed the case, particularly with respect to cross section areas; in the case 
of the Ohio River these were determined from storage volumes obtained by — 
balancing flows observed in past floods over the intervals, called reaches, be: 
tween the gaging stations shown in Fig. 4, which have lengths of from 60 up t 
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1early 100 miles. In the other two cases, accurate average cross sections 
vere obtained from topographical survey data. The fact that fairly reasonable 
-esults for the Ohio River were nevertheless obtained indicates that it is 
ossible to go rather far with the use of average properties over large 
listances; however, our coefficients for the Ohio River should be revised and 
nade more accurate. The accurate results obtained in Kentucky Reservoir 
or the flood of 1950 were checked by making calculations for the flood of 
1948, with good results; this was to be expected since the observed flood of 
(950 was reproduced accurately by the calculations with practically no chang- 
Ss in the coefficients obtained from an analysis of the basic data. It is thus 
ndicated that the method of numerical calculation will in general require no 
‘reat revision of the coefficients if average cross sections are obtained from 
Opographical surveys. 


CONCLUSION 


We have seen that the method of numerical computation is flexible. That 
S, it is simple to vary the physical parameters since all that is necessary is 
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to change some instructions to the calculating machine. For instance, even 
cross section areas might need some revision after comparisons of calculated 
with observed floods since the topographic surveys might be so old that the 
river channel had changed. Some rivers, such as the Missouri River for ex- 
ample, can change their characteristics fairly rapidly—because of scouring © 
of the bed, or conversely by silting up, even to the extent that changes in the 
course of the river occur. Such changes would manifest themselves in dis- 
crepancies between observed and calculated river stages; revision of the basic 
data could then be made. In other words, constant back-checking of calculat- 
ed with observed floods could lead to steady improvement in the accuracy of © 
the basic physical data for a given river. Another feature in which numerical 
methods are flexible concerns the method of dealing with the local runoff — 
the flow from tributaries. In the numerical method there is no difficulty of 
any kind involved in feeding flows in at any net points whatever. Furthermore 
once the coding of the basic data for a river has been completed, the neces- A 
sary instructions to the computing machine are placed on a few small mag- ; 
netic tapes which are available at any moment to solve a given problem, 
These remarks should not be regarded as adverse criticisms of model 
studies in hydraulics. For example, it would be difficult if not impossible to 
deal numerically with the details of the flow close to Kentucky Dam itself an 
through its turbines and over its spillways: for such things model studies 
are needed. Even for such a case as a large reservoir the method employing 
models was the only one available until the recent development of high speed 
digital computing machines, and current methods of numerical analysis, ma‘ 
possible methods of the sort described in this paper. 
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FOOTNOTES 


1. Numbers in parenthesis refer to the bibliography at the end of the paper. 


Although this work was initiated by the office of the Chief of Engineers, 

U. S. Army and was administered by the Ohio River Division of the Corps 
of Engineers, the views expressed here are solely those of the writers and 
do not necessarily represent those of our sponsoring agency. 


For a derivation of these well-known equations, see, for examples, refer- 
ences (5) and (7). Also, detailed descriptions of the methods described 
briefly in this paper can be found in the reports of the writers (2) and in 
(7). 


This is what electrical engineers would call a matching impedance con- 
dition, the purpose of which is to replace appropriately the effect of the 
remainder of the river below Cincinnati. 


. In early attempts to make weather predictions numerically, such situations 


actually occurred. 


BIBLIOGRAPHY 


. A. Craya, Calcul graphique des reg imes variables dans les canaux, La 


Houille Blanche, 1946. 
E. Isaacson, J. J. Stoker, A. Troesch, Numerical Solution of Flood Pre- 


diction and River Regulation Problems. 
Report I: Derivation of Basic Theory and mocmdlation of Numerical 


Methods of Attack, 1953. 


Report II: Numerical Solution of Flood Problems in Simplified Models of 


the Ohio River, 1954. 


Report III: Results of the Numerical Prediction of the 1945 and 1948 Floods 


in the Ohio River, of the 1947 Flood through the Junction of the 
Ohio and Mississippi Rivers, and of the Floods of 1950 and 1948 
through Kentucky Reservoir, 1956. 


1810-18 HY 5 October, 1958 — 


7. 


. H. Rouse, Engineering Hydraulics, New York and London, 1950. (Especial 


RELEASES FOR 1950 FLOOD 


’ 
‘ 
‘ 
! 
' 
1 
i] 
' 
‘ 
' 
i 
‘ 
U 


350 discharge 
in 1000 CFS 


s 


300 


250 § 


‘ 
i] 
‘ 
' 
' 
‘ 
. 
' 
’ 
' 
‘ 


Kentucky Dom 


a ae 


YN 4 § 6 7 8 9 10 1 2 13 14 «18 16 «617 18 19 20 2 22 2 26 


Fig ake 


Published as reports of the Institute of Mathematical Sciences at New York 
University under the terms of a contract with the U. S. Army, Corps of 
Engineers, Ohio River Division. 


J. Massau, Mémoire sur Vintég ration graphique des eq uations aux derivees 
partielles, Annales des ingenieurs sortis de Gand, Tome XII, 1889. , 


. E. Preiswerk, Anwendung gasdynamischer Methoden auf Wasserstroémung- 


en mit freier Oberflache, Mitteilungen aus dem Institut fur Aerodynamik, ; 
Technische Hochschule, Ziirich, 1938. P 
a 


rt 


Applied Mathematics, Vol. I, p. 1, 1948. 


J. J. Stoker, Water Waves, Interscience Publishers, New York, 1957. 


8. H. A. Thomas, The Hydraulics of Flood-Movements in Rivers, ee 


Institute of Technology, Pittsburgh, Pennsylvania, 1937. 


ly the article by B. R. Gilcrest, p. 635.) 
. J.J. Stoker, The Formation of Breakers and Bores, Communications on 


Paper 1811 October, 1958 HY 5 


a SSD 


Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 
—————— a 


QUEUING THEORY AND WATER STORAGE 


W. B. Langbein,! A. M. ASCE 
(Proc. Paper 1811) 


ABSTRACT 


A method for determining the amount of holdover storage for regulating 
streamflow is presented based on analogies with queues. Where discharge 
varies linearly with reservoir contents a unique expression is derived. The 
queuing analogy leads to a general solution by “probability routing”. The ef- 
fect of sequential correlation is also discussed. 


Reservoir storage is usually calculated by the mass-curve analysis intro- 
duced by Rippl (1883). This method, or modifications of it, is based on a de- 
tection of periods of the part record when flows were less than required for a 
given need. The cumulative difference between desired flow and the actual 
flow during the period of flow deficiency equals the needed storage to supple- 
ment the available flow. The mass-curve method appears deficient to the ex- 
tent that it yields unique answers that may be deceiving in accuracy. Reser- 
voir design involves some chance that the final decision may prove too small 
or too large. These alternatives, although not often expressed in engineering 
reports, are implied in any water development that is linked to a variable re- 
source like river discharge. There is need, therefore, for improved methods 
of storage design that give tangible expression to the probabilities of shortages 
and excesses. The theory of queues (waiting lines) offers a new approach to 
the determination of storage on a probability basis. 

All of us have experience in waiting on line, and waiting offers a good time 
for thinking about the logic and operations of a queue as it is sometimes 
called. A line waiting for service at a ticket office, for example, is an obvious 
queue. There are other kinds of queues, more subtle in form, but quite simi- 
lar in operation. Airplanes stacked over an airport, store inventories, 
persons waiting for a clear telephone line, all are queues. A dam placed 
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across the river with gates or openings to control the outflow impounds water. ; 
The impounded water is a queue. The inflow to reservoir represents the ar- 
rivals, the regulated outflows represent the departures, and the impounded ( 
water represents the queue. The analysis of queues is a problem treatedin 
the new science of “operations research.” (Churchman, et al, pp. 389-471, 
1957.) However, few of the developments in queuing theory, as it is called, 
are directly applicable to storage analysis. A paper by Moran (1954) recog- 
nized the analogy between queues and water storage and introduced probability 
theory in storage analysis. This paper is intended to develop the subject 
further. 

One of the most useful applications of queuing theory is that it provides a 
classification of the factors that can be used to describe water storage. 


1. Arrival rate. This is generally expressed by a frequency or time 
distribution of the arrivals, the demands for service, or in the case of 
reservoirs, the inflow. 

2. Queue discipline: This is the rule establishing priority of servicing 
among those on the queue. Ordinarily the rule is first come, first , 
served. But special rules may be set up or observed. Quick jobs on 
my desk get attention before long jobs. One acre-foot of water ina i 
reservoir is usually as useful as another; but selected drafts may be i 
made from high or low levels because of a desire to draw off water that 
is cold or warm, turbid or clear, or water differing in salinity. 

3. Service function. This is the rule defining the rate of service of items 
in the queue. The service might be automatically controlled like a 
traffic light or a fixed reservoir outlet. Or the control might be manual 
like a traffic cop or a gate tender ata reservoir. In either case the 
control action is described by what is:called the “service function”. 
The traffic light controls the flow of traffic in time sequence; the reser- 
voir orifice controls the outflow in proportion to the head and hence to 
the storage; the cop controls the flow of traffic in time sequence, but 
hopefully he gives weight to changes in relative density at the cross- 
roads; the gate operator receives instructions from a dispatcher who 
takes account of demands and the volume of storage on hand and in 
prospect. 

4. Attrition rate. The longer a line awaiting service, the greater isthe 
likelihood that persons will either prematurely depart or refuse to join — 
the queue. The attrition rate depends on the kind of service. People hi 
will wait longer for a doctor’s attention than to buy a pack of cigarettes. 
There is a parallel with reservoirs. Here the loss is that due to evapo- 
ration which is also proportional to the amount of storage that exists at © 
any time. The rate, of course, depends on the regional climate and the 
geometry of the reservoir. Decisions on tolerable queues or optimum 
size of reservoir must consider loss through attrition or evaporation. 


4 


> 


’ 


The parallels appear sufficiently close to warrant consideration of this 
method of classification in the design of those storage reservoirs, built to 
store water for subsequent use. Reservoirs may be classified in different 
ways—capacity is one way, detention period (the ratio of capacity to average 
river flow) is another. In queuing theory, reservoirs are classified in the way 
we have just described queues—principally by the inflow distribution and the 
service function. The inflow distribution is generally expressed by the fa- 
miliar duration curve, and queuing theory provides the long-sought means for 
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relating the duration curve to storage. Service functions are less formally 
defined in customary reservoir design. Some are depicted in Fig. 1. Type IA 
includes those where the outflow discharge is linearly related to storage; the 
reservoir is intended to maintain a minimum low flow b, and to decapitate 
floods. Type IB is the orifice type, and is typical of reservoirs used for flood 
control. Type II is the usual type of regulation assumed in most analyses of 
storage by the mass curve or Rippl diagram method (Rippl, 1883). Discharge 
is maintained at a uniform draft unless it is empty or full. The “uniform 
draft” may vary seasonally. In Type IIA, where the uniform draft is less than 
mean (low water regulation), capacity is established such that frequency of an 
empty reservoir is limited. In Type IIB, where the uniform draft is greater 
than mean (flood control), capacity is established such that frequency of uncon- 
trolled spills are limited. Type III shows a seasonal type of rule curve. The 
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Figure 1.--Three types of service functions. 
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discharge varies seasonally and with contents. All of these general classes 

and combinations of them can be handled by probability analysis. The princi- 
pal technical problems of hydrologic design of a reservoir are usually met by 
determination of the frequency with which a proposed reservoir will contain 
differing amounts of water (called the frequency distribution of storage). 

Other factors of interest are the frequency of spill, the frequency that the 
reservoir may be empty, and the frequency distribution of reservoir discharg- — 
es. Some of these applications are developed in this paper. 

Two kinds of solutions are given. The first example is an algebraic so- 
lution applicable only to a linear service function, such as Type IA. This 
special solution is of interest because of the inferences regarding storage that 
may be drawn from it. However, direct algebraic solutions must necessarily 
be restricted because some service functions or arrival rates become cumber-_ 
some when expressed algebraically, a more flexible method, termed “proba- 
bility routing” is given in the second part of the paper. } 

In all the examples to follow, storage is given in units equal to flow in one 
unit of time. Thus, if flow is in acre-feet per year, storage is acre-feet; if 
flow units are in cfs and monthly units are used, then storage is in cfs- 
months. 


Example 1 
Service Function: D = b + kS in which D is the outflow discharge (1) 
and S the storage. 


A change in storage (length of queue) occurs whenever there is a difference © 
between rates of inflow (arrival) and discharge (departures), hence f 


AS=I-D in which I is arrival rate (inflow) (2) 
from (1) and (2) and Fig. 1 (Type IA function) 
AD/k =1-D 


If time intervals are sufficiently close 


De- D, = kIle— kDe (3) 


where subscripts 1 and 2 refer to successive time intervals. 
From Eq. (3) we obtain 


Peghewur 


k | 


De Tek et TKO (4) 
We may also write : 
a ae . eat (5a) 
Dg let T,Ds (5b) 
DJ= “1, * nD, (5c) 


etc. 
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Inserting the value for D from Eq. (4) into Eq. (5a) and for D from Eq. (5a) 
into (5b), etc. we obtain 


\ 
tess + —— + ——— Se 
D, I+k Ine (1+k)? In-2 (l+k)® In-3 | 


nm i+k [Ia 
Hence, each discharge (departure) is a weighted average of the current and 
preceding inflows (arrivals). Since each outflow discharge is made up of suc- 
cessive proportions of the inflows, the standard deviation of discharges, o d 
can be computed from the sum of the squares of the standard deviations of the 
components, hence: 


Ree er Ne) aR 


where @ is the standard deviation of the inflows. 

The distribution of the discharges will therefore be defined by the mean X 
which would necessarily be identical with the mean inflow and a standard 
deviation og. 


(6) 


Storage Capacity Required 
The average capacity follows immediately from the service function: 
D=b+ks (1) 


Since this is a linear relationship, we can insert means and solve: 


Sz hie (7) 


where X is the mean rate of inflow, and S is mean storage. 

But greater capacity than average is required, otherwise the reservoir 
would frequently spill uncontrolled. We can set capacity equal to any desired 
probability that it would be full and spill. Hence, let 


S= S+to, (8) 


Where t is standard measure of probability of a normal distribution, and og 
is standard deviation of reservoir contents. Since Ge 4 we can also 


write 


(9) 


Z 


But the service function can be set up such that minimum probable dis- 
charge from the reservoir 


Serle (10) 
Substituting tog from Eq. (10) into Eq. (9) yields 


Cue 
S == (X-b) (11) 


1811-6 HY 5 October, 1958 
Substituting S from Eq. (7), 
$=2S (11a) 


Under the conditions given, there are some pertinent restrictions to be 


satisfied. (1) No storage would be needed if b is equal to the minimum flow of ~ 


the stream. In this case, the natural flow of the stream would provide the re- 
quired supply. (2) Since the entire flow is to be regulated, there is inverse 
linkage between b andk. For example, if a high value of b is to be set, the 
value of k must be low. These conditions are met by the following consider- 
ations. 

The minimum probable regulated discharge, b, can be expressed 


b=X-ta-~/ a 


The minimum probable natural inflow, m, is 


m=X-to 


Taking equivalent probabilities for minimum inflow and regulated discharge © 


and solving for k we obtain 

. 2(%=ni 

(X-m)?- (X-b)? (12) 
Note that the greater the value of b, the lower is the value of k. 
We can combine Eqs. (11) and (12) 
ei (Xeml- (X-b)? 

MERE BP a 

We can also write 


2 2 
FSS idk! 
oa tam (14) 
These are general equations for storage for given values of minimum 
natural flow, and minimum regulated discharge. Note that storage is zero 
where b = m, and that storage increases rapidly as b approaches X, and be- 


comes infinite when b = X. There is a close analogy with ordinary queues. If . 


attempt is made to equate service capacity to mean arrival rate, length of 
queue becomes infinitely long. (Churchman et al, p. 398 and p. 401, 1957.) 

Although Eqs. (13) and (14) appear quite general in form, in that they con- 
tain the pertinent factors generally considered in reservoir operation, the 
equations cannot be applied indiscriminately. They apply only under the con- 
ditions for which they were derived, viz. 


1. Random inflow, normally distributed. 

2. Service function D = b + kS. 

3. Value of b, minimum probable discharge (at which reservoir is empty) 
has the same probability of occurrence as m, minimum observed 
unregulated (natural) flow. 


Modification for Non-Random Inflows 


The preceding analysis is based on random occurrence of inflows to the 
reservoir. The particular implication of a random model that concerns us is 
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that there is no relationship between sequential flows. In terms of stream- 
flow, this means that a drought is as likely to follow a high year as another 
drought. However, it has been observed that there may be a sequential re- 
lationship—something in the nature of persistence. Abnormalities tend to per- 
sist. Without attempting to prove or disprove the reality of persistence, it is 
possible to include such effects in the storage analysis. 

The first step is to describe the persistence or sequential relationship. A 
simple expression would be r, = r™ where rp is the correlation coefficient be- 
tween items n apart. A more general expression of the sequential effect would 
be r, = rat-1 where r is the coefficient of correlation between successive 
items (n = 1), ra is the correlation between alternate items (n = 2) and so on. 
The values of a and r are less than 1. 

The standard deviation of the discharges as given by Eq. (6) presumes 
independence between successive items. If, however, successive items are 
correlated by r, alternate items by ar, and that between n items by a"r, it is 
possible to revise the equation to include the effect of the linkage. The re- 


vised equation is 
‘, /_k / 26 
So ae Ikan 


The expression ~/ | + Sos is therefore a correction factor (always 


greater than unity) to allow for the effect of serial correlation. Introducing 
the correction factor, c, into Eqs. (13) and (14) for storage we obtain 


(6a) 


_ & (Xm)? (X-b)? 


S= ih (13a) 
and 
2.2 
(Oe — 19 f5| 
Neat (14a) 


It will be observed that the effect of non-random distribution of inflows is 
to increase storage requirements, 


Probability Routing 


Example 1 illustrates an application of discharge routing to a queuing 
problem. This application is a simple one of normal inflow and a linear 
service function. But wherever the frequency distribution of inflows is non- 
normal and service functions are non-linear, a more flexible method is need- 
ed. These are based on a general solution of the queuing equation that is 
adaptable to all cases, including that in the preceding example, if one wished 
to handle it by this general method. The procedure involves a kind of “proba- 
bility routing” in lieu of the more familiar discharge reservoir routing. The 
answers are obtained from a solution of the single-step queuing equation by 
the method of finite differences. The single-step queuing equation seeks 
answers to the question, what is the probability that there will be n units in 
the queue at the end of a given time interval? This probability is the sum of 
the independent probabilities that (1) there were n units in the queue at the 
end of the previous interval and there were no arrivals or departures during 
the time interval, (2) that there were n units in the queue at the end of the 
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previous interval, and that arrivals and departures were equal in number, (3) 
there were n + 1 items in the queue at the end of the preceding interval but 
that departures exceeded arrivals by one unit, and (4) there were n - 1 items 
in the queue at the end of the preceding interval but that arrivals exceeded de- 
partures by one unit. Formal and general solutions are possible from this 
statement, provided several restrictions can be made—the time interval is 
infinitely small, the probability of more than one unit arriving in a time unit 
is negligible, and that servicing and arrival rates can be expressed formally 
and simply. None of these restrictions are applicable in the case of reservoir 
design. However, none of these restrictions are involved in the basic : 
statement of the queuing equation, which can be generalized so as to apply to a 
wide array of problems. ; 


—) © en he egies 


Example 2 


The procedure is best shown by an example. The service function for the 
selected example is defined by two lines shown on Fig. 2. The distribution of 
inflows is shown on Fig. 3. The distinguishing feature of this example is the 
irregular service function and the non-normal distribution of inflows. The 
problem is the determination of the distribution of storage and discharge. 
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FIGURE 2.--SERVICE FUNCTION, EXAMPLE 2. 
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oe 3.--FREQUENCY DISTRIBUTION OF 


IG 
INFLOWS AND DISCHARGES, EXAMPLE 2 


I. Probability of storage equal to or less than 1.0 at time tp. 


1. 


Storage of 1.0 or less can exist at time t,, if storage at time tp-; was 
zero and inflow during this interval exceeded discharge by 1.0 or less. 
Since discharge for a storage of zero is zero (see Fig. 2), and since 
discharge for a storage of 1.0 is 2.05, average discharge during an 
interval in which storage increased from 0 to 1 would be 1.02. Hence, 
inflow during this interval must be 1 plus 1.02 or 2.02 to produce an 
increment in storage from zero to 1.0. The probability of an inflow of 
2.02 or less according to Fig. 3 is 0.27. The partial probability for this 
combination is 0.27 P, which means that if a storage of zero exists, then 
there is a 27 per cent chance that the storage will be 1.0 or less at the 
end of the next interval of time. 

Storage of 1.0 or less can exist at time t,, if storage at time tp-1 was 
less than 0.5 but more than zero, and if inflow during this interval ex- 
ceeded discharge by a sufficient amount. Since discharge for a storage 
of 0.25 (mid range value) is 0.5, and since discharge for a storage of 
1.0 is 2.05, average discharge during an interval in which storage in- 
creased from 0.25 to 1.0 would be 1.28. Hence inflow during this inter- 
val to produce an increment in storage from 0.25 to 1.0 must be 

0.75 + 1.28 or 2.03. The probability of an inflow of 2.03 or less 


-—“ 
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according to Fig. 3 is 0.27. The partial probability for this combination 
is 0.27 Pg-0,5, which means that if storage in the range 0 to 0,5 exists, 
then there is a 27 per cent chance that the storage will be 1.0 or less at 
the end of the next interval of time. ; 
3. Similarly computations are carried out for various combinations. For — 
example, the computations for initial storage between 2.5 and 3.0 are in 
the following paragraph. 
4. Storage of 1.0 or less can exist at time t, if storage at time tp-j was 
less than 3.0 but more than 2.5 if inflow during this interval was less 
than discharge by a sufficient amount. Since discharge for a storage of 
2.75 (mid-range value) is 2.32, and since discharge for storage of 1.0 
is 2.05, average discharge during an interval in which storage de- 
creased from 2.75 to 1.0 would be 2.18. Hence inflow during this inter- 
val must be -1.75 + 2.18 = 0.43 or less. The probability of an inflow of 
0.43 or less according to Fig. 1 is 0.004. The partial probability for 
this combination is 0.004 Pg 5-3,9, which means that if storage is in the © 
range from 2.5 to 3.0, then there is a probability of 0.004 that storage 
will be 1.0 or less at the end of the next interval of time. ; 
5. The results of this series of computations can be written 


a 


Rig? O27 Pt 0.27R_,.«+'0.27P,5_,5¢tO.20R.,<+0.06P.. 5+ O.0166., - O.O04R. ae 


° 
in which Pg is the probability that the storage at the beginning of the 
period is zero. This equation states that the total probability that ~ 
storage at the end of any period (therefore at any time) will be 1.0 or 
less is the sum of the several partial probabilities. 


Il. The same computations are carried for other total probabilities Pp to P7 9 © 
being carried high enough to assure that all plausible amounts of storage 
are included. The equations are listed in Table 1. Note that vertically each — 
coefficient increases from zero to unity. If storage is between 3.0 and 3.5, 
the probability is zero that storage will be 1.0 or less at the end of the next 
interval. It is virtually certainty (= 1.00)that it will be 4.5 or less at the 
end of the next interval. 

Ill. The problem is to solve this set of equations for Pg9.... P7. Recognizing 
that these are cumulative probabilities, we can compute the expressions for 
storage in the successive ranges by successive subtraction, thus 
P9-2.5 = P2.5 - P2.0. The results are given in Table 2. Note that coef- 
ficients add up to 1.0 vertically. Each horizontal line is an equation: there 
is an implied equal sign between the first and second columns. (For ex- 
ample, the second line reads 
P0-0.5 = -02P0-0.5 + .02P0,5-1.0 + .01P1.0-1.5 + .003P1,5-2.0. Thus we 
have a set of simultaneous equations to be solved. The simplest method ; 
appears to be successive approximation. Any reasonable set of values can 
be used for the first approximation. F 

If no other initial estimate can be made, all probabilities can be as- : 
sumed to be the same. The first approximation values are inserted in each ; 


¢ 


equation, to derive the second approximations. These are reinserted to 
derive a third approximation, and the process is continued until the differ- 
ence between inserted and derived figures becomes negligible. Care must 
be taken in the computations that the sum of the probabilities after each 
trial equals one. The values in the column labeled “Ans.” of Table 2, were 
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obtained after a third approximation using a slide rule for the multipli- 
cations. Values of zero in this column must be recognized as nominal 
values. The last column of Table 2 gives the cumulative probabilities in 
the values of Pg, Pg 5, etc. 

IV.Fig. 4 shows a plot of the calculated frequency distribution of storage as 
defined by the results in the last column of Table 2. This gives the answer 
sought. The frequency distribution of discharge from the reservoir can be 
readily determined from Fig. 4 by use of the service function on Fig. 2. 
For example, storage is equal to or less than 2.7 for 30 per cent of the 
time. The corresponding discharge according to Fig. 2 is 2.3. The 
frequency distribution of discharges on Fig. 3 was built in this way. 


Example 3 


In this example, a reservoir of capacity 2 is to be operated so that the dis- 
charge rate does not exceed 1.15. The service function is therefore of Type 
IB. The problem is to define the frequency distribution of storage and dis- 
charges. 

The procedure is exactly as before, except there are upper and lower 
bounds to the storage. The probability equations are given in Table 3. The 
frequency distribution of storage is defined on Fig. 5. The distribution of dis- 
charges, also shown on Fig. 5, requires some further explanation. The dis- 
charges equal 1.15 while there is storage and the reservoir does not spill—in 
this example the probability of spill is very small. The probability of 


© POINTS DEFINED BY CALCULATIONS 


IN TABLE 2. ; . 
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Figure 4.--Calculated frequency distribution of storage, 
example 2 
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Figure 5.--Frequency distributions of inflow, storage, and 
discharge, example 3. 


discharges less than 1.15 are computed as follows: The probability of a dis- 
charge of 0.5 or less is the sum of the partial independent probabilities that 

there is no storage and inflow is 0.5 or less; that storage is in the range 0 - 

0.25 and inflows are 0.37 or less; that storage is in the range 0.25 - 0.50 and 
inflows are 0.12 or less. The calculations are given below: 


Storage Prob.# Inflow Prob. Partial prob. 
0 0.598 0.5 0.11 0.066 
0 - 0.25 .183 “3b 05 .009 
(0.25 - .50 SLIL 11 005 0 
075 


a. From Table 3. 


The probability of a discharge of 0.5 or less is therefore 0.075 which es- 
tablishes one point on the discharge frequency curve. Other points are com- 
puted to define the graph shown on Fig. 5. As a final check the average dis- 
charge must equal the average inflow. 


Example 4 
(Storage Analysis Considering Monthly Regimen) 


The analysis previously described is based on the simple case of uniform 
flow and draft within each unit of time. However, if the time unit is long—say, 
a year—the condition may no longer be met. The general solution might be to 
develop the queuing equations for each month. In the case of years, the proba- 
bility distribution of storage at the beginning of years is the same as at the 
ends. In the case of months, this is no longer true, as there is a seasonal 
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progression in the probability distribution of storage. The matter can be 
solved by considering the probability distribution for the end of one calendar 
month is the same as the probability distribution of storage for the beginning 
of the next following calendar month, Thus, 12 separate sets of queuing 
equations are set up, to be solved as a single system of equations. Tables 
can be set up for each month similar to Tables 1 and 2, except that the so- 
lution now involves a simultaneous solution of all 12 months. The solution 
tends to be long, hence a simpler method permitting use of the year as the 
unit of time for queuing would be useful. 

The record of Cedar River near Landsburg, Wash., from 1896-1950 is used 
to illustrate a procedure. Fig. 6 shows the duration curve of annual flows. 
Cedar River, like streams generally in the Pacific Northwest, has a highly 
seasonal regimen; flow reaches a maximum in the winter and a low in the late 
summer, Fig. 7 shows the monthly median discharges. To emphasize the 
seasonal contrasts, demands, as shown in Fig. 7, are assumed to be minimum 
in winter and maximum in summer, such as might be expected by combined 
water supply and irrigation requirements. Thus, a reservoir would have a 
twofold job—to store water in winter for summer use, and to store water in 
“wet” years to supplement the flow of “dry” years. A reservoir of 250,000 
acre-feet capacity is assumed. The problem is to determine the frequency 
distribution of storage, probability of spilling, and probability of being empty ~ 
for each month of the year. : 

The principal modification is in the basic equation: inflow (I) equals re- 
quired outflow (D) plus change in storage (AS). This equation cannot be applied 
to long units of time during which there is apt to be considerable variation in 
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Figure 6.--Duration curve of annual flows Cedar River. 
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Figure 7.--Median monthly flow and monthly demand, 
Cedar River example. 


inflow and demand, with the likelihood of spillage or an empty reservoir dur- 
ing the interval. Excess flow during spillage is not available to offset subse- 
quent deficiencies in flow nor can deficiencies while the reservoir is empty 
be compensated for by subsequent excesses in flow. The procedure therefore 
is to prepare graphs such as on Figs. 8 and 9 that for Cedar River show the 
net annual change in storage in relation to the annual inflow and the initial 
storage. These graphs are prepared by first determining the average distri- 
bution of the monthly flows during years with various amounts of total flow. 
These monthly flows are then routed through a reservoir of 250,000 acre-feet 
capacity to determine the net annual change in storage. The routing is carried 
out for twelve-month periods ending with each calendar month. Fig. 8 applies 
to the year ending October 31, in the low-water period, and Fig. 9 to the year 
ending May 31 in the high-water period of the year. These diagrams are used 
in lieu of the equation I = D + AS to determine annual inflows required to pro- 
duce a specified net annual change with a given initial storage. The subse- 
quent procedure is exactly the same as previously described, leading to a 
frequency distribution of storage for each calendar month. However, we have 
not yet considered the effect of sequential correlation, so we therefore con- 
sider these frequency distributions as “trial solutions”. 


Sequential Correlation 


The preceding analysis is based on the condition that inflows are random 
in occurrence—that the probability of an inflow as given by the duration curve 
is independent of preceding flows. Without in any way attempting to prove or 
disprove the existence of non-random tendencies in stream flow (other than 
that due to carry over between days and months), probability routing pro- 
cedure can be applied to the determination of storage even where streamflow 
data are correlated in time sequence. 

When inflows are random the probabilities of inflows are independent of 
anticedent discharges or of the amounts of storage on hand. Where inflows 
are related, then it is necessary to recognize that the probability of inflows 
depend on what went on before. If a reservoir is near empty as a result of a 
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NET ANNUAL CHANGE IN STORAGE, IN THOUSANDS OF ACRE-FEET : 


Figure 8.--Net change in storage in relation to initial contents, 
year ending Oct. 31 


period of low inflow, then the probability of high inflow would be less thanif 
the reservoir is full. As was observed.before, the effect of serial correlation ~ 
is to increase the range in storage. | 
The problem therefore is to obtain a clue as to the proper frequency distri- 
bution in the probability routing. It is necessary to define duration curves 
that are conditional on certain flows having occurred. If the coefficient of 
correlation, r, between successive items is computed, then the mean and 
standard deviation of the duration curves can be determined as follows: 


er ee 


Xm=X(I-r) +rX ne 
o'savi-r® 


When Xm is the mean of the duration curve of the flows that follow in the 
time interval immediately after a flow X, o' is the standard deviation of the . 
prospective flows, X is the general mean, and o is the standard deviation of — 
the whole array of flows. Since r is generally small, say 0.3 or less, o' 
equals o very nearly. f 
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NET ANNUAL CHANGE IN STORAGE, IN THOUSANDS OF ACRE- FEET 


Figure 9.--Net annual change in storage in relation to initial contents, 


year ending May 3l 


Another procedure which was followed in the Cedar River example was to 


equal number of years: viz. an upper third (those above 540,000 acre-feet), 
a middle third (those between 540,000 and 450,000), and a lower third (those 
below 450,000 acre-feet). Separate arrays were then prepared of the annual 
discharges in the years following those in the upper third, middle third, and 


divide the array of annual discharges into three groups, each containing an 
| 


lower third. The resulting frequency distributions are shown on Fig. 11. If 


there were no sequential correlation, the three frequency distributions would 
nearly coincide, but as it is, there is a noticeable difference. The frequency 
distribution of the annual flows following those in the upper third is above the 


others, and vice versa for the frequency distribution of the flows that follow | 
those—the lowest third. The dispersion is equivalent to a sequential 


correlation coefficient of 0.28. 
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Figure 10.--Results of storage computations Cedar River, Wash. 


Having completed trial solutions previously mentioned using the single 
duration curve as given in Fig. 6, the probability routing is now repeated us- 
ing the conditional duration curves shown on Fig. 11, as follows: 

The trial solutions of the frequency distribution of reservoir contents are 
used as guides to the selection of the proper conditional duration curves (Fig. 
11) to be used in the probability routing, based on the fact that low, medium, 
or high reservoir contents are associated with low, medium, or high ante- 
cedent inflows. Because in this example the conditional duration curves on 
Fig. 11 were defined by dividing the range of antecedent discharges into equal 
thirds, the range in storage contents is similarly divided. For example, the 
trial solution for years ending on October 31 gave results as follows: 


a ee 
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Figure 11.--Variation of annual flow of Cedar River depending on 
antecedent flow. 


Contents Probability of contents less 
(1,000 acre-feet) than indicated 

0 0.014 

50 .044 

100 .105 

150 .261 

200 -807 

250 1.000 


A graph of the above data indicates that on October 31 contents are less 
than 160,000 acre-feet for one-third of the time, and are less than 190,000 
acre-feet for two-thirds of the time. Therefore, in the probability routing for 
October 31, for initial contents less than 160,000 acre-feet, the lower duration 
curve is used, for initial contents between 150,000 and 190,000 acre-feet, the 
middle curve is used, and for initial contents above 190,000 acre-feet the 
upper curve is used. 

For May 31 the trial solution indicates that contents are less than 250,000 
acre-feet for one-third of the time, and the reservoir is spilling for two- 
thirds of the time. Thus the lower duration curve is used in the probability 
routing where initial contents are less than 250,000 acre-feet, and a duration 
curve intermediate between the middle and upper curves when initial contents 
are 250,000 acre-feet. 
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These computations lead to a second solution of the frequency distribution 
of reservoir contents, If this differs greatly from the first or trial solution, 
a third solution may be indicated. In this example this was not done. A com- 
parison of the results of the first and second solutions for October 31 is as 
follows: 


‘ 


Contents Probability of contents less 
(1,000 acre-feet) than indicated 
First solution Second solution 
0 0.014 0.045 
50 .044 .085 
100 .105 .158 
150 -261 .309 
200 -807 .795 
250 1.000 1.00 


The results of the second solution have been graphed on Fig. 10. The upper ; 
graph shows the probabilities of the reservoir spilling during the several 
months of the year. The probability of spilling reaches about 70 per cent dur- 
ing March, April, or May and is virtually zero during September or October. 
The central graph shows the probability of reservoir contents of indicated 
amounts during the course of the year. For example, reservoir contents in ‘ 
25 per cent of Novembers will be less than about 135,000 acre-feet. Reser- 
voir contents in about 25 per cent of the Mays will be less than about 230,000 
acre-feet. The lower graph shows the probability of the reservoir being empty 
during the several months of the year. The probability of an empty reservoir 
is virtually zero from April through July, and is about 5 per cent during 5 
October. A decision whether a 250,000 acre-foot capacity reservoir is satis- 
factory or not depends on the economic consequences of spilling or of the 
reservoir being empty. Spillage may be the cause of flood damages. While 
the reservoir is empty, discharges will be less than needed. 

The difference between the first or trial solution and the second solution is 
in the direction of increasing the frequency of spilling during the wet season _ 
and increasing the frequency of an empty reservoir in the dry season, in other 
words adding to storage requirements. But this effect can be offset by using 
the correlation tendency as a forecast tool in adjusting reservoir releases. 

In a sense this is already done in practice whenever reservoir releases are 
curtailed during a period of dwindling reserves and increased as reservoir 
contents approach the capacity of the reservoir. 


* 
Evaporation | 
5 


Evaporation draft is a Type III service function (see Fig. 1). It therefore — 
can be handled by combining the intended operational service function with a 
Type III function. If the operational service function is, say, Type I, then the 
combined service function is a series of graphs showing relation between 
storage and discharge varying seasonally. 

If the operational service function is Type II, the combined service function 
is a family of lines varying seasonally. The frequency of storage can be com 
puted directly. 

The frequency of discharge, computed as previously explained, will BAR 
total discharge plus evaporation. As the losses due to evaporation are 
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included in the service function they must be subtracted from the computed 
discharges. Given the season and probability of various amounts of storage, 
the computed discharges can be corrected for the evaporation loss. 


Discussion 


The method described was prompted by considerations of the operations of 
queues. Queuing theory is used as the model for synthesis of the essential 
characteristics of streamflow (such as the frequency distribution and the se- 
quential correlation) in the solution of reservoir storage problems. Analyses 
and synthesis according to the queuing theory offers some advantages. The 
procedure seems to avoid some of the difficulties of storage analyses made by 
use of the traditional mass curve or by trial routing of historic discharge 
through reservoirs. The results of such computations are biased by the length 
of record available. The longer the record, the greater the storage that will 
be indicated—for example, Hurst (1950) showed that storage needed to average 
out a fluctuating discharge increases as the 0.72 power of the length of record. 
The system of “probability routing” described in this paper also avoids the 
vexing question of what initial storage is to be used in the customary discharge 
routing or mass-curve study. Probability routing, moreover, is an exact pro- 
cedure, sometimes called “nonparametric” by which is meant that the combi- 
nation of probabilities is unaffected by the kind of frequency distribution back 
of the probabilities. 

Sequential correlation is a problem not only in this application of queuing 
theory, but in all methods of storage analyses. The problem is more evident 
in queuing theory because sequential correlation must be clearly expressed. 
Sequential correlation is obscured in storage analyses by the mass curve be- 
cause the available record is taken “as is”. Where storage planned is but a 
small fraction of the annual flow, such tacit expression of sequential corre- 
lation may be satisfactory. However, for large and important storage develop- 
ments, the simplicity of this treatment may be misleading and may represent 
inefficient use of the available record which contains only one of the many 
possible combinations that may exist among the river flows. It appears more 
efficient to analyze the record to discern the underlying factors in its compo- 
sition, of which the sequential relation is one, and to use this information in 
a model upon which to erect all possible combinations, rather than to limit e 
oneself to the single combination offered by the limited record. 

The probability routing automatically assures that all possible combi- 
nations of inflows and discharges are reflected in the results. In this way it 
is tantamount to results that would be obtained from an infinite array, subject 
only to the sampling error of the duration curve and the desired refinement 
in selecting limits for the computations. It surmounts the difficult selection 
of a “critical” period for storage design. Reservoir design on the basis of a 
selected “critical” low water period is likely to be quite variable from one 
stream to another. The probability method permits establishment of design 
criteria on a basis of uniform risks. 

The utility of obtaining results of a storage analysis directly in terms of 
probability becomes most evident when operational decisions must be made. 

It will be remembered that the probability curves on Fig. 10 are the resultants 
of all possible combinations of initial storages and prospective flows. In actu- 
al operations, contents at a given time are known and thus the probabilities of 


future contents can be calculated starting with the known storage as a base, 
If forecasts of prospective flows are available, then these can be incorporated, 
provided they are coupled with a statement of the standard error, to enable 
one to construct the probability distribution of the prospective flows. This s 
forecast distribution is then used in lieu of the distribution defined by the 
whole array of past flows. The results of this analysis leads, for example, to 
more definite estimates of the probabilities of the reservoir spilling or being 
empty. These estimates in turn may dictate operational revisions of the 
service function in order to minimize undesirable risks. , 
Once the procedure is understood, the computations are fairly rapid—they 
take longer to describe than to carry out. It is only essential to remember 
when to multiply and when to add probabilities. The basic rules are that the 
probabilities of contingent events are multiplied, those of alternate events are ~ 
added. 
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MECHANICS OF SEDIMENT-RIPPLE FORMATION 


Discussion by John L. Bogardi 
Closure by Hsin-Kuan Liu 


JOHN L. BOGARDI.!—The paper by Mr. Liu should be considered pioneer- 
ing in all respects. The only general remark that can be made in connection 
with Mr. Liu’s paper is that the reader cannot help but perceive a slight feel- 
ing of incompleteness in noting the mere reference to the study of Kalinske, 
without, however, pointing out directly that parameters U,/w and U,4a/¥ are 
identical to 1/t and P respectively although the latter were used over ten 
years ago by Lane and Kalinske. It should be noted at this juncture that as is 
commonly known a number of attempts have been made during the historical 
development of the theory of sediment transportation to introduce dimension- 


Pac Crd 

p f sD 
introduced by Einstein(1) and referred to as flow intensity is essentially 
identical with the channel stability factor d/DS in use since 1942 in Hungary(2) 
and is an invariant proportionate to the square of parameter t. A correlation 
is indicated—and certainly would not have presented difficulties to Mr. Liu-- 
of the multitude of different parameters. Also, a definition is needed which 
includes the interrelation and the application of the parameters to sediment 
movement as clearly as our present theoretical understanding of the problem 
will permit. 

Inspection of Fig. 10 reveals that along the transposed Shields’ curve with- 
in the validity of Stoke’s law, where w~d2 


less parameters. To quote a single example: the parameter 


os dt AEN (1) 


in which N is the approximate slope of the curve showing the variation of 
U,/w with U,d/y whence 


[-N 2+N (2) 


For N = -2 
Uy = eps Constant (3) 


2. Proc. Paper 1197, April, 1957, by Hsin-Kuan Lit. 
l. Prof. of Hydraulic Eng., Univ. of Budapest, Hungary. 
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which means that the movement of very small particles is not affected by the 

shear velocity —that is, by the magnitude of the shear. The correctness of 

this assumption is obvious if it is remembered that neither shear velocity i 

nor shear have an affect upon the movement of small particles. ; 
The fall velocity within the intermediate range from d = 0.1 to 2 mm ; 

varies according to different powers of particle size. Assuming again 4 

Wd) co 43/4 


s 
Lstallan eee (4) 
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if N is chosen as -1/2 as an average 


Uy 7 Cr d @ (5) 
For large particles 
we qd@ j 
mn (6) 
I-N zt N . 
acre 
A 
Since for this range N = 0 i. 
4 
eh 4 
= (7) i 
Uk =C, 42 ; 


Therefore, for large particles, and even for fractions within the intermediate — 

range, shear velocity is proportionate to the square root of particle size. 
Equation 7 is in full agreement with the frictional drag theory of sediment : 

movement since, as it is known that 


} 
r 
; 


Tex We Cid (8) 


where T, is the critical tractive force. 
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HSIN -KUAN Liu,1 A.M. ASCE.—The writer welcomes all the preceding 
discussions. In view of the fact that various comments made by the discus- 
sers are related to two major subjects: (a) the relation between the ripple- 
formation and the transition of the flow near the boundary from laminar to 
turbulent and (b) the comparison between the Shields’ criterion of the begin- 
ning of motion and the writer’s criterion of the beginning of ripple-formation, 
these two subjects will be discussed first. Other comments will be discussed 
according to the order by which they appeared in the journal. 

As can be seen from Fig. 9, the ripples form approximately in the transi - 
tion stage, that is, the shear-velocity Reynolds number is greater than 3 and 
less than 70. Within this range the transition of the boundary from hydrauli- 
cally smooth to hydraulically rough can be understood clearly, if the shear- 
velocity Reynolds number U,d/y is interpretated as the ratio of the thickness 
of the laminar sublayer to the size of the sediment. If the thickness of the 
laminar sublayer §' is assumed to be equal to 11.6 ¥/U,,,(11) at U,d/y = 3, 
the laminar sublayer is about 4 times the height of the sediment; therefore 
the sediment is well covered by the laminar sublayer and the boundary is 
hydraulically smooth. At U,d/Y = 70, the size of the sediment is about 6 
times the thickness of the laminar sublayer, the boundary is hydraulically 
rough, and the flow near the boundary is turbulent. Although this concept is 
very useful in explaining the hydraulic condition of a plane boundary, the 
shear-velocity Reynolds number loses this physical significance once the 
Sand waves such as ripples and dunes appear on the alluvial boundary. 
Therefore it would be better if the shear-velocity Reynolds number is inter- 
preted as an index of instability from the point of view of sediment-ripple 
formation. According to Dryden,(18) as the boundary changes from smooth to 
rough, the laminar sublayer becomes wavy first. This leads to the writer’s 
hypothesis that as the laminar sublayer becomes wavy and eventually disin- 
tegrates into small eddies, at the same time ripples appear on the alluvial 
bed. It is further hypothesized therefore that both the turbulence originated 
from the boundary and the ripples are all caused by the instability of the 
laminar sublayer. Furthermore, the laminar sublayer is considered as an 
intermediate zone having a high velocity gradient between the main flow and 
the alluvial bed. The instability is directly related to (a) the velocity gradient 
represented by the shear-velocity, (b) the state of the disturbance represented 
by the sediment size d, and (c) the degree of damping due to viscosity repre- 
sented by the kinematic viscosity ?. Therefore U,d/Y is important in classi- 
fying the bed configuration. 

The writer has proposed that ripples can form under three conditions (a) at 
the transition of the boundary layer from laminar to turbulent flow, (b) at the 
transition of the boundary from hydraulically smooth to hydraulically rough, 
and (c) the instability of vortex sheet at the interface in case the flow near 
the boundary is already turbulent. 

Fig. 9 indicates that laboratory data collected by the writer are mostly 
related to condition (b). However, a re-examination of the writer’s experi- 
mental data Runs 1-1 to 1-3 shows that these three runs may be related to 
condition (a). To illustrate this point, data of Run 1-3 and Fig. 8 will be used. 

According to Schlichting’s neutral instability curve,(19) the critical length- 
Reynolds number Uyx/Y is 1.12x 105 at the point of instability for flow along 


1. Assoc. Civ. Eng., Colorado State Univ., Ft. Collins, Colo. 
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a flat plate at zero incidence, where the flat plate has a very smooth surface, — 
and the approaching flow has a low degree of turbulence. The corresponding ~ 
wave-length is A = 22.4 6*, where 6* is the thickness of the laminar bound- — 
ary layer at x, and is equal to 1.73 / xy/Ua, where x is the distance from the 
leading edge, and Uq is the velocity at the edge of the laminar boundary layer. 
The critical length-Reynolds number for flow in an open channel having an 
alluvial boundary may be smaller than 1.12 x 105, The effect of boundary 
roughness on the instability which is not fully known is neglected at present. 
Ripples shown on Fig. 8 were about at a distance of 8 inches from the en- 
trance. If the corresponding velocity Ux is taken as the mean velocity of the 
flow which is estimated to be 0.78 fps, at the place where the first ripple was 
Unt _ 0.78 x%% _ 4300 


found Pe nen eee 


which is in the same order of mag- 


nitude as Schlichting’s value. The thickness of the corresponding laminar 
Fl - 
boundary layer is a. 2472 Ve = 5.54 xo 3 f¢ the corres- 


ponding wave-length is = 22.4 6 gi 0.124 ft. = 1.5 in. As shown in Fig. 8, 
the measured wave-length is about 2 inches. 

If the effect of roughness on the transition from laminar to turbulent flow ~ 
is considered, according to Schlichting (with reference to some Japanese f 
data) the critical shear-velocity Reynolds number at the stage of instability 
is U, k/y = 15 where k is the height of roughness. Data of Fig. 8 listed as 
Run 1-3 indicate that ripples begin at U,d/y = 3 where d is considered to be ' 
equal to k. The calculated shear Reynolds number although smaller than 15, © 
is of the same order of magnitude. On the other hand U,a/Y = 3 also means 
the beginning of the boundary transition from hydraulically smooth to hydrauli- 
cally rough. Hence there is a question as to whether ripples of Fig. 8 are : 
caused under condition (a) or condition (b). The writer thinks that this ques- 
tion cannot be settled till the theory of instability for open channel flow havin 
an alluvial boundary becomes fully known. Nevertheless it is rather clear 
from this discussion and also Fig. 9, that ripples are associated closely with ~ 
the transition from laminar flow to turbulent near the boundary. j 

It is interesting to know that Smoltczyk(37) of Germany has also demon- _ 
strated in 1955 that the Schlichting’s instability concept canbe usedinthe 
prediction of ripple length caused by the currents. He found that the wave 4 
lengths computed as an average of the upper and lower branch of Schlichting’s 
curve was of the same order of magnitude as the ripple length obtained by him 
in an experiment. The approach of using the instability of laminar boundary © 
layer to demonstrate the interfacial wave has gained popularity in recent - 
years, for example Wuest(36) of Germany related Schlichting’s theory of in- _ 
stability to the generation of surface waves caused by wind. | 

The third case that the interfacial wave is caused by the vortex sheet phe- 
nomenon cannot be demonstrated as yet by laboratory data, because ordinari-— 
ly when the Reynolds number is high there are some other factors such as 
surface waves, which can also cause sand waves. 

The comparison of the beginning of bed-load movement, and the beginning 
of ripple formation raises some very interesting questions. For example, 
Drs. Vanoni and Brooks propose that the ripple-curve should be moved down 
ward to coincide with the Shields’ curve; Messrs. Albertson, Simons, and 
Richardson propose that the ripple-curve should be modified so that it has th 
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same trend as that of the Shields’ curve. The writer feels that both the 
Shields’ curve and the writer’s curve are obtained from experimental data. 
Therefore they can be modified if there is new evidence from laboratory 
data. Since the proposed changes are not based upon experimental evidence, 
the writer thinks that such changes are not justified at the present time. 

Most observers have obtained, from a statistical point of view, the dis- 
tinction between the beginning of motion and the beginning of ripples; there- 
fore these two curves should not coincide, when U,¢a/Y is less than about 90. 
Furthermore if it is assumed that the ripples will form as soon as the move- 
ment begins, and since the beginning of motion is very difficult to be defined, 
especially for mixture, the writer would propose for practical purpose that 
the curve of beginning of ripples should be used in the study of sediment 
transport rather than that of the beginning motion. Although there is certain 
relationship between the beginning of motion and the beginning of ripples, it 
is not necessary that the curve for the beginning of ripples should have the 
same trend as the beginning of motion as suggested by Messrs. Albertson, 
Simons, and Richardson. 

Before comparing the parameters of Shields’ and those of the writer’s, 
the methods of dimensional analysis will be used to demonstrate all the im- 
portant variables pertinent to this phenomenon. The beginning of movement 
and the beginning of ripples depend upon the forces of the flow on the sedi- 
ment particle and upon the resistance of the particle to the flow. The vari- 
ables for the flow are: D - depth of flow, S - slope of the flow, Tp - boundary 
shear, uy - local velocity, g - gravitational constant, wt - dynamic viscosity, 
P - density of the flow, and k - roughness height of the boundary. Since the 
slope of the flow is defined if the depth of flow and the boundary shear are 
given, slope can be omitted. Since the velocity distribution uy is a function 
of To, P, HM, k, it can be omitted also. The roughness height in this case is 
represented by the sediment size d. The variables of the sediment are: d, 
sediment size, shape factor, and the Ate. the difference of specific gravity 
between water and the sediment. For given sediment size, given specific 
gravity, the shape factor can be evaluated if the fall-velocity of the sediment 
w is known; therefore, the shape factor can be represented by the fall- 
velocity. In summary we have the following variables pertaining to both the 
flow and the sediment. D, Top, g,P,u,d,w,A¥g. Through the use of 7 
theorem, we have the following dimensionless term if To, d and P are used 
as the repeating variables: 


(To ek ats 
Pte AE) RE AEE, 1) NE 


In which d/D is called the relative roughness; , is called channel stability 


factor by Dr. Bogard;(43) VT/,P d/y, the shear-velocity Reynolds number, 


ay 
Bor w, the movability number; To/AY st, an inverse of the flow intensity. 


Bice that the inverse of the flow intensity for the beginning of motion is 
called the coefficient of critical tractive force by Shields. If T/A 7 sq is used 
for water flowing over sand bed, it is proportional to the inverse of the chan- 
nel stability factor d/ Rps (iE ed/ Uy 2), where Rp is the hydraulic radius of the 
sand bed. 
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The first two terms can be neglected if we re-examine the variables for _ 
the flow. It is generally known that the total depth of the flow does not effect 
the velocity distribution near the boundary according to Karman-Prandtl’s 
theory. If the effect of secondary circulation is neglected, the ratio of the 
flow-depth to the width of the channel can be neglected; therefore the depth D 
can be omitted in the consideration. The factor g can be omitted if the sur- 
face wave of the flow is not important which is generally true for the case of 
beginning motion and beginning of ripples. If we omit the D and g as the ’ 


Ee 


variables, we have the following dimensionless parameters 4 
pd, ik _)=0 2) 
p ( oa , sinha Aa cal / 


where U, = {Ts ° i 
Equation 2 is the general equation for the beginning of motion or the beginning 
of ripples. It should be noticed that both D and g may affect the velocity dis- 
tribution near the bottom, if the surface waves are present in a flow of shal- — 
low depth. Under this condition, Eq. (1) should be used. If the sediment is 
assumed to be spherical, Eq. 2 can be further simplified. For given sediment 
size and specific gravity at a given temperature, the fall-velocity is deter- s 
mined; therefore only one of the two: the fall velocity w or the difference of : 
specific weight 47g is necessary. If w is used, we have the following dimen- — 
sionless parameters for the beginning of ripples or the beginning of motion, — 


p( Ud, Ue )aeo 9 
v wv 

which was used by the writer. If the AY, is used, we have the following 4 
equation : 
a ri : 
ri es ae | 
q 
which was used by Shields. It can be shown that U,/w can be computed if the i 
shear-velocity Reynolds number U,d/Y and T/AY sdare given, or TAY sd i 
can be computed if the shear-velocity Reynolds number and U,/w are given. 2 
The following was suggested to the writer by Dr. Iwagaki + 
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oat t Ao aay 


Ue. Ugd wed (8) 
vw v ye 


The fall velocity can be expressed as shown in Eq. 5; therefore, Eq. 6 can be 
obtained accordingly, where Cp is a function of the fall-velocity Reynolds 
number of the particle, if the shear-velocity Reynolds number is divided by 
/T,/A¥ 4d as shown in Eq. 7, the ratio is a certain function of the fall- 
velocity Reynolds number (see Fig. 12). Therefore, for a given ratio 
U,¢/Y / /T)/AY gd , there is a definite value of the fall-velocity Reynolds 
number. By knowing the fall-velocity Reynolds number and the shear-velocity © 
Reynolds number, U,/w, the movability number can be obtained by Eq. 8. ; 
Although it has been shown that Eq. 3 and Eq. 4 are interchangeable for 
spherical particles, it can be demonstrated, however, that Eq. 3 is preferred. 
If we recall that Shields derived his criterion for the beginning of motion by 
basing upon the concept of laminar sublayer. According to the concept of 
laminar sublayer, when the shear-velocity Reynolds number is approximately — 
greater than 70, the laminar sublayer does not have any more effect on the 
velocity distribution. However referring to Shields’ criterion of the beginning © 
of motion, we can find that when the shear-velocity Reynolds number is i 
greater than 70, it still has effect on the coefficient of critical tractive force 


iy . This fact is a mislead to the concept of boundary layer which was 


AYsd 
used in the derivation of Shields’ criterion. Shields has recognized such a 
descrepancy, and attributed to the non-quadratic resistance of the grain in 
the bed. In other words the drag coefficient pertaining to the grain at this 
stage is still a function of the fall-velocity Reynolds number. } 

The coefficient of the critical tractive force is changed to movability num- ; 
ber U,/w as shown in Fig. 10. When the shear-velocity Reynolds number is 
greater than 70, the transposed Shields’ curve is independent of U,d/ Y, which — 
is consistant with the theory of boundary layer. This means that when the 
boundary is fully turbulent (or completely rough), the critical tractive force 
for the beginning of motion is dependent on the fall velocity of the particle. 
It implies that the viscosity may still have a slight effect on the determina- 
tion of the tractive force. For example when the shear -velocity Reynolds : 
number is 70, the corresponding fall-velocity Reynolds number is around 550 
of which the drag coefficient is still a function. 

To demonstrate further that the coefficient of critical tractive force is not — 
a very suitable parameter for use, we can examine again the Shields’ plot, t 
Fig. 13. When the shear-velocity Reynolds number is less than 3, Shields 
extends his curve to a straight line, where the equation is: 


Te 2 Def (9) 
Atvadse XG 


v 


or 


r¢ y, i 
7c = gl” [o/ v (ar)] i; (10) 
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which means when the sediment particle is very small, the shear is independ- 
ent of the sediment size. However, White(31) and Iwagaki(32) have shown that 
for small shear-velocity Reynolds number, the critical tractive force is pro- 
portional to the sediment size as shown in Eq. 11 


eae. F d tanh (1) 


where C is a coefficient, ¢ is the friction angle pertaining to the grain. The 
difference of results between Shields’ and White’s is very significant, yet can- 
not be explained according to Fig. 13. If the movability number is used as 
shown in Fig. 10, the conclusions of White’s and Shields’ can be obtained 
separately, depending upon the slope of the extended curve at very low 
Reynolds number. 


If 
ee eee eS oe 
Ww va (Beh 2 (12) 


Where U,,¢ is the critical shear-velocity 
Therefore 


3 
ae a, ae 
Within the range of the Stoke’s law 7 


1 ated ® (14) 


idee e C, Ass y (15) | 


Therefore 
mp 
Ge)? Cg ater] (16) 
the coefficient of 0.1 in Eq. 10, corresponds a value of Cz = 1.8. 
On the other hand if 
C2 


hed ch) 


v 


= 


Upe 
Ww 
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ie BS aed (18) 


Substituting Eq. 14 into Eq. 18 


ve Cae A 
‘7 ig eS See 
7 i 


ay 
ne 


16 
l (19) 
ie Salve! 
/8 
according to Fig. 10 Cg = 1.3, therefore 
Be Ge73 (20) 


The movability number U,/w is essentially proportional to the square root 
of the coefficient of tractive force TAY sd, therefore the scatter of data is 
reduced if U,/w is used. However as demonstrated above that U,/w is better 
than T),/AY sd for several reasons. Furthermore, recent evidence by Ippen(30) 
and Durand(38) also demonstrates that the fall-velocity is a very suitable 
parameter in the study of bed-load movement. It is not clear however how 
the fall velocity is a better parameter to be used as far as the mechanics is 
concerned. Therefore further research on this parameter is urgent. 

In order to use Fig. 10 for determining the beginning of motion or the be- 
ginning of ripples, it is necessary to use trial and error methods; because 
the shear velocity is involved in both coordinates. However such a difficulty 
can be overcome if the shear-velocity Reynolds number is changed as in the 
following 


Ud _ 4. wa (21) 
vp » 


Therefore Fig. 10 can be changed as shown in Fig. 14 from which the criterion 
can be obtained very easily if the sediment size, fall-velocity and tempera- 
ture are known. Similarly Shields’ coefficient of critical tractive force is 
plotted in Fig. 15 according to the corresponding fall-velocity Reynolds num - 
T 
c 
ber. Comparison of Fig. 15 with Fig. 13 will illustrate why Aysd is not a con- 


stant when the shear-velocity Reynolds number is greater than 70. At this 
point the shear-velocity Reynolds number is about 550 which is not the 
Reynolds number at which the drag coefficient is constant. 

Since the change of water temperature is not appreciable under ordinary 
conditions, the effect of viscosity on the criterion of the beginning of ripples 
is not very important. For engineering practice the criterion for the begin- 
ning of ripples and that for the beginning of motion have been reduced to be a 
function of the sediment size only as shown in Fig. 16. It is shown that when 
Sediment size is about 2 mm, ripples will not form, although some other type 
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of interface waves may occur when U,d/Y > 10~ 80 or approximately 
d>2mm. 
Plate,(28) based upon the present knowledge of velocity distribution near 
the boundary in a turbulent flow and the mechanics of sediment transport, h 
proposed some emphirical formulas for the beginning of motion and for the 
beginning of ripples as shown in Fig. 18. The basic procedure of his appro 
is as follows: 
since 


- ; 2 (22 
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Transition Boundary 


Fig.| 7— Definition sketch for Velocity 
Distribution of a Turbulent Flow 


Near a Boundary 


where at y = yo, u = 0 (see Fig. 17), and k is the universal constant. 


At ps ae 4 md Uzw (23) 
therefore 

w —_ 

Tee ae a P. (24) 


‘a) For rough boundary, if the local velocity being equal to w is assumed to 
be at the grain level, that is 


Yd = A- Y, (25) 
und 
= 3 of (26) 


Substituting Eqs. 25 and 26 in Eq. 24 yields 


Ww = 3S (27) 
Ue 


b) For transition boundary, m is assumed to be proportional to Uy /” ics; 


f= F7" Us ee 


Vv 


ubstituting Eq. 28 into Eq. 24 yields 
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Hence based upon the assumption that when the motion of sediment begins, 
the local velocity at certain level y, + md is about equal to the fall-velocity, 
and using Karman-Prandtl velocity law, Plate found some empirical equations 
for the beginning of ripples, and the lowest possible case of the beginning of 
motion as shown in Fig. 18. 

The writer agrees with Mr. Maddock that it is important to study the tran- 
sition of the boundary from hydraulically smooth to hydraulically rough, be- 
Cause it is the prerequisite for the understanding of mechanics of sediment 
transport, although most of the sediment transport in the laboratory and part 
of the sediment transport in the natural rivers are connected with this tran- 
sition region, our present knowledge of fluid mechanics is not sufficient to 
deal with this phenomenon. Furthermore the writer agrees with Mr. 
Maddock’s statement, 


“The lengthening of the ripple form as the flow changes from laminar to 
turbulent in terms of the Reynolds number as noted by Shields might be 
explained in an over-simplified fashion as being due to the changes in 
the cycle of formation of vortices in the boundary layer.” 


This simplified fashion well serves the purpose for explaining the mechanics 
of ripple formation for practical engineers. 

Mr. Maddock noticed that in Fig. 10 the ration U,/w is proportional to 
/ Cp; this can be shown as follows 


Ue / 


Ww Ud N (30) 
( = ) 
(4 ) |4N / 
wi pees 
For flow within the Stoke’s law 
| 
é. nice er (31) 
¥v 
U i J4+N 
therefore et N (32) 
() SG, 
aN = 1, 
Ji. 


This is the case which White stated that the shear is proportional to the size 
of sediment for the beginning of motion. 
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On the other hand if N = 2 


P xy, (34) 
ik 


which corresponds to the Shields’ result. 

Mr. Maddock is correct in stating that when there is considerable move- 
ment of sediment, the Karman-Prandtl Equation for smooth boundary may not 
be applicable, because the laminar sublayer might be destroyed due to the 
sediment movement. 

Mr. Tinney questions the meaning when the shear-velocity Reynolds num- 
ber exceeds 120, the ripples will form as soon as the motion of sediment be- 
gins. There seems to be confusion due to the writer’s definition of ripples. 
As shown in Fig. A, dunes can form at 100<U ~a/Y < 400, at the beginning of 
motion. Those which form at U,d/Y > 400 have not been classified clearly. 
The reason that interfacial waves will be formed as soon as the motion begins 
at U,d/y >400 might be due to other factors such as surface waves which are 
associated with high values of U,d/ y. Recently Barton(39) has attempted to — 
classify the bed configuration according to the appearance and geometry of the 
interfacial waves. 

Drs. Vanoni and Brooks have raised some very vital questions pertaining 
to the writer’s paper. These questions are probably most typical in the minds 
of engineers. Excluding those questions discussed previously, the answers to 
the remaining questions are as follows. 

They question the basic assumptions that (a) a sand bed behaves like a 
dense fluid, and (b) sand waves and interfacial waves are basically caused by 
the same fluid forces. The writer recognizes that the reaction of the sand 
bed is probably not the same as a denser liquid. However, the kind of reac- 
tion may be different yet the basic principle, such as Newton’s law of motion 
remains the same. For example, an ideal fluid is non-viscous and a real 
fluid is viscous, one may state that therefore the study of an ideal fluid flow 
might be meaningless and should not be linked with a real fluid flow. How- 
ever, the study of the ideal fluid flow has contributed tremendous progress to 
those of the real fluid flow. The writer’s assumptions in dealing sand waves 
may be considered over-simplified, nevertheless they may be very useful in — 
the study of sediment motion. The writer emphasizes the points of the simi- — 
larity; Drs. Vanoni and Brooks emphasize those of dissimilarity. : 

They question the necessity of using stability theory for deriving the 2 
criterion for ripple-formation and propose to use the dimensional analysis. — 
The method of dimensional analysis has been known to hydraulic engineers 
for some time; however, it is not possible to use the dimensional analysis to 
get all the answers the engineer wants. In order to use the dimensional an- 
alysis, it is necessary to have some knowledge about the physical meaning 
the mechanics of the phenomena. Otherwise we will get many unnecessary 
dimensionless parameters. 

They have questioned whether there is an interval between the beginning 
motion and that of the beginning of ripple formation. Mathematically it may 
difficult to prove that such an interval does exist. It is true that the bed con-. 
figuration will be very stable or is changing once the ripples are formed. 
However, the data used by the writer are at the point of beginning of ripples 
and are not after the ripples have begun. Therefore, the corresponding shea: 


ASCE DISCUSSION 1832-21 


will not be much more than the incipient shear which is necessary for the 
incipient motion. Tison(6) has found that sediment may be transported along 
1 plane bed if the flow is laminar. This is another evidence that ripple- 
formation should not coincide with the beginning of motion. The writer ob- 
jects to moving the curve for the beginning of ripples down to that of the be- 
sinning of movement as suggested by them. 

Drs. Vanoni and Brooks do not concur with the writer that dunes will dis- 
appear when the flow conditions fall below the curve in Fig. 9. They use the 
field observation as an evidence, that is, in a natural river if the flow or 
stage were reduced below the beginning of motion, the dunes would remain 
inchanged. This situation may happen when the flow is suddenly reduced be- 
low the critical stage of sediment movement. However, if the low flow which 
is capable of transporting sediment is maintained for some long period of 
ime and at the same time all dunes are affected by this low flow, the writer 
delieves that eventually all the dunes and ripples will be smoothed out by the 
flow. This can be shown in Fig. 9 by the laboratory data of the U. S. Water- 
way Exoerimental Station and Kramer’s data designated as falling stages. 
These data are included in Table 3. 

They seem to agree that the mechanics of ripple formation can be ex- 
slained by the theory of Sir Inglis, who favors the theory of turbulence effect. 
Messrs. Albertson, Simons, and Richardson correctly stated about the effect 
of turbulence on the formation of ripples, 


“----Such turbulence in turn is directly related to the instability of the 
laminar sublayer. In other words, turbulence and instability of the 
laminar sublayer are intimately related and cannot be separated--” 


The basic concept of instability also implies that when the eddies caused by 
he accumulation of sediment are initiated, and when the flow reaches certain 
mstability, such an accumulation of ridges or sam waves will be amplified. 
Jn the other hand if the flow does not reach the degree of instability, any such 
*xistence of wave-formation will be damped out eventually. Therefore the 
sffect of the eddies behind the sediment grains may act as an instigator, but it 
is up to the flow near the boundary to amplify or to damp its existence. Ac- 
sording to Drs. Vanoni and Brooks, ripples are mainly caused by turbulence; 
1owever, they noticed that in the writer’s experiments, except Run 1-4, the 
Reynolds number was rather low. The lack of fully developed turbulence for 
he ripple formation does not necessarily indicate that the writer’s results 
ure opened to question as stated by them, but rather indicates that the theory 
9y Drs. Vanoni and Brooks is questionable. 

They think the ripples shown in Fig. 8 should be called antidunes; those of 
Fig. 6 should be called oblique waves. According to them the reasons that the 
sand waves of Fig. 8 should be called antidunes are (a) the ripples are round 
n profile and (b) the corresponding Froude number is 1.3. The writer thinks 
hat the conditions for the formation of antidunes are (a) sand waves moving 
oward upstream, (b) flow of shallow depth with surface waves (c) high Froude 
umber, (d) fully developed turbulent flow, and (e) constant supply of sediment 
ransported. For the case shown in Fig. 8, there are three of these five con- 
litions missing: (a) ripples did not move upstream, (b) there was no equilib- 
‘ium supply of sediment and (c) the flow was in the transition region, and was 
ot fully developed. The fact that Fig. 8 shows the round crest of the ripples 
ire not necessarily the conditions of antidunes. These round crested waves 
‘an also be caused by undulation of the flow pattern near the bottom. 
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TABLE IV | Data of Exp 


Investi- Sediment Specific 
gator Size, Gravity 
+ Pe 
8.8 
7.6 
ray 
6.8 
5.0 
4.6 
2.4 
46.0 
14.5 
27 6-41 do 2265) T5.5 14,5 
at 42 15.4 
27 a9 do 2.65 15,9) 14.5 
27 10 15.9 ] 
27 8-45 do 2.65 17.4 14.5 
27 46 17,3 
27 9-9 do 2.65 15.0 14.5 
27 10 14.9 
7X 10-34 do 2.65 14.5 14.5 
27 35 14.6 
27 Ha-20=5 0.908 1.85 25.9 9.6 
27 6 26.0 
27 H3-21-4 29 1,74 14,9 12,4 
27 5 15.7 
rH 22-6 do 1.74 21.0 12.4 
oat | is ZLeU 
27 23-5 do 1.74 22.1 12.4 
at 6 22.2 
7A} C,- 24-9 3.300 1535 25.6 16.7 
27 10 25.8 
27 C,- 26-4 1,097 1335 B35 3 6.1 
ah 5 se Re 5) 
27 Zt =2 do Pgs 10.9 6.8 
27 ‘ cS 11.0 ; 
27 C3-28-9 3.197 Ie 32: e350 15.3 
27 10 23.0 
27 C3-29-6 ie As ir é Wy. 23:0) 15.3 
27 7 . 23.0 


- *10-5-10 indicates Sand No. 10, Table 5, Run No. 10, etc. 
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3 on Ripple-Formation 


Av. Sheer 
Velocity 
Ux 
cm/sec 


ace3 0, 235 15.23 66 4.96 0.0523 
2.24 0.255 14.2 56 5.01 0.0571 
oni 0.281 10.3 38 4,54 0.0535 Average 
2.16 0.300 a2 34 4.66 0.0569 incipient 
2.01 0.296 8.6 30 4,04 0.0517 ripple 
L56% 0.331 5.8 17.23 2.79 0.0498 formation 
1,70 0.370 Onze Lies 2.89 0,0573 
1,50 0.654 2.8 4,5 2a2> 0.0660 
5.99 0.130 250.0 1860 S529 0,0545 
9 2.44 0.168 20.9 122 5.95 0.0376 Smooth 
9 local waves 
1 2.54 0.176 rae) 123 6.45 0.0417 local waves 
7 Smooth 
2 2.61 0.180 oie 127 6.80 0.0430 Smooth 
0 local waves 
1 2.58 0.178 23.4 134 6.66 0.0421 local waves 
6 Smooth 
rA 3001 0.208 257 125 9.05 0.0571 Smooth 
0 local waves 
7 2.65 0.183 22.4 120 7.01 0.0450 local waves 
4 Smooth 
3 Zee 0.235 23.4 100 5.06 0.067 Smooth 
2 waves 
5 1.92 0.155 22.7 136 3.69 0.0384 Smooth 
a . waves 
4 2.25 0.182 30,3 162 5.06 0.0526 Smooth 
2 waves 
" VAR AS, 0.182 3154 166 5.06 0.0526 Smooth 
J waves 
2, PAR B 0.134 80.5 630 4,62 0.0408 Smooth 
3 waves 
j 1.54 0.253 14.4 5365: PAREN 0.0632 Smooth 
L waves 
: 1.57 0.258 L356 48.6 2.46 0.0637 Smooth 
5 waves 
| 2.12 0.136 1240 504 4,50 0.0448 Smooth 
J waves 
} 2.32 0,151 79 503 Cee! 0.0535 Smooth 
i waves 
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Refer - 
ence 


Investi- 
gator 


Run No. 


C,-30-5 
6 


4 
Casey h-1/200-14 2,455 
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TABLEIV Data of Ex 


Sediment 
Size, 


15 
1/400-11 2.455 
12 
1/800-17 2,455 
18 
Jorissen,A.L, I-I-2 0.720 2.67 22,2 lites 
3 yey 
I-I-4 0.720 2.67 22.8 lie 
5 20.5 
Il-ll-14 0.930 2.67 22.8 15. 
15 2303 
W-U-17 0.930 2.67 22.8 15.0 
18 20.0 
l-0-18 0.930 2.67 20.0 15 
19 210 
W-l-23 0.930 2.67 20.5 15 
24 
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nts on Ripple-Formation -- Con't. 


Tractive 


0.0457 Smooth 
252 waves 
54 1,62 0.184 252.6 125 2.62 0.0562 Smooth 
70 waves 
e360 1,45 0.234 14,3 53 (500) 0.0593 Smooth 
E53 waves 
S51 1.60 0.258 16,0 56), 1 2.56 0.0723 Smooth 
.69 waves 
Gye 1.58 0.255 17.6 57 2.50 0.0856 Smooth 
265 waves 
10 1,16 0.150 30.7 170 1,35 0.0450 Smooth 
peo waves 
94 1.01 0.269 11.9 37.8 P02 0.0719 Smooth 
08 waves 
.03 1,05 0.185 SL 200 be Mo) 0.0403 Smooth 
L500 waves 
76 0.85 0.411 6.2 10.1 OLTZ 0.118 Smooth 
293 waves 
61 Smooth 
69 waves 
90 0.80 0.212 21.9 85 0.64 0.0726 Smooth 
00 waves 
.60 4.70 0.142 96 665 Ble 0.054 Smooth 
ea? lt ripples 
.70 3.74 Onres 765 7300 14,0 0.034 Smooth 
78 lg waves 
94 4.01 0.122 82 665 16,1 0.039 Smooth 
.09 lt waves 
19 23 0.281 24,4 87.5 10.4 0.087 Smooth 
27 ripples 
40 3216 0.277 23.4 S575 ot 0.083 ripples 
93 Smooth 
94 2.99 0,200 (Ag ey 148 8.93 0.058 Smooth 
03 mod ripples 
.30 3.00 0.200 28.8 142 9.00 0.058 Smooth 
atl ripples 
val 2.94 0.196 27.6 140 8.64 0.0556 Smooth 
7 sml ripples 
94 2.97 0.198 27.8 138 8.81 0.0567 Smooth 


.00 sml ripples 
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TABLE IV Data of Exper 


Sediment 
Size, 


Specific 
Gravity 


Chyn, 5.D. 3-I=I1V-1 


27 2 21.2 

27 10 0.847 2.65 22.2 13.6 
27 ll 21.8 

29 Kramer I 0.530 ana. * iz 7.4 
29 0 0.510 2.70 12 TA 
29 OI 0.550 2.70 ae 7a 
29 0.550 2.70 12 7.3 
29 0.550 2.70 12 7.7m 
27 Boragardi, Yen I-6-5 10.56 2.63 Lise 95.2 
27 7 925 95. 2m 
27 6-21 10.56 2.63 15.0 95.2 
27 22 15.0 
27 I-8-57 7ah2 2.64 19.0 64 

27 59 19.0 

27 8-61 Key ys 2.64 19.0 64 

27 3 19:2 

27 8-91 vin ¥3 2.64 19.8 64 

27 92 19.6 

27 Lived JY 5. I-L7f=26 4.4 2.66 26.5 49 

27 =27 26.5 

27 l7e-6 4.4 2.66 30.0 49 

27 ii 29.0 

rAd | Il-17i-8 3.4 2.66 28.5 42 

27 9 28.5, 

27 17j-5 3.4 2.66 27.5 42 

27 6 27.0 

27 OI-171-13 2.3 2.66 30.5 33 

ra bey 14-14 30.5 

27 17m-6 2.3 2.66 29.5 33 

27 7 29.5 

27 IX-17p-8 1.4 2.66 25.5 22 

27 9 25.5 

27 V=Live=9 3.6 2.66 24.5 aa 

27 10 24.5 ‘a 
27 VI-17w-12 a8 2.66 26.5 27.5 
27 is 26.5 

27 17x- 11 1.8 2.66 26:55 27.5 
27 12 26.5 

27 l?y-7 1.8 2.66 24.0 27.5 
27 8 24.0 
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on Ripple-Formation -- 


Av. Sheer 
Velocity 


Us 


cm/sec 


0.288 
ee 
0.300 
0.295 
0.320 
0.117 


0.127 


0.134 
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Con't. 


25.0 119 
25.0 120 
9.01 52.5 
9.00 50 
10.05 61 
LOat 35.5 
10.8 34.7 
905 6000 
1.110 6700 
590 4150 
647 4250 
650 4210 
294 2400 
255 2550 
192 1670 
192 1630 
128 950 
107 920 
52.5 347 
218 1670 
83.5 540 
WAC) 560 
69.3 520 


Tractive 


Force 
Tb 


Dynea/cm? 


pee at) ht 
(6s - é)d 


0. 


0. 


oocccococ 


0. 


0. 


0. 


0. 


On 


0. 


0. 


(OE 
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iy Bed 


0595 Smooth 
dune 0.5cm 
0566° Dune 0.5cm 
Smooth 


0534 
0405 
0472 
. 0563 
0661 
.0725 Frly Smooth 


scour 


.0857 Smooth 


Scoured 


0643 Smooth 


Scoured 


.0730 Frly Smooth 


Scoured 


.0756 Smooth 


rough 


0465 Smooth 


genripples 


0309 Smooth 


localripples 


.0395 Smooth 


lt local rips. 
0415 Smooth 
local ripples 
0516 Smooth 
lt local rips 
0379 Smooth 
local ripples 
0482 Smooth 
light dunes 
0511 Smooth 
local ripples 
0548 Smooth 
local ripples 


0471 Smooth 
dunes 
0422 Smooth 
local ripples 


te Ty, Peal 


1832-28 HY 5 October, 1958 


According to the writer, ripples can be formed by various kinds of disturb- : 
ances. It is only when the flow reaches the stage of instability that the ripples — 
can be maintained. Fig. 6 probably shows that ripples were caused by the 
side wall at the entrance. Therefore it does not contradict the writer’s ex- 
planation of ripple formation. 

It is interesting to notice that earlier Dr. Brooks(49) stated that the shear 
force or the shear velocity cannot be used in the study of sediment transport, 
because of the change of bed configuration. In his discussion of Fig. 9, he 
and Dr. Vanoni propose to use shear velocity in the dimensional analysis of 
ripple formation. In correspondence to the writer, Dr. Brooks stated that the 
shear velocity is a suitable parameter for studying the low rate of sediment 3 
transport up to dune regime, but it is not suitable for flow having flat bed (a 
transition from dune to antidune) or antidunes. 

The discussion by Messrs. Albertson, Simons and Richardson can be con- 
sidered as an excellent continuation of the writer’s paper on the mechanics 
of ripple-formation. By the use of Fig. A, most of the types of bed- } 
configuration can be estimated. Therefore, the classification of flow in 
alluvial channel has been clarified considerably. The space is not available — 
here to emphasize all the important contributions. A complete understanding ~ 
of Fig. A is recommended for those who are interested in the hydraulics of 
alluvial flow. The writer only brings up the points which need to be discussed — 
excluding those already mentioned earlier. i 

, 


In Fig. A, there are three kinds of abscissa used. Their relationship can 
be stated as follows: 


assume is die v {11) :} 

5 Us i 

Ag ese: :! 
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¥ 
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Considering that there is considerable information presented on Fig. A al- 
ready, for the sake of clarity we need not use different kinds of abscissa. 
Furthermore, the definition of the thickness of the laminar sublayer 6' is 
rather arbitrary, for example, Rouse has reported §' = 11.6 ¥/U,, 
Rotta(40) defines 5' = 6.83 ¥/U,, and according to Einstein and Li(42) that 
the laminar sublayer cannot be well defined. Furthermore, the ratio of sedi- 
ment to the thickness of the laminar sublayer does not have physical meaning 
in case of dune bed, therefore the parameter U,d/Y seems to be preferable. — 

In order to apply Fig. A for estimating the required shear for certain kinds 
of bed configuration, a trial and error method has to be employed. This step 
can be avoided if Fig. A is replotted accordingly U,/w against wd/Y as shown . 
in Fig. 14. The parameter wd/Y is also known to be the parameter pertaining 
to the straight line shown as for constant sediment size, if the fluid viscosity 
and shape factor of the sediment remain constant. 
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The criterion for the formation of antidunes is doubtful, because the 
Froude number is very important for antidune flow. Therefore, more labora- 
tory data on the formation of antidunes are needed. 

Through private correspondence with the writer, Dr. Bogardi made some 
very interesting comments on the discussion presented by Messrs. Albertson, 
Simmons, and Richardson. He thinks that the discussion completes the 
writer’s paper very well. With regard to Fig. A, he stated: 


“For constant temperature and constant d, the nature of regime is 
defined by the values of U,. 

If U, is assumed as the factor which quantitatively defines the con- 
figuration of the bed, the missing third parameter may be expected to 
depend uniquely upon the water depth D, the slope, as well as on the 
particle size d, or upon the combination of the above three quantities. 

A logic combination of these three quantities is the channel stability 
factor d/DS introduced in Hungary as early as 1942 to sediment investi- 
gation. This factor is proportional to ed/U,2, the inverse of the shear- 
velocity Froude number of the particle U,2/ed. 


Dr. Bogardi indicated that along the straight line of d = constant in Fig. A, the 
bed configuration changes as the channel stability factor changes. 

From Fig. A, and by assuming a temperature of 20° C, Dr. Bogardi has 
reduced the criterion of various bed configuration into a series of empirical 
equations in the form 


(35) 


U=Kdf 


as shown in Fig. 19, where o and # are determined empirically. 

The suggestions for future research made by Drs. Albertson, and Simons, 
and Mr. Richardson are mainly for studying the variation of bed-configuration. 
The writer is particularly interested in seeing item (5) being studied in 
the near future. Namely, the effect of the drag coefficient of the particle on 
the mechanics of sediment movement. It has been shown in the previous dis- 
cussion that the correlation of data by introducing the fall-velocity becomes 
less scatter and more logic, and from dimensional analysis, the effect of 
particle shape can be incorporated in the analysis if both the fall-velocity W 
and the intersity of the submerged weight of the particle AY d are used. 

_ Dr. Bogardi pointed out that the parameters used by the writer are similar 
to those proposed by Lane and Kalinski. By reviewing the works of Lane and 
Kalinski,(42) the writer could not find the similarity between the factor P and 
the shear-velocity Reynolds number U,d/ Y. The factor P has been defined by 
Lane and Kalinski as the ratio of the average sediment-concentration in the 
vertical to the concentration at the zero-level, and is a function of U,d/Y and 
n/d 1/6, where n is the Manning’s roughness factor. 

His suggestion to discuss various parameters related to the sediment 
transport is very inspiring. Those parameters related to the mechanics of 
ripple-formation have been discussed earlier in connection with dimensional 
analysis. A general discussion of various parameters is intended by the 
writer in connection with other phase of sediment transport at a later 
occasion, 
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A HIGH HEAD CAVITATION TEST STAND FOR HYDRAULIC TURBINES2@ 


Closure by W. G. Whippen and G. D. Johnson 


W. G. WHIPPEN! and G. D. JOHNSON,2 A.M. ASCE.—The writers appre- 
ciate the high quality of the several discussions, especially the significant 
differences in viewpoint that were raised. 

The repeatable accuracy of * 1/4% parallels the writers’ previous ex- 
perience with a well-designed pump-turbine test stand in that it also applies 
to a unit that has been completely dismantled and reinstalled at a later date 
» +... even with a completely new runner casting of identical design and 
finish. Although this cavitation stand does not lend itself to tests with con- 
trolled air content such as that at the California Institute of Technology, it is 
possible to determine whether or not the extremes of high or low air content 
have any measurable effect upon performance. 

Professor Hooper should note that the model of Figure 5 shows even 
higher efficiency than that of Figure 4, which had been increased some 25% 
in capacity without increasing overall dimensions. Since all instruments have 
been calibrated accurately in place or under installation conditions, it is felt 
that these results are accurate to within 1% or less. The difficulties in ob- 
taining accurate field test results have been pointed out and are freely admit- 
ted. The establishment of a reliable efficiency step-up formula will require 
equally accurate field tests, or the rational analysis suggested by the late 
Dr. Knapp. 

Mr. Vander Boegh will understand that, as a result of the commercial re- 
quirements for our testing facilities, it has not been convenient to test a com- 
pletely homologous (or identical) model in both the high and the low head 
cavitation stands. However, almost identical Kaplan units indicate an effi- 
ciency about 2% higher in the new stand. This is attributed to a poor venturi 
location and difficult alignment conditions in the older stand. 

_ It was stated in the paper that the primary function of the new stand is the 
testing of homologous models of prototype installations under prototype con- 
ditions. The high test heads and resulting high speeds require a special shaft 
that does not lend itself to the oil head and other equipment required for the 
Jetermination of Kaplan blade operating forces. This type of testing, as well 
aS visual observation, is carried out with sufficient accuracy in the older low 
nead stand. 

Results to date indicate that water temperature does not affect cavitation 
fest results as long as the corresponding vapor pressure is used in the 


4. Proc. Paper 1201, April, 1957, by W. G. Whippen and G. D. Johnson. 
l. Hydr. Engr. and Supervisor of Hydr. Lab., S. Morgan Smith Co., York, Pa. 
2. Chf. Hydr. Engr., S. Morgan Smith Co., York, Pa. 
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formula for sigma, but we must admit that commercial testing requirements ; 
have precluded any conclusive basic research in this as well as other areas 

of interest. The 100-ton cooling system permits the testing of models at the ~ 
water temperatures expected to prevail at the prototype installations. 

Although it was not emphasized in the paper, the head tank was sized to 
permit installation of a complete model inlet structure for units in semi- 
spiral case with short intakes; a modeled length of three or four diameters 
with a bell-mouth intake is installed upstream to reproduce the inflow condi- 
tions to a full spiral case. 

As mentioned in the paper, a separate pump-turbine test stand has been in 
operation for more than ten years. The construction of a new high-head 
pump-turbine stand has been authorized recently; upon its completion it will 
be possible to carry out test programs in both stands concurrently, if desired. 

Naturally, we were pleased by the late Dr. Knapp’s agreement that accu- 
rate laboratory tests on completely homologous models under prototype con- 
ditions could be substituted for less accurate, expensive and time-consuming ~ 
field tests. We agree with him that the model and the prototype drawings 
should be identical ..... except for size..... and that both the model and 
the prototype should be manufactured in accordance with the drawings ..... 
within reasonable tolerances, as is our usual practice. In view of the fact 
that a sincere attempt has been made to duplicate prototype conditions, the 
use of field testing techniques seems quite appropriate to us. As pointed out — 
in the paper, we are interested primarily in accuracy and simplicity; more 
elaborate instrumentation will be developed in due course. 

It does not seem unusual to us, in spite of the number of people who sus- 
pect otherwise, that we have found no significant changes in the calibrations 
of venturi meters over a period of years ..... provided that they are in ‘ 
first-class condition when calibrated and are maintained properly. 

The writers would like to mention that cavitation has two important mani- i 
festations: (1) rough operation accompanied by reduced efficiency and output — 
and (2) “pitting” or removal of metal from critical areas of the runner. 
Operation at critical sigma does not reduce efficiency and/or output signifi- 
cantly, but it could produce noticeable cavitation pitting damage over a 
period of time, unless the critical areas have been protected by stainless 
steel or some other resistant material. 

The writers hope that this paper, together with the points raised by the 
discussers, constitutes a useful addition to the technology on hydraulic 
laboratory practice. In closing, we would like to point out that the increasing © 
use of laboratory tests as contract acceptance tests indicates the desirability — 
of a universally accepted test code for models. ; 
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THE EFFICACY OF FLOOR SILLS UNDER DROWNED HYDRAULIC JUMPS 


Closure by Ahmed Shukry 


AHMED SHUKRY,! M. ASCE.—The writer wishes to thank all discussers 
of the paper for their contributary comments or adverse criticism. In fact, 
it was conceived that the paper would receive due discussion because the con- 
clusions drawn are opposite to some of the common-practice procedures for 
the following points: 


a- Laboratory scour tests may not give definite results when the jump 
formed below a river barrage is deeply submerged. 

b. The best location of an anti-scour sill is not at the floor edge, but rather 
upstream of that edge. 

c- The hydraulic advantages of dentated sills over solid ones, having the 
same peripheral dimensions are slight. 


Long before the present tests were carried out, points b and c were raised 
by the writer himself in a discussion to the results obtained by the Hydraulic 
Research Laboratory of the Delta Barrages, Egypt. This was in 1945 when 
the old Esna Barrage was under re-construction according to the final scheme 
shown in Fig. (I-f). Sand-scour tests on a model of this barrage showed that 
when the floor sill was placed at the middle section of the solid floor, the 
best results were obtained. This conclusion was surprising to engineers 
working on the design of regulators and barrages in Egypt. It had been the 
general practice to locate anti-scour sills at the floor edges, or even far 
downstream from those edges for more safety, as shown in Figs. (I-c) and 
(I-d). 

In 1947, the design features of the new Edfina Barrage were under investi- 
gation. Scour tests on a model of this barrage were also carried out in the 
Delta Barrages Laboratory. The results of these tests were similar to those 
of Esna Barrage regarding the better functioning of mid-floor sills. Guided 
by the fresh results of Esna and Edfina models, the sills of some of the minor 
regulators in Egypt were constructed at their mid-floor sections. This new 
practice which was starting to replace the old common one of using edge sills 
raised considerable discussion in 1948. The tests of Esna and Edfina barrage 
models were condemned although the two years of actual experience with the 
reconstructed Esna Barrage showed that with a midfloor sill, the downstream 
movable bed was practically scourless. Following the same testing procedure 
of sand-scour tests, which is explained in the paper, new tests were made on 
the same model of Edfina Barrage and in the same laboratory. The results of 


a. Proc. Paper 1260, June, 1957, by Ahmed Shukry. 
1. Prof. of Irrig., Faculty of Eng., Alexandria Univ., Alexandria, Egypt. 


. 


5 


1832-36 HY 5 October, 1958 


a 
these latter tests, which are presented in Fig. (2) in the paper, were more 
surprising due to their complete lack of agreement with the first ones. This — 
contradiction of the test results added more complexity to the problem, and 
the location of the sill in the prototype, which was under construction, was 
postponed until other tests, with other experimental technique, would be car- 
ried out in another laboratory. } 

The first thing the writer did was to review the results of most of the 
laboratory scour tests which had been performed on the Nile Barrages and 
major Egyptian regulators since 1928. It may be of interest to give a general 
description of the main features of this type of structures! Nile Barrages 
are open dams with vents varying between 5.0 and 8.0 m in width (16.40 and 
26.24 feet) and piers varying between 2.0 and 2.5 m in thickness (6.56 and 
8.20 feet). The working head on the barrage is generally about 4.5 m (14.76 
feet). The river bed material being a mixture of silt and sand, extending to 
considerable depths, a certain length of solid apron is always required down- 
stream of the piers to provide safety against percolation. The tailwater depth 
is much bigger than the conjugate depth corresponding to the flow issuing 
from the vents according to equation (1), and a deeply drowned jump is always 
formed on the downstream solid floor. The hydraulic conditions of these 
structures are, therefore, different from those of most of the American spill- 
way dams, where perfect jumps are usually formed on their downstream 
aprons, and for which the Froude number is the main influencing parameter 
of the flow. The adhering property of the drowned jump is that it does not 
dissipate as much energy as that of a perfect jump. The more the jump is 
drowned, the more will be the part of the excess energy which is to be dissi- _ 
pated by viscous action downstream of the jump. Accordingly, the Reynolds 
parameter may become as much important as that of Froude. This explains 
the inconsistency of some of the previous scour results as pointed out in the 
paper. 

From the discussions by Mr. Peterka and Dr. Leliavsky, it seems that 
there is a general misconception of the conclusion given in the paper regard- 
ing the validity of sand scour tests. It is only when the jump is deeply 
drowned that the results of these tests may become uncertain. The writer 
wishes to write the suggestion of Dr. Leliavsky as follows: For important 
structures operated under deeply drowned jumps, velocity tests should be 
used to check any uncertainty in the results of scour tests. In fact, this pro- — 
cedure was followed in the design of Edfina Barrage. However, since the ob- © 
ject of the velocity tests is to predict the erosive tendency of water, the tests — 
should be made on models with fixed beds. Velocity patterns recordedina 
movable-bed model may not indicate the actual erosive forces in water be- _ 
cause part of the excess energy may have been dissipated in forming the 
scour pool. 

When the jump is perfect or slightly submerged, there will be only one 
effective parameter. In this case qualitative similarity of scour tests is 
feasible. Professor Vito A. Vanoni and the writer have recently used scour 
tests to design the optimum proportions of some of the bed stabilizers re- 
quired for some of the flood control channels in Southern California. The re- 
port of these tests is not yet published. 

Mr. Thomas raised a question about the difference in action of an end sill _ 
when the jump is drowned and when the jump is clear, provided that in both — 
cases measures are taken to normalize the flow upstream of the floor end. 
The writer believes that if such normalization has taken place by the provisi 


\SCE DISCUSSION 1832-37 


f sills, baffle piers or blocks, there will be no need for the addition of an 
-dge sill in both cases of jumps. Such an addition will not secure more 

safety to the floor, but on the contrary, an edge sill may cause re-disturbance 
f the normalized velocity distribution and may give rise to bed shearing 
orces due to the momentum changes at the floor edge. This is clear in Figs. 
IIa) and (IIb) of the paper. 

The statement given by the writer that “The efficiency of any sill against 
ed erosion is indicated by the rate of adjustment of the flow to the normal 
listribution in the downstream channel,” is not an assumption as commented 
y Mr. Peterka but it is rather a definition of the efficiency of a floor sill. 
This sill is intended to deflect the high-velocity currents upwards in order to 
sstablish normal distribution of velocities within the shortest length from the 
lier edges. 

The scour hole is mathematically defined in equation (2) of the paper just 
0 illustrate the various parameters involved in the problem of drowned 
umps, under which conditions qualitative similarity between models and pro- 
otypes may not be feasible. The writer agrees with Mr. Peterka that identi- 
al erosion patterns may be obtained for the prototype, if models of the same 
structures are constructed with different scales, and if the same bed material 
Ss used in all models. Such testing procedure may be, however, very long and 
‘xpensive. 

The results of the present tests show that the improvements gained by the 
low when using toothed sills instead of solid ones of the same peripheral 
limensions, are relatively small compared with the expenditure of maintain - 
ng toothed sills against the more severe cavitation actions. 

The writer agrees with Professor Ingersoll that due to the part of the ex- 
ess energy which is dissipated by viscous action downstream of the model, 


lifferences between the model and the prototype values of B2 may be expected. 


\ccordingly, the values of To, as computed by equation (6), may be less for 
he model than for the prototype. However, due to the relatively short length 
f the solid apron over which the normal velocity distribution is established, 
he expected differences in the values of To may be small. 

The application of the energy equation between Stations I and VII for the 
onditions presented in Figs. (6a), (7a) and (7b), as given by Professor 
ngersoll, is an illuminating example which illustrates numerically the wake 
ffects of the piers. Due to this wake, the energy losses or the shear 
tresses at the bed cannot be numerically determined between successive 
ections, unless the distribution of velocities in the lateral direction has been 
Iso recorded. In the present work, qualitative analyses of various designs 
yas only required. This was achieved by velocity measurements along the 
lume centerline. 

From his experience in West Pakistan on works similar to Egyptian bar- 
ages and regulators, Mr. Mushtaq Ahmad advocates that velocity distribu- 
ion and comparative studies of scour should be taken as criteria for the 
roper design of scour prevention devices. Mr. Ahmad states the general 
onclusions of his very interesting experimental study on inclined weirs of 
ne “Glacis” type. Fig. (I-c) in Mr. Ahmad discussion shows that for a deeply 
rowned jump, the depth of the scour pool is relatively larger than what the 
elocity distribution diagrams might indicate. This may explain the writer’s 
Oint of view that velocity distribution studies should be carried out on models 
ith fixed beds. Velocity distribution, recorded after the excess energy in 
rater has been expired in scour, gives no indication of what might take place 
1 nature. 
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The discussion by Mr. Ahmad under the heading “Mechanics of Energy 
Dissipation in the Stilling Basin and Function of Scour Control Devices” is 
most interesting. It throws much light on the hydraulic jump as an energy 
converter. 

Questions were raised by Mr. Thomas and Mr. Peterka about the data to 
which Figs. (8) to (II) relate. In all tests, the discharge was passed below the - 
lower gate. The super-critical depth of the flow between the piers was 
1.20 m (3.94 ft.), while the tailwater depth was 4.50 m (14.76 ft.). The bar- 
rage was assumed to work under a hypothetical head of 4.00 m (13.12 ft.). 

The maximum rate of discharge through the three barrage vents was fixed 
from the flood conditions and was taken as 302 cubic m per sec. (10670 cubic © 
ft. per sec.). The width of one vent is 8.00 m (26.24 ft.) and the thickness of 
one pier is 2.50 m (8.20 ft.). The length of the downstream solid floor, ex- 
cluding the block extension, is 20.00 m (65.62 ft.). All these dimensions re- ¢ 
late to the prototype. F 

In closing, the writer wishes to mention that although research on erosion ~ 
problems have started earlier than most of the works in other problems of 
hydraulic engineering, there are still many aspects in the scour phenomenon 
requiring further analyses. Rationalization of the old empirical rules in the 7 
light of modern fluid mechanics seems to be still lacking. ; 


a 


\ 
CORRECTION. —Page 6 in the paper, conversion of 302 cubic ms. per sec. 
yields 10670 cubic ft. per sec., rather than 106700 cubic ft./sec. é 
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The fifth reason presented has to do with the better dissipating action of a 
slightly drowned jump. While this is true, the point is somewhat academic 
when applied to a large flood control or multi-purpose reservoir project with 
basin designed for a discharge three times or so larger than any previously © 
recorded flood. 

With regard to abrasion damage caused by bed material or other debris 
circulating on the apron, the writer agrees with Mr. Peterka that “dead” 
areas should be eliminated and that ground rollers can sometimes be too 
strong. 

Mr. Peterka’s theory concerning the cause of damage at Lucky Peak Dam 
is an interesting one. Whether the damage was caused by impinging jets or 
vortices (the model did not indicate either), the writer agrees that further in 
vestigation is desirable. 

In conclusion, it was hoped that this paper would contribute to the profes- 
sion a little more knowledge regarding the difficulties and shortcomings of 
actual prototype hydraulic structures and to some extent the possible limita 
tions of hydraulic models. It was further hoped that the paper would serve t 
open discussion of experiences of others in this field. 
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STILLING BASIN EXPERIENCES OF THE CORPS OF ENGINEERS2@ 


Closure by R. H. Berryhill 


R. H. BERRYHILL,! A.M. ASCE.—Since this paper was presented, a rare 
pportunity to observe stilling basins in action presented itself in the form of 
1e Southwest Flood of 1957. This flood ranks as one of the six most destruc - 
ve floods in America since 1900. At least 30 dams of the Corps of Engi- 
eers were affected to varying degrees in Texas, Oklahoma, Arkansas, 
ansas, and Louisiana. Many of the spillways in this area were operated for 
e first time during this flood. The experiences encountered here and in 
everal other parts of the country recently should probably be told in a sub- 
equent paper. 

With regard to the discussion by Mr. Peterka, which incidentally was very 
uch appreciated by the author, comments in the following paragraphs are 
fered. 

The first reason given by Mr. Peterka for believing that a basin floor 
10uld be set at full conjugate depth or more is that the safety of the entire 
sructure may hinge on whether the baffle piers remain fully effective after ' 
drasion or erosion. It should be noted that the model tests conducted by the 
orps for hydraulic-pump-type basins invariably include tests with the baffle 
ers removed to insure that spray action does not occur even with the ex- 
‘emely rare flow termed the spillway design discharge. An additional safety 
ctor lies in the fact that even on those projects where some rounding and 
jrading of the baffle piers has occurred over the years, the resulting stilling 
fect of the baffle piers certainly could not be termed “ineffective,” although 
is granted that some reduction in their efficiency had taken place. Neither 
‘the two examples given by Mr. Peterka were the result of reliance on baf- 
e piers. 

The second and third reasons given by Mr. Peterka concerns the future 
*gradation of the river channel and lowering of the tailwater as well as the 
yssible incorrectness of the tailwater curve. It should be noted that for 
ajor projects, considerable investigation precedes the determination of the 
ilwater curves, both present and ultimate degraded. A safety factor is in- 
uded in the determination of the latter, and the stilling basin design is made 
lequate for both tailwater conditions. 

The fourth reason presented has to do with the lag of actual tailwater 
pth during periods of increasing discharge. With the usual large flood- 
mtrol project, this feature is usually of little consequence whereas on a 
in-of-river project it is sometimes of major importance. Discretion in this 
atter during the planning stage of such projects is absolutely necessary. 


Proc. Paper 1264, June, 1957, by R. H. Berryhill. 
Chief, Design Branch, U. S. Army Engineer District, Albuquerque, N. Mex. 
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The fifth reason presented has to do with the better dissipating action of a 
slightly drowned jump. While this is true, the point is somewhat academic 
when applied to a large flood control or multi-purpose reservoir project with 
basin designed for a discharge three times or so larger than any previously 
recorded flood. 

With regard to abrasion damage caused by bed material or other debris 
circulating on the apron, the writer agrees with Mr. Peterka that “dead” 
areas should be eliminated and that ground rollers can sometimes be too 
strong. ; 

Mr. Peterka’s theory concerning the cause of damage at Lucky Peak D 
is an interesting one. Whether the damage was caused by impinging jets or 
vortices (the model did not indicate either), the writer agrees that further in- 
vestigation is desirable. 

In conclusion, it was hoped that this paper would contribute to the profes- 
sion a little more knowledge regarding the difficulties and shortcomings of 
actual prototype hydraulic structures and to some extent the possible limita- 
tions of hydraulic models. It was further hoped that the paper would serve to 
open discussion of experiences of others in this field. 
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A STUDY OF BUCKET-TYPE ENERGY DISSIPATOR CHARACTERISTICS@ 


Closure by M. B. McPherson and M. H. Karr 


M. B. McPHERSON, |! A.M. ASCE and M. H. KARR,2 J.M. ASCE. —After the 
paper had been completed it was realized too late for rectification, that all 
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a. Proc. Paper 1266, June, 1957, by M. B. McPherson and M. H. Karr. 
1. Research Engr., Philadelphia Water Department, Philadelphia, Penna. 
2. Investigations Engr., State Water Resources Bd., Salem, Oregon. 
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test runs had been made with a falling tailwater. Thereupon a special series 
of tests was conducted using the Greensboro model shown in Figure 2, with a 
horizontal solid floor at the elevation of the bucket invert. The results of 
these tests, for a q-parameter of 0.055 and an h;/R of 3.4, are plotted in 

Fig. C-1. The free-jet definitely occurred at an hp less than 0.2 hg with a 
falling tailwater, but with a rising tailwater roller action was not unequivoc- 
ably established until hp was about 0.3 hg. 

Plotted in Figure C -2 is the ratio of the hg’s for the free-jet limit with a 
rising and falling tailwater versus h1/R for the tests of Fig. C-1, comparable 
tests on the Sakuma Dam model! and for unpublished design data from the 
Army Engineer’s Waterways Experiment Station. The data points are in- 
dorsed with the corresponding q-parameter values. The comparative tail- 
water requirements appear to be a function only of hj/R, but this limited evi- 
dence is hardly conclusive. 

The statement given under “Recommendations for design”: - - “Free-jet 
performance can be expected whenever hp is less than 0.2 hg, in general,” is 
therefore applicable to a falling tailwater only. Due to the small changes in 
h2/h1 at an hp less than about 0.3 hg, and considering the limited evidence 
presented in Figure C-2, a limit of hp equals 0.3 or slightly higher is indi- 
cated for hj/R greater than about 3 to marginally avoid a free-jet with a ris- 
- ing tailwater. This tentative conclusion diminishes somewhat the competitive 
status claimed for the solid bucket over a hydraulic jump in a horizontal 
channel (See Figure 5). 

The data presented by Mr. Dougherty from the Weiss Dam model studies 
do not and should not agree with the graph of Fig. 3 for several reasons. The 
head above the crest is as much as 41.5% of hj. Mr. Dougherty recognizes 
the probable effect on the data by the resultant roller distortion in stating 
that “This is attributed to the effect of the flow over the crest.” Because of 
the spillway crest piers the flow width at the crest is only 83.5% of the bucket 
width. Only the exit angle of 45° should be compared with Fig. 3. (The data 
for the modified 45° exit is not identified; no correlation for a non-radial exit 
should be expected). The q-parameter curves for 0.12 and 0.15 are Mr. 
Dougherty’s extrapolations. Whether or not they are representative is diffi- 
cult to say. The three curves from Fig. 3 in his graph have been accurately 
reproduced, Since the hj/R ratios for the 45° exit represented a good range 
(1.9, 3.1 and 5.1), it is unfortunate that the Weiss model characteristics were 
so dissimilar as to preclude a proper comparison with the data by the authors. 

It is not clear why Mr. Dougherty is hesitant to accept locating the down- 
stream bed at the elevation of the bucket invert. The only reason the sloping 
floor tests represented in Fig. 2 were given was to show that despite sloping 
the bed beyond a point where the surge roller action is dissipated, the rela- 


fr 


tion of hg in terms of hp is unaffected even with moderate modifications of the — 


bed adjacent to the bucket. In Fig. C-1 is shown a plotted point for an hj/R 

of 2.3 (25'R rather than 17'R of all other points) for which the 41' horizontal 
portion of the solid floor of Fig. 2 was reduced to 32' through replacement 
with the larger bucket in the original model tests. The slight displacement of 
this data point could be attributed to a small degree of crowding of the expand- 
ing surge by the nearer sloping bed. The data for the three cases detailed in 


1, “Hydraulic Model Tests on Flood Spillway of Sakuma Dam,” by Atsuya 
Okada and Takeshi Ishibashi, Central Research Institute, Electric Power 
Industry, Japan, June 1, 1956. 
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Fig. 2 are in good agreement with the solid and movable bed horizontal floor 
data of Fig. C-1. 

In the Sakuma Dam model tests! it was necessary to provide a substantial 
secondary dam downstream from the bucket to insure adequate tailwater 
depths. When this dam was placed too close to the bucket, the surge expan- 
sion and the general action was disrupted considerably. For a certain criti- 
cal position, and beyond, the action above and at the bed was unaffected 
and unchanged. 

Mr. Dougherty’s statements concerning erosion are entirely true. How- 
ever, in the low to medium head flood control and water supply spillways for 
which the study data is particularly suited, scouring below the main structure i 
can often be tolerated. The worst scour would occur at the design head, 
which is of course rarely attained. To answer the question on bed material 
size, the 1:40 scale Greensboro model bed material was “No. 1-B aggregate, 
light on No. 8 with an analysis as follows: 100% passing a one-half inch sieve, 
75-100% passing a 3/8", 10-30% passing a No. 4 sieve and less than 1% pass- 
ing a No. 8 sieve.” It is quite likely that prototype erosion would be more ex- © 
tensive with a continuous horizontal bed than with a downstream sloping bed 
similar to the one shown in Fig. 2. ' 

In partial answer to Mr. Dougherty’s question concerning setting wall 
heights lower than hg, the training walls of both the Greensboro and Penn 4 
Forest Dams were set below hg with a rip-rapped berm immediately above — 
the walls in the splash area. How much the walls can be lowered is related 
to the position of the end of the wall relative to the surge location, and to the ~ 
bed and foundation conditions beyond a longitudinal wall in the case of a short 4 
training wall. A training wall must be at least high enough to contain hp. It 
appears doubtful if training walls could be lowered much more than 2/3 to i 
4/5 of hg without causing transverse surge action with a consequent disrup- y 
tion of the bucket roller. é 

The Sakuma Dam 1:90 scale model tests! were performed with a maximum 
q-parameter of around 0.03 and an hj/R of between about 6.0 to 6.7 (using the - 
net hj). Only one series of runs were made with an R of 17 1/2" but 46 series” 
were performed with an R of 8 3/4". Unfortunately, because of a failure to 
detail the many tests, it is not possible to correlate test results with test 
conditions except in a few isolated instances. The authors, in consideration 
of the above and after careful review, have failed to find any sound basis for 
the interpretation of Mr. Elevatorski that “Investigations by Okada and 
Ishibashi substantiate the functional relationship given by (Mr. Elevatorski’s) 
equation (3)---.” The authors could not locate the quotation alluded to in the 
remainder of the sentence. The authors have failed repeatedly in attempts to 
identify those overall test results which were for a 450 radial exit without 
extensions or modifications of the lip. The approach slope was constant at 
about (10 on 8). 

The U. S. B. R. has made tests on a slotted type of bucket! as shown in 
Fig. C-3, a reproduction of a page from their report. The report acknowl- 
edges fully the difference between the solid type and slotted type bucket: 
“The hydraulic action and the resulting performance of the two buckets are 
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1, Bureau of Reclamation, Progress Report III, Research Study on Stilling 
Basins, Energy Dissipators, and Associated Appurtenances, Section 7, 
Slotted and Solid Buckets for High, Medium and Low Dam Spillways, Hy- 
draulic Laboratory Report No. Hyd.- 415, Denver, July 1, 1956. 
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A. GRAND 


See Figure 5 for 
tooth detail 


B. ANGOSTURA TYPE SLOTTED BUCKET 


FIG. C-3: REPRODUCTION OF FIG.I, PAGE 3, FROM 
US.BR. HYDRAULIC LABORATORY REPORT 
NO. HYD.- 415, JULY |, 1956. 
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quite different.” All the data presented appears to be for a falling tailwater 
according to the manner in which the results were tabulated. Test results ; 
reported include only “seeepout” conditions (same as free-jet), and “diving 
flow” conditions, with no data for any conditions in-between. In the slotted 
bucket tests, at a higher tailwater, a point is reached where the jet leaving 

the bucket is depressed downwards with consequent severe scouring action; 
this phenomenon is probably impossible of occurrence with a circular solid 
bucket. The values of hp are not given for the tests. In order to compare 
performance, values of hg/h, and the q-parameter were computed-from their — 
tabulations of test results and plotted in Fig. C-4, using the headwater pool 
waterlevel to bucket invert as hy (gross head). Also shown are lines of hp in — 
terms of hg taken from Fig. 3, for comparison. Apparently the slotted bucket 
requires a smaller hg to form a free-fet with a falling tailwater. The “diving 
flow” regime appears to be a serious restriction. The ranges of variables 
tested with the slotted bucket appear in Table C-1. if 

In the U.S.B.R. report are mentioned (no test data) tests with three solid 
buckets with a conventional radial 45° exit. The hj/R tested was only about — 
1.5, 2.4 and 3.6 at a maximum q-parameter of around 0.10 for which the head 
on the crest was about 27% of hj. 5 

Mr. Elevatorski’s equations (3) and (6) are modifications of the relations 
set forth on pp. 39-45 in the U.S.B.R. report. The report relations are inte 7 
pretations based solely upon the somewhat limited tests with a slotted type 
bucket, which appear to be inconclusive in many respects. 

It is difficult to see why Mr. Elevatorski would prefer to interpret the re- 
quired minimum tailwater and minimum bucket radius for a solid bucket on 4] 
the strength of results from radically different slotted bucket tests. More- — 
over, the ranges of variables represented in Table 1 and Table 2 are much r 
more comprehensive than those in Table C-1. The authors feel that the mini- 
mum bucket radius requirements for a solid bucket have been fairly well de- 
fined through their tests. The solid bucket tests mentioned in the U.S.B.R. 
report were obviously too fragmentary to permit any generalized comparison - 
with those for a slotted bucket. ‘ 

Mr. Elevatorski’s interpretation of the U.S.B.R. report in stating that the — 
slotted bucket approach slope does not, “within ordinary limits, affect the § 
bucket size” is unfounded since only one approach slope was used in the tests 

Table 1 provided by Mr. Elevatorski is a welcome contribution, but is not 
by any means a complete inventory. Two of the Greensboro Dam figures " 
given in his table are incorrect: the maximum discharge per foot of bucket 
dataset 137 cfis/ft. and the maximum gross head is 58 feet, rather than 540 
and 72. 

The discussions have contributed one very important point: Comparison 
with the authors’ results cannot be made unless the comparative tests are 
for a solid bucket with a radial exit (without exit modifications), flow uni- 
formly distributed over the channel width and with the height from the spill- 
way crest to the bucket invert being more than about 75-80% of hy. Modifica 
tions of the exit may have profound effects upon performance. 

Acknowledgement is made to Mr. Fred R. Brown for providing the author. 
with the W. E. S. data, to Mr. Elevatorski for a copy of reference d and to 
Mr. S. H. Poe of the U.S.B.R. for a copy of reference e. Also, Mr. 
Elevatorski discovered a serious omission in an explanation given in the 
original draft, in the form in which it was presented in the ASCE Pittsburgh 
Convention on October 18, 1956. 
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TABLE C-1 


IGES OF VARIABLES INVESTIGATED WITH AN ANGOSTURA TYPE SLOTTED 
“KET WITH (10 ON 7) ENTRANCE AND 45° MODIFIED EXIT. 


sing gross head for h))- U-S.B.R- REPORT®. (See Figure C-4). — 


Readius Range Range Range of Range 
R, of h,, of q, q of 
Feet Feet efs/ft. Ve hy / R 
0.50 5.20 0.31 0.0046 10.4 
to to to to 

5.68 2.25 0.029 11.4 

0.75 5.42 1.05 0.015 ae 
to to to 

5.88 3-52 0.043 7:8 

1.00 5.54 1.58 0.022 5.5 
to to to to 

6.22 6.08 0.069 6.2 

1.50 5.63 2.00 0.026 3.8 
to to to to 


(h,VS.h, CURVES ARE FROM FIG.3, PROC. PAPER 1266 ) 
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JULY 1, 1956- 
SLOTTED BUCKET 


SLOTTED BUCKET VS. SOLID BUCKET 
CHARACTERISTICS 


FIGURE C- 4. 
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SYSTEMATIC CHANGES IN THE BEDS OF ALLUVIAL RIVERS@ 


Closure by Walter C. Carey and M. Dean Keller 


WALTER C. CAREY,! M. ASCE and M. DEAN KELLER,?2 J.M. ASCE.—The 
vriters are gratified by the interest aroused by their paper. They are also 
leased with the nature of the comments by Messrs. Blench, Vanoni, and Liu 
ind the evidence that they have reported phenomena in a field relatively unex- 
lored and one that may have an important bearing on river hydraulics. The 
vriters hope that it will be possible to explore this field further, in a more 
systematic way. The remarks below deal with each of the discussions in turn: 


T. Blench.—The observed rate of travel of bed waves of 250 feet per day is 
lighly significant as a rough indication of the possible quantity of material 
hat may be moved in this way. It may be possible to initiate, on the lower 
Mississippi, a set of sand wave observations conformable with the require- 
nents set forth by Mr. Blench. 


Vito A. Vanoni.—As to the increase in size of dunes with increase in stage, 
he writers were vaguely aware that “stage” was more important than 
velocity.” Although unable to explain why, they think it reasonable, and by 
he same token, they think it reasonable that the “increased velocity” in Mr. 
Janoni’s flume, without any increase in stage, resulted in decreased size of 
lunes. The writers agree on the desirability of similar sand wave observa- 
ions at flood stage; in fact, “all around the calendar” and through all varia- 
ions in stage. Such observations might also throw light on what causes the 
‘ating curve at a given station to be a “loop.” The writers think it possible 
hat the sediment actually forming the waves is not measureable by any exist- 
ng procedure. All other sediment, they suspect, may not be too important in 
‘overning the hydraulic properties of a stream. They would not be surprised 
o find that “waves” at different stages may be constituted of somewhat dif- 
erent materials (as to size). 


Mr. Hsin-Kuan Liu. —Concerning the transport of bed load around a bend, 
he writers are prepared to find (when the matter is studied) that bends have 
igher transporting power than crossings; in fact, that may be why there are 
crossings.” As to the remark “.. . bed configuration is not an independent : 
ariable.... ,” the writers suspect it is controlled by the river itself for 
‘easons that may be understood some day as a result of further study. This 
oncept includes all systematic variations in the river bed with stage, not on- 
y as to systems of sand waves, and heights of crossing, but possibly the 


. Proc. Paper 1331, August, 1957, by Walter C. Carey and M. Dean Keller. 

. Asst. to Chf. of Eng. Div., New Orleans Dist., Corps of Engrs., U. S. Dept. 
of The Army, New Orleans, La. 

. Structural Engr., Bedell & Nelson, Cons. Engrs., New Orleans, La. 
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nature of the material on the bed surface at a given moment. With respect to 
Mr. Lane’s thoughts as to the influence of the sea on slopes in the reach un- — 
der discussion, the writers’ observations still further up the river revealed 
similar phenomena. Furthermore, the systems of waves below Baton Rouge, 
where the river runs through a deep trough (sometimes in the Pleistocene) 
did not apparently differ from those above Baton Rouge, where the river is 
pretty much the same for 600 or 700 miles. 

The writers were pleased to note that all of the discussions point to the 
need for more study of systematic river bed variations with respect to their 
possible effects on water surface slopes and flow conditions. 

An added point, the writers wish to report a recent technological advance 
that makes it easier to profitably study and compare Fathometer rolls, for 
such studies as those under discussion: Ordinarily the Fathometer roll is 
most inconvenient to handle and, being opaque, cannot be compared conveni- 
ently with another roll, even when to the same scale. A method has been 
developed for continuously reproducing Fathometer rolls on translucent vel- 
lum by the so-called “electrostatic” process. The reproduced rolls can re 
ly be superimposed on one another, or traced, or used for making 
reproductions. 
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CHARACTERISTICS OF FLOW OVER TERMINAL WEIRS AND SILLS@ 


Discussion by F. Paderi 
Closure by P. K. Kandaswamy and Hunter Rouse 


F. PADERI.!—The writer has read with pleasure this paper which demon- 
strates such explicit recognition of the persistent theoretical and practical 
nterest in the subject of Basin weirs. He also appreciates the comprehen- 
sive portrayal of the problem for the entire range 


bee | y 
h+w 


ind he finds of particular merit the care and breadth of view with which the 
iuthors revise and refine their own prior results to approximate more close- 
y the true phenomenon. 

Along these same lines, in many of his own works from as early as 
937, (1,2) the writer has studied this important phase of hydraulics both 
inalytically and experimentally. These works are also to be considered as 
steps toward a progressive deepening of knowledge on the subject. 

One such item of accomplishment was contained in the paper(3) presented 
o the Sixth Congress of the International Association for Hydraulic Research 
it The Hague in September 1955. This paper placed the fundamental contri- 
yutions of Boss, Rouse and Ippen(4) in evidence and indicated the conservative 
ature of the limitations, for practical purposes, of the writer’s non- 
limensional expression for the discharge coefficient 


2 
sed Ontos 0.302 1- (4 ny (2) 
Maw y 
n the relationship 
ray 3/; z 
+ = ; (t Ls 4 h 
or the calculation of the capacity of overflow dams over the whole of range 
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. Proc. Paper 1345, August, 1957, by P. K. Kandaswamy and Hunter Rouse. 
. Ass’t. Prof. of Hydraulics, Univ. of Pisa, Pisa, Italy. 
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As was pointed out in the foregoing paper, the use of Eq. (2) in the case of 
overflow dams of the proper profile is justified by experiments on Basin 
weirs which confirm the existence of such a trend. 

The numerous experiments of Meyer and See(5,6) and those more recenay Z 
performed in the hydraulics laboratory of the U. S. Bureau of Reclamation(7 
should also be cited. 

In his studies on the Basin type of weir, the writer started with a deter- ; 
mination of the merits and limits of validity —as recognized by the authors—of — 
the empirical expression of Rehbock for Cg and of the theoretical treatment 
of von Mises, (8,9) and he warned about the controlling effect of the phenomena ~ 
of the second critical stage. : 

The writer dealt, in effect, with a homogeneous perfect fluid starting from 
rest under the sole action of gravity (and hence in irrotational motion). He 
applied the theorem of Bernoulli and also—after Boussinesq, (10,11) who found 
the constant value 2Cg/3 = 0.4342, and De Marchi(12) —the theorem of mo- 
mentum, and arrived at an analytic expression for the discharge coefficient 
Cqg of the Basin type of weir. 

Following the inverse procedure, he then began with the theoretical and 
experimental condition of the second critical stage (Cg = 3/2 Y2 = 3 x 0.707/2 
= 1.06) and, without experimental aid, determined the value Cg = 3 x 0.4048/2 
= 0.607 at the other extreme h/(h + w) = 0 of range (1), as was then expressed 
by Eq. (2). The distribution of velocity at the section of measurement was 
assumed practically uniform. 

The writer does not consider it worthwhile to dwell upon details of the 
subject; however, he does believe it to be of general interest to note that, 
starting with the study of 1937, these results were obtained by a procedure 
analogous to that which he has presented in his discussion of Proceedings 
Paper 1038.(13) This procedure consists of the introduction, at a certain 
point in the analytical treatment, of a coefficient K (of the class Kg), which 
resulted analytically in a value equal to 2 for the theoretical and experimental 
conditions of the second critical stage. P 

The last word still remains to be said regarding the experimental aspect 
of this problem; the writer—after Rouse—has long concerned himself with 
this aspect of the free overfall for both horizontal and sloping channels.(13,14 

Laboratory practice warns against experiments at too small heads, be- 
cause of related phenomena which complicate the laws of similitude. Perhap 
certain phenomena are also brought into existence by small heights of weir. 
These phenomena could possibly be connected with the corresponding distri- — 
bution of velocity. 

One particular phenomenon has been perceived by the writer in the course 
of his long experimentation on the free overfall: the accentuation of the last 
undulation just before the final drop in surface profile as subcritical flow ap- 
proaches the brink (channel of mild slope). This is related to local disturb- 
ances of the flow at the inlet of the test channel. Consequently, the theoreti- — 
cal critical section of depth _ 


is shifted slightly away from the brink, and hence the absolute value of x/h 
becomes greater than x,/Yc. 
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The writer observes, however, that the second critical stage has a par- 
cular constant physical characteristic of stability. This limits the conse- 
uences of such disturbances. It therefore becomes feasible, because of the 
mited slope of the free surface in the critical region, to use this character- 
stic to realize an “overfall-brink discharge meter.” The corresponding 
xperiments are in progress. 15) 

That all these fields of research offer great possibilities of development 
5 attested by the authors’ indications, which are as valuable from the ex- 
erimental side as for their theoretical repercussions. 


P. K. KANDASWAMY,! A.M. ASCE and HUNTER ROUSE,2 M. ASCE.—Dis- 
ussions of the writers’ presentation by Messrs. Engel and Paderi are much 
ppreciated for the sidelights that they throw upon the presently incomplete 
nderstanding of the phenomenon under consideration. Mr. Engel states that 
1e writers have not explained when and why it occurs; and Mr. Paderi, in his 
wn work on the problem, finds rather different results. Hence the question 
aised seems to involve not only when and why but also whether. 

Mr. Engel! prefers the earlier empirical idea of a base coefficient (which 
e believes to decrease as the head-depth ratio increases) and an apparently 
idependent velocity -of-approach coefficient. As explained in the paper, both 
1e contraction coefficient and the velocity -of-approach coefficient (which to- 
ether determine the discharge coefficient) are functions of, and increase 
ith, the head-depth ratio throughout the weir range. Mr. Paderi, proceeding 
n the latter basis, has derived a discharge relationship that increases con- 
nuously to the free-overfall limit. Although he has no experimental data to 
onfirm this relationship in the vicinity of its upper limit, several series of 
ata in the lower range duplicate its trend. As Mr. Paderi has pointed out, 
owever, other series of data exist which duplicate the trend of the writers’ 
ndings. 

If one approaches the zone in question from the other direction, it must be 
oted that the original observation by BOss—that a very low sill can be placed 
t the brink of a free overfall without appreciably increasing the depth of the 
pproaching flow—is readily verified. It follows that the discharge coefficient 
just attain a maximum value at some intermediate sill height. Evidence of 
uch a reversal in trend can also be seen in the basic difference in nappe- 
rofile variation for the weir and sill ranges (compare Figs. 4 and 5). 
fhether the change in function occurs and when it occurs thus seem to be 
nswered satisfactorily by the writers’ diagrams. Why it occurs is probably 
ubject to the same sort of answer as why a sluice gate has the effect that it 
oes as it begins to penetrate the surface of flowing water; such gravitational 
henomena are more readily subject to analytical evaluation (for example, by 
elaxation) than to descriptive explanation. 

Mr. Paderi has pointed out still another aspect of wave formation in its 
elation to the profile measurements, and both he and Mr. Engel have men- 
oned the possible effect of surface resistance. Each of these factors un- 
subtedly has bearing upon the discrepancies that are evident not only between 
\e independent measurements of the two writers but to a much greater extent 
mong sets of data from wholly unrelated investigations. As the writers 
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showed, the two influences can be minimized—but never wholly eliminated— 
by using as large a radius of floor curvature and as short a subsequent re- 

gion of uniformity as possible, and by measuring the head consistently only 

sufficiently far upstream to avoid the final zone of curvilinearity. 

Mr. Engel appears to disregard the effect of curvilinearity upon the 
critical-depth relationships developed for parallel flow, since he uses these 
to relate the depth over the sharp-crested weir to that over the broad- 
crested weir and to the head upstream. He also raises questions as to the 
magnitudes of the actual depth ratios, which—if the writers understand his — 
query —can all be read directly from the several diagrams that they provided 
The critical depth for curvilinear flow is a matter that still remains open to 
generalization; the writers will hence be interested in following Mr. Paderi’ 
development of an overfall-brink discharge meter. 
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TURBULENCE IN A DIFFUSER BOUNDARY LAYER@ 


Closure by J. M. Robertson and G. L. Calehuff 


J.M. ROBERTSON, M. ASCE and G. L. CALEHUFF.2—This paper was 
ntentionally kept reasonably short and as uncomplicated as possible in the 
lope that a larger number of people could follow and understand it. The basic 
-elations were therefore presented in the simplest fashion. As Mr. Strelkoff 
oints out, this resulted in some inexactitude of statement in regard to the 
urbulence production terms. The expression calculated for P; represents 
mly part of the energy transformed from the mean flow into turbulence. The 
-omplete turbulence production relation is the product of Reynolds stresses 
ind mean velocity gradients. In tensor notation this may be written as 


‘or plane two-dimensional flows the nine terms implied by the above expres- 
sion reduce to the following four: 


ee ge tr ay TOY Gee de 


this is Equation 4 of Schubauer and Klebanoff’s(2) report. In their evaluation 
f the several production terms, the third term was found to contribute 
oughly 80 percent of the total while the fourth term was negligible. From 
his indication the P, values presented in Fig. 9 should not be in error by 
nore than 20 percent. 

‘The flow situation reported on in the present paper did not represent a 
imple two-dimensional system, since the boundary layer developed along a 
onical diffuser wall. However in view of the geometry? of the experimental 
rrangement, it would seem that the above simplified production relation in 
erms of an x-y coordinate system located at the wall is approximately cor- 
ect. On this basis the four terms have been calculated and are presented in 
"ig. 10. The values of re) U/6x utilized to calculate the first term in the P, 
xpression were found by numerical differentiation of the pitot tube data(9, 10) 
etween station 12 and stations 5 inches upstream and 2.5 inches downstream. 


. Proc. Paper 1393, October, 1957, by J. M. Robertson and G. L. Calehuff. 
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. The virtual apex of the conical surface was 50.7 inches upstream of 
‘Station 12. 
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FIG.10 VARIATION IN TURBULENCE PRODUCTION 
TERMS ACROSS BOUNDARY LAYER AT 
STATION 12 


The second term was similarly evaluated with the aid of the continuity rela- 
tion oN =m Ou, 
= A OF Ox 
transverse integration of the second difference of U with respect to x. The 
result was very approximate, but as in the case analyzed by Schubauer and 
Klebanoff, this turbulence production term was negligible, being the order of 
one-half a percent of the total at any point. The general variation of the fou 
turbulence production terms is quite similar to that found by Schubauer and 

Klebanoff. The major term in this case is nearer 90 percent of the total. 

A complete presentation of the “history” theory for the shear stress was 
not included since this is adequately covered in the literature.(13) Mr. 
Strelkoff’s discussion of the relation between the lateral shear stress gradi 
and the longitudinal pressure gradient is in agreement with the theory. As — 
noted in the paper, this relation applies near the wall and is so included in t 
theory. In the outer region the “history” effect enters and the lateral shear 
stress gradient is assumed to be that occurring at the start of the pressure 
gradient. All of the terms in Fig. 6 are defined in the paper except Two, 

e 


. The last term involves OV ox which was estimated from 


This represents an effective (hence the last subscript e) wall shear at the 
initial station. It was taken as 1.2 times the actual value of the shear in ord 
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to define the initial slope of the shear stress in the outer portion of the 
boundary layer. It is this slope which the history theory assumes invarient 
as the boundary layer progresses into the pressure gradient. 

The turbulence intensities can be presented as any of four dimensionless 
ratios; three of these were employed in the paper. Mr. Strelkoff suggests 
that if we had utilized the fourth w' / U, that the peak would have appeared at : 
each station. This peak is clearly evident in the contour plot of Fig. 5 and is — 
seen to drift to a location near the middle of the boundary layer as the layer 
develops. Nothing new in this regard would result from a presentation of 
u' /U, instead of U' /Ujg. The ratio u' /U-r was utilized (Fig. 8) only for the © 
data a station 12. As Mr. Strelkoff suggests, such a presentation of the ne 
at all stations is quite informative, as shown in Fig. 11. A progressive devia= 
tion from the pipe results of Laufer is apparent even in the wall region (inner 
tenth of the boundary layer). The measurements for the first station seem in 
good agreement with the pipe data. The turbulence levels in the outer 
boundary layer region appear very high indeed since end bear no relation to — 
the local wall shear stress and hence U7. ' 

The zero, or wall reading, of the hot wire was obtained by adjusting the 
traverser until the wire was a small, but known distance from the surface, a 
noted with a scale with a least reading of 0.05 inches. No corrections to the 
hot wire readings were applied due to wall proximity. The heat loss to the 
wall probably had some effect on the mean velocity readings. Effects on 
turbulence indications are possible but the authors are not aware of any in- 
formation on this subject. Laufer’s measurements were in a brass pipe and 
a wooden channel while Klebanoff’s boundary layer was on an aluminum plate, 

The comments by Mr. Strelkoff are appreciated. The authors are in 
complete agreement with him as to the desirability of further work along th 
line. The flow in a similar ten degree diffuser is under study at the Univer- 
sity of Illinois but no hot wire anemometer measurements have been 
attempted in it yet. 
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100 FREQUENCY CURVES OF NORTH AMERICAN RIVERS# 
Closure by E. Kuiper 


E. KUIPER,! M. ASCE.—The author appreciates very much the discussion 
his paper by Messrs. Riggs, Alexander, Kohler, and Powell, and believes 
at the points they have brought out have added materially to its value. The 
ithor was particularly gratified that the main thesis of the paper “that fre- 
lency curves of maximum annual flood peaks can be drawn with greater 
scuracy when due attention is paid to other frequency curves on the same or 
ymparable streams” was found acceptable. 

Mr. Riggs is to be thanked for his analysis of 17 stations, deriving their 
*incipal floods from rainfall, and pointing out that a very close relationship 
ay exist between the one per cent flood as an independent variable and the 
ean flood and average discharge as dependent variables. This would seem 
indicate that drainage basin characteristics, including the prevailing cli- 
ate, are sufficiently expressed in the mean flood and average discharge, 

id do not have to enter the relationship as separate parameters. This would 
mplify the problem considerably. The author hopes to find time in the near 
ture to pursue this thought and test it on other streams, possibly with dif- 
rent functional relationships. 

Mr. Alexander has attempted to make a theoretical analysis of the thesis 
hat an increase in the size of the drainage area tends to flatten the fre- 
lency line” and finds that actual analysis of existing frequency curves is 
quired. He further states that “such analyses now being undertaken by the 
, S. Geological Survey do support this thesis. The author and some of his 
ssociates, at one time, made a similar attempt to provide theoretical proof, 
‘combining assumed frequency curves by graphical procedure, and intro- 
icing correlation and coincidence of flood peaks. The conclusion of this 
alysis was that frequency curves do flatten indeed with an increase in drain- 
te area. This analysis was not included in the paper since it is rather in- 
ilved, and since there is ample and conclusive proof in the data presented in 
able 7. 

Mr. Kohler made a very valuable contribution to the paper by presenting 
s chart, relating the one per cent flood, the mean annual flood and the size 

the drainage area of the 100 river stations under discussion. The author 
us particularly impressed by the close relationship of the mean annual flood 
d the one per cent flood, regardless of the size of the drainage basin. For 
proximate estimates of the one per cent flood in cases where the mean an- 
al flood is known, the line shown in the chart will serve a very useful pur- 
se. In cases where the mean annual flood is not known, the lower part of 


Proc. Paper 1395, October 1957, by E. Kuiper. 
Chief Engr., Manitoba Water Resources Investigation, Manitoba, Canada. 


1832-60 HY 5 7 October, 1958 


the chart would provide guidance. One could select a line representing river 
stations with similar drainage basin characteristics, enter the diagram with © 
the given size of the drainage area, find the mean annual flood and the one 
per cent flood. Such a chart would undoubtedly gain in value when more in- 
formation is added. 

Mr. Powell is quite right in pointing out that the median annual flood is 
less susceptible to an unbalanced sample than the mean annual flood. How- 
ever, most of the work on frequency curves that has been done in the past 
been in terms of the mean annual flood. In order to facilitate comparison of — 
past and present data, the author has adhered to this parameter. Mr. Powell 
suggests that a formula may be developed for estimating the one per cent 
flood directly. This suggestion could very well materialize if one pursued 
valuable contribution of Messrs. Riggs and Kohler and tested them against 
available data in various drainage basins. 
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THE HYDRAULIC DESIGN OF STILLING BASINS: 
HYDRAULIC JUMPS ON A HORIZONTAL APRON (BASIN I)# 


Closure by J. N. Bradley and A. J. Peterka 


J. N. BRADLEY! and A. J. PETERKA2 —In discussing the length of jump 
is measured by the writers, Mr. Thomas and Mr. Elevatorski raise several 
yuestions. Mr. Thomas feels that the writers used an arbitrary point for 
ump length determination and that the point used in each test was not speci- 
ied, i.e., whether the point on the roller or the point where the jet began to 
eave the floor was the determining factor. In the section “Forms of the 
dydraulic Jump,” the determining factor is discussed in terms of the Froude 
lumber ranges illustrated in Figure 9. Thus, for any particular test the de- 
ermining factor for the length of jump criterion is reflected in Fy, Column 
10, Table 1. 

Mr. Thomas feels that it would have been preferable to record the length 
yf roller as the length of jump in all cases. Mr. Elevatorski voices some- 
vhat the same opinion when he states that he prefers the definition” * * * the 
0int where no return flow is visible.” Had these points been used, the jump 
engths for Froude numbers less than 4.5 would have been too short since for 
Froude numbers up to 4.5 the point at which the high velocity jet begins to 
eave the floor is the controlling factor and the roller is relatively short. 
}erious downstream erosion might result if a stilling basin were constructed 
mly as long as the roller. The bed was smooth at all times so the measure- 
nents could not have been affected by the roughness of the bed as Mr. Thomas 
issumes. The writers believe that the length determinations given in the 
laper are not based on arbitrary conditions, as stated by Mr. Thomas, but 
‘rather provide a more realistic evaluation of the length of the jump than has 
ieretofore been presented. The writers devoted several days to observing 
umps of various sizes and Froude numbers before formulating criteria for 
letermining the length. The criterion suggested by Messrs. Thomas and 
tlevatorski applies to only one type of jump. The writers wish to emphasize 
gain that the length of jump was chosen from a practical viewpoint; the end 
f the jump represents the end of the concrete floor and training walls of the 
shortest stilling basin possible. In view of the practical nature of the tests, it 
yas not considered important that the jump length was not” * * * easily ob- 
erved, and a characteristic parameter of the jump * * *” as stated by Mr. 
thomas. It was considered important, however, that the jump length be 
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chosen to represent the shortest stilling basin possible, consistent with good 
performance and a reasonable factor of safety. 4 

The experiences of Mr. Thomas with canal falls in India are interesting 
and it is possible that rapid side expansions in the stilling basin were the ; 
cause of lateral oscillations resulting in poor performance. In the writers’ 
experience, wall angles of 6° to 8° provide safe divergence. Larger angles { 
may produce unstable conditions caused by a contraction of the flow width in : 
the basin rather than an expansion. However, baffle piers help to increase 
the allowable angle of divergence but no data are available as to the degree. J 

Mr. Thomas states that insufficient evidence is presented to conclude that 
baffle piers are of little value in a stilling basin when Fy is between 2.5 and 
4.5. Paper 1404 covers this problem thoroughly and Figure 21 shows a few 
of the many arrangements which were tested in an attempt to use baffle piers. 
In this low Froude number range, conventional baffle piers do more harm ' 
than good, and they cannot be used to stabilize the flow or reduce the length 
of basin required. Conventional baffle piers were found to create instability 
and it was necessary to lengthen the basin to obtain flow conditions equivalent — 
to those in a basin without them. The only effective arrangement found by the 
writers is shown in Figure 22. Flow currents passing over the tops of the 
oversized chute blocks tend to reinforce the roller to provide a more stable 
jump. 

Mr. Elevatorski has indicated other methods of expressing the writers’ 
data, but the writers do not see an advantage in altering the method of pre- 
sentation. Before a designer completes the design of a stilling basin, he will 
have need for the value of Vj, and calculation of the Froude number is as ea 
as the calculation of the critical depth. The Froude number is used several 
times in the design of a basin using the writers’ methods. The writers also 
prefer expressing the jump length in terms of Dg rather than in terms of D2 
Dj as suggested by Mr. Elevatorski and Mr. Thomas, because in practice it 
has been found to be convenient. The writers were aware of the relations 
described by Mr. Elevatorski since they were first worked out by Riegel and 
Beebe in 1917.* 

Mr. Elevatorski’s concern over the lower part of the curve in Figure 5 is 
not understandable by the writers since, obviously, there is very little choice 
in the manner in which the curve may be extended from Fy = 2 (the last test 
point) to Fj = 0. Further, there is no mention or inference that an undular- 
type jump will be formed ‘at Froude numbers less than 2. On the contrary, 
Figure 9 shows that the prejump type of action will occur. The writers sho 3 
have stated, however, that the portion of the curve from Fj = 1 to F, = 0 has 
no particular significance. 

The writers concur with Mr. Elevatorski’s comment that the submitted 
information should not be used to design trapezoidal stilling basins, unless, 
course, the increased area introduced by placing a batter on the side walls is 
a very small percentage of the total cross-sectional area. 


i 


* The Hydraulic Jump as a Means of Dissipating Energy, Technical Reports 
Part III, the Miami Conservancy District State of Ohio, Ross M. ee an 
John C. bores 
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Corrections to Proc. Paper 1401 


Correction 
Change forms to form in line 5 


Line 4 of the Introduction—-Change Matske to Matzke 


Item 4 should read “No particular difficulty is encountered in the form 
shown in Figure 9c. Arrangements of baffles and sills will be found 


valuable as a means of shortening the basin. This is discussed in 


Papers 1402 and 1403. 
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HYDRAULIC DESIGN OF STILLING BASINS: 
HIGH DAMS, EARTH DAMS AND LARGE CANAL STRUCTURES (BASIN II)@ 


Closure by J. N. Bradley and A. J. Peterka 


J. N. BRADLEY! and A. J. PETERKA.2—The fact that the titles of Papers 
1402 and 1403 contain suggested uses of the basin should not influence the de- 
signer in his choice of basins. The only reason for using a Type IJ basin 
rather than a Type III basin is the possibility of cavitation damage resulting 
from the use of baffle piers. 

As Mr. Thomas suggests, there is no doubt that a row of baffle piers is a 
more effective energy dissipator than a row of chute blocks. However, if baf- 
fle piers are to be used there is no advantage in installing them in a Type II 
oasin. Rather, a Type III basin should be used and the resultant savings from 
reducing the basin length and simplifying the end sill can be realized. In 
Paper 1403, Introduction and Recommendations, the limiting discharge and — 
velocity for use of a Type III basin are given. 

The length of basin in terms of scour downstream can be clarified in the 
original paper by pointing out that Table 2 (which was printed in Paper 1403 
jut which should have been printed in Paper 1402) lists some 36 basins that 
were model tested using a sand bed downstream in which extensive erosion 
ests were made. Since the recommended Type II basin is patterned after 
hese, and since some of the verification test models contained a sand bed, it 
iS known, but was not stated in the writings, that the recommended structures 
will be satisfactory for installation where the scour problem is of average 
lifficulty. Figure 12 shows that the majority of existing basins similar to 
Type Il are shorter than the recommended Type II basins. Since the scour 
matterns were satisfactory for the short basins, the scour patterns for the 
onger basins will also be satisfactory. 

The drawing in Figure 10, criticized by Mr. Thomas, is not a recommended 
lesign. Rather, it is, as the drawing title states, a “Definition of Symbols” 
ind was intended to be used with Table 2. Some of the early stilling basins 
iad the end sill upstream from the end of the apron. This is not now con- 
sidered good practice. Figure 14 shows the recommended proportions of the 
[ype Il basin. 

The direct solution proposed by Mr. Thomas in Figure 1 for determining 
pron elevation may be used satisfactorily. However, the trial and error dif- 
iculties mentioned by Mr. Thomas in connection with the use of Figure 11 will 
ye found to be minor in nature when solving practical problems. The change 
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' overall error in not recalculating after repositioning the apron would be less 
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in Dy, Vy, and Fy occasioned by changing the estimated elevation of the apron 
over the usual possible range will be found to be relatively small. A single 
correction is sometimes necessary, but this is usually done only to be cer- 
tain that the first assumption was correct. For example, if a spillway dis- 
charges 14 cfs per foot with a 2.5-foot head on the crest over a fall 50 feet 
high, Fj = 14.13. If the apron elevation was assumed to be 5 feet too high and 


it is necessary to recalculate for 55 feet of fall, Fy = 14.28. would then 


TW 
Di 
change from 19.4 to 19.8 and Dg would change from 6.05 feet to 6.12 feet. 
than 1 inch. In the same problem, if the discharge per foot of width is in- 
creased to 41 second-feet the error would be less than 2 inches. Both errors 
are negligible. 


Corrections to Proc. Paper 1402 
Page Correction 


4 Under Verification Tests there is a reference to Table 2. Table 2 wa 
printed, however, in Paper 1403. Table 2 should be printed adjacent 
this first reference. 


7 Inthe section Aids in Computation there is a reference to Figure 15. 
Figure 15 should be printed here rather than in Paper 1403. 
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HYDRAULIC DESIGN OF STILLING BASINS: SHORT STILLING BASIN FOR 
CANAL STRUCTURES, SMALL OUTLET WORKS, 
AND SMALL SPILLWAYS (BASIN III)a 


Closure by J. N. Bradley and A. J. Peterka 


J. N. BRADLEY! and A. J. PETERKA.2—The discussion by Mr. Thomas ° 
follows two main trends. In one he has preferred to express the writers’ 
variables in other terms, perhaps more familiar to him. This is a justifiable. 
procedure and the choice is up to the individual. In the other he has inter- 
preted certain criteria contrary to the meaning and use intended by writers. 
In these, the writers will explain the original concepts further to insure that 
other readers will understand the true significance of the data and analysis 
presented. 

Mr. Thomas feels that the use of his Figure 1 simplifies the vertical place- 
ment of the apron by eliminating the cut and try process. As demonstrated 
in a sample problem in the writers’ closing discussion of Paper 1402, the 
writers believe that in most installations the errors involved in a first rea- 
sonable assumption are small and may be eliminated with one relatively sim- 
ple correction. They prefer this method, when using the overall procedures 
given in the paper, to the one explained by Mr. Thomas which requires the 
calculation of additional values not otherwise used. 

The necessary depth, Dg, for a Type III basin is given in Figure 11 and 
this is shown to be 15 to 18 percent less than the conjugate depth computed 
for a hydraulic jump, depending on the Froude number. In the paragraphs 
under the heading “Tail-water Depth” it is recommended that full conjugate 
depth be used for a Type Il basin. The reasons are stated. Mr. Thomas 
states in Item 1 of his discussion that the basin floor may be raised because 
of the greater factor of safety against jump sweepout in a Type III basin. If 
only the arithmetic in the problem is considered, the basin floor may truly be 
raised. However, in the light of past experience with prototype structures 
where operating difficulties are often experienced because of a lack of tail- 
water depth and rarely, if ever, because of too much depth, the writers de- 
cided to recommend full conjugate tail water. The basin then contains a factor 
of safety against phenomena of which the designer is often unaware such as 
retrogression of the stream bed or a poorly derived tail-water curve. In 
addition, full conjugate depth will provide improved performance including 
general appearance, velocity distribution, wave and surge formation, and 
scour. Safe, but minimum, performance will result if the designer chooses 
to use the minimum tail-water value shown in Figure 11. 
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Mr. Thomas states that the length of jump is more directly a function of 
D; = Dg - Dy, than of Dg. The writers chose D2 because they found it more 
convenient to work with one variable rather than two. Either method is ac- 
ceptable however. 

Mr. Thomas states that the length of basin required is not necessarily the 
same as the length of jump as measured by the writers. The basin length 
measured by the writers, as stated in Paper 1401, is as follows: 


“It was the intention to judge the length of the jump from a practical 
standpoint; in other words, the end of the jump as chosen would represent 
the end of the concrete floor and walls of a conventional stilling basin.” 


The jump length in Basin I and the basin lengths for Type Il and Type Il 
were all judged on the same basis. For similar flow conditions, Basin I, ; 
Basin II, and Basin III will provide the same degree of protection against q 
scour and in many cases will provide nearly identical overall performance : 
based on the resulting downstream flow conditions. Representative erosion t 
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tests were made (but not reported) in sufficient numbers.to provide assurance © 
that the recommended basin lengths were sufficient to prevent undermining 
of the end sill and excessive scour downstream. In fact, scour was one rea- 
son for not reducing the basin lengths still further. Complete scour test data 
were not recorded because of the difficulties in obtaining and presenting 
meaningful erosion data. The writers are confident that the scour problems 
which are encountered at the ends of training walls and other areas which 
cannot be evaluated in a sectional model will be of greater concern toa de- 
signer than the scour encountered downstream from the recommended stilling” 
basins. On the other hand, the writers believe, and have advised in several 
places in the group of papers, that model tests should be performed to check ; 
the recommended designs whenever unusual conditions are encountered on a 
particular structure. Certainly, an easily eroded downstream channel bed 
should be included in this category. 

The position of the baffle piers on the apron is given by the writers as 
0.8 Dg. It was found that if the piers were placed closer to the toe of the 
jump, the individual jets of water rising vertically from each pier penetrated 
the covering tail water resulting in incomplete stilling action. As the piers 
were moved downstream from the 0.8 Dg position, the usefulness of the piers 
in creating turbulence decreased rapidly. If the basin entrance velocities are 
sufficiently high, cavitation at the baffle piers is most likely to occur when thi 
piers are in the upstream position. However, data taken on one set of piers 
since Paper 1403 was written indicate that moving the piers downstream does 
not reduce the possibility of cavitation as much as might be expected. For 
Fy = 10, q = 240 cfs, and Vq = 100 feet per second, the lowest pressure mea- 
sured on a sharp-edged baffle pier placed 0.64 Dg downstream from the start 
of the jump was -20 feet of water. When the pier was moved to 0.98 Dg, the 
pressure increased but only to -10 feet of water. In the downstream position, 
the effectiveness of the piers was reduced and the basin did not perform as 
well as a standard Type III basin. 

Moving the baffle piers farther downstream resulted in greater loss of per: 
formance. Streamlining the piers and placing them in the upstream position 
resulted in considerable loss in effectiveness. In fact, when the piers were 
streamlined sufficiently to produce acceptable pressures, the turbulence- 
creating qualities were so reduced that it was doubtful whether the use of the 
streamlined piers could be justified. The piers were less effective than 
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ndicated by the area of the vertical face. Because of the many possibilities 
resent as illustrated above, no flexible arrangement of baffle piers could be 
‘iven in the original paper. 

The limiting factors concerning the use of baffle piers placed 0.8 Dg down- 
itream from the toe of the jump are given in the paper as Vj = 50 feet per 
econd and q = 200 cfs per foot of basin width. These values contain a judg- 
nent factor by the writers, but until more data are available the stated limits 
thould be used and the baffle piers should be placed at 0.8 Do from the chute 
locks, regardless of the absolute size of structure. Mr. Thomas suggests, 
‘* * * in a large work, blocks must be placed at a greater relative distance 
rom the toe than in smaller works with equal Fj.” This is not so. If hy- 
lraulic models can be accepted as a means for determining stilling basin 
limensions, the absolute size of structure has nothing to do with the position 
if baffle piers relative to the toe of the jump. The writers’ tests apply to 
arge structures as well as small ones. 

Mr. Thomas has failed to realize the significance of Figure 19. The pro- 
ile represents the approximate water surface profile for a Type III basin 
yperating at conjugate depth. It does not show, as Mr. Thomas states, “* * * 
he tail water had been lowered to the stage when the beginning of the jump 
iad receded from the toe and reached the blocks * * *” 

Mr. Thomas feels that the writers dismissed the possibility of using a 
second row of baffle piers without due consideration. As stated in the paper, 
ests with two rows of baffle piers were made for each arrangement of baffle 
jiers tested. Although the second row of piers did improve the performance, 
he improvement was not sufficient to allow reducing the basin length further. 

There is no doubt that in the writers’ tests the chute blocks created tur- 
uulence in the upstream portion of the basin, helping to make the second row 
wf baffle piers relatively ineffective. The writers also believe that where an 
srosion problem is particularly difficult, a second row of piers might help in 
‘educing scour. It would then be necessary, however, to lengthen the basin 
yy the distance between the rows of piers to maintain the proper distance be- 
ween the piers and the end sill. For the general case, however, in the judg- 
nent of the writers a second row of baffle piers in Basin III cannot be justi- 
ied. 

-Mr. Thomas expresses a desire for more supporting data in order that the 
lesigner may judge the effects of deviations from the standard Basin III. 
‘ince this is a new design and had not been verified by prototype experiences, 
he design was kept rather rigid. However, studies are being continued to de- 
ermine the pressures on the baffle piers for variations in the 0.8D2 dimen- 
ion, the effect of baffle pier streamlining on the baffle pier pressures, and 
ther features of the Type III basin presently covered in a general way. How- 
ver, the data presented are sufficient to provide a first design which will 
rovide satisfactory performance with a minimum structure. Refinements or 
urther reductions in dimensions should be accomplished with hydraulic 
nodel tests. 
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Corrections to Proc. Paper 1403 
Page Correction 
5 Figure 18 in the last line should be changed to Figure 17. 


12 and 
13. Table 3 should be printed in Paper 1402. 


16 In the last sentence on Baffle Piers, chambers should be changed to 
chamfers. 
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HYDRAULIC DESIGN OF STILLING BASINS: 
STILLING BASIN WITH SLOPING APRON (BASIN Vv)4 


Corrections to discussion by A. Rylands Thomas, ! M. ASCE 


CORRECTIONS—The accompanying figure was inadvertently omitted from 
Mr. Thomas’ discussion of Proc. Paper 1405 in Proc. Paper 1616 (April, 
1958): 
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. Proc. Paper 1405, October, 1957, by J. N. Bradley and A. J. Peterka. 
. Consultant, London, England. 
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THE HYDRAULIC DESIGN OF STILLING BASINS@ 


Closure by J. N. Bradley and A. J. Peterka 


J. N. BRADLEY! and A. J. PETERKA.2—Mr. Nimmo describes Sommer- 
et Dam which has a relatively complicated discharge system. The data of 
jasin II, therefore, should not be expected to apply in the sense that a basin 
esigned from the Basin II rules will operate properly with nonuniform flow 
istribution. On the other hand, if model tests were used to provide sills, 
tc., to distribute the flow so that the overall effect was the same as uniform- 
7 distributed flow entering the basin, Basin II criteria may then be used to 
heck the overall size, depth, and other dimensions of the basin. This, Mr. 
immo has done and found that Basin II would have been slightly larger than 
1e basin developed from individual model tests. This helps to substantiate 
1e authors’ statements that there is a factor of safety inherent in the recom- 
1ended basins, I-VI. 

On Tinaroo Falls Dam, the stilling action described by Mr. Nimmo re- 
ulted from impact of the flowing water on the apron and on the sills. It is 
ften possible to provide impact-type structures which are smaller than hy- 
raulic jump basins. However, impact-type structures must be tailored to fit 
1e exact flow conditions encountered and judgment and experience are also 
equired to provide a satisfactory stilling structure. 

Mr. Rand feels that more flexibility should have been provided so that the 
esigner could vary the dimensions of the basin to fit individual cases. An 
nalysis which would completely evaluate each variable would provide an 
ileal solution but would have required a testing program too involved and 
sngthy to complete. 

To obtain the most usable data from a reasonable expenditure of time and 
inds, the authors evaluated the major variables, basin length and depth, in 
rms of fixed arrangements of blocks, piers, and sills placed to provide the 
mallest stilling basin consistent with good performance. The arrangements 
10Sen resulted in the least tail-water depth and/or basin length. Thus, if 
ore flexibility in the appurtenances had been provided, a larger stilling 
isin would have been necessary for each arrangement. The major basin 
imensions, depth and length, usually determine whether the designer can 
sonomically justify the use of a hydraulic jump stilling basin. If he can, he 
as only to place the recommended appurtenances in the basin to obtain opti- 
ium basin performance. Thus, there is no primary need for alternate baffle 
er sizes or locations or other minor variations. 


_ Proc. Papers 1401 through 1406, October, 1957, by J. N. Bradley and A. J. 
Peterka. 

_Hydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly 
of U. S. Bureau of Reclamation, Denver, Colo. 
Hydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 
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As stated in the papers, the dimensions and spacing of the blocks, piers 
and sills are not critical. Departures from exact dimensions may be made 
to avoid construction joints, fractional blocks, etc. Thus, a complete analysis 
of the basin would have served mainly to prove that the blocks, piers and 
sills were of optimum size and in the best arrangement. As the material is 
presented, however, it is necessary for the reader to rely on the judgment of ‘ 
the authors to some degree, although the data contained in the tables show the 
range of the investigation concerning other arrangements of basin appurte- 
nances. g 

Despite Mr. Rand’s objections to the authors’ use of the word 
“generalized,” the authors feel that they have generalized the design of 
several types of stilling basins, since data are presented which make it pos- ~ 
sible to design a stilling basin for most any combination of discharge and : 
height of fall. To provide generalization to the degree suggested by Mr. Rand 
would require an almost impossible amount of testing and data preparation ' 
and would not provide more usable designs. A more reasonable approach to 4 


: 
Ly 
€ 
; 
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the problem is to generalize the design of other types of basins which ex- % 
perience shows to be in demand. When more of these have been completed, 5 
the designer will have a greater variety of basins from which to choose. ; 


Mr. Rand poses several questions concerning the jump length and the basin : 
length as measured by the authors. These can probably best be answered by 
reviewing the basic objectives of the test program. In Basin I, the length of 
paved apron is the same as the length of the jump. In Basin II, chute blocks 
were added to the basin to promote stability in the jump and to increase the 
factor of safety against jump sweepout at tail-water depths less than conju- 
gate: The end sill was then moved upstream on the paved floor until the flow 
conditions downstream from the dentated sill were just beginning to become 
objectionable with regard to scour and water surface roughness downstream. 
In this position, the high point of the boil was at or upstream from the end 
sill; in other words, the water surface did not continue to rise beyond the end © 
of the basin training walls to allow inflow into the basin from the true tail- . 
water surface behind the training walls. The jump length was found to be les 
for Basin II than Basin I, although the action at the ends of both basins is sim: 
lar. The end sill was, therefore, not subjected to a particularly high velocity 
current and its primary purpose was to lift upward the remaining bottom cur- 
rents to reduce scour. The jump is therefore “forced” only to a slight degree 
by the end sill. Some of the reduced length is produced by the chute blocks in 
the upstream end of the basin. 

For Basin III, the same procedure was followed with chute blocks and baf- 
fle piers in the basin. After the baffle piers had been moved upstream as far 
as possible, without jets of water rising from the piers to penetrate the water 
surface, the end sill was moved upstream to the position which produced the 
same degree of performance at the end sill obtained for Basin II. In this cas 
the baffle piers “force” the jump to occur in a shorter distance. The end sill 
in Basin III is therefore not a necessary item in the sense that it forces the 
jump to occur. Only when the tail water is lower than conjugate depth does — 
the sill help to hold the jump in the basin. 

Since the sill is subjected to greater flow action in Basin II than in Basin 
III, the dentated sill dimensions are more critical than the flat sill of Basin 
Ill. The height of the sill in Basin III is not critical; however, the sloping fac 
should be sufficiently long to provide directional effect to the flow leaving the 
basin. A steep sill tends to produce a more violent ground roller which may 
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ft bed material into the current. Fine material will then be carried away by 
le upper water currents. A flatter sill produces a mild ground roller and 
ear neutral conditions in the bed material just downstream from the end sill, 
esulting in minimum scour and surface waves in the downstream channel. 
he recommended end sill therefore was chosen to provide optimum per- 
yrmance, a smooth water surface, and minimum scour. 

In Figure A, Mr. Rand shows a range of basin lengths obtained by the 
uthors and others and relates them to a coefficient K. He summarizes the 
esults by stating “* * * there is no basin length that is the only correct one.” 
he authors believe that there is only one correct length—the shortest one 
iat can be used safely. This the authors have shown in their length charts. 
.§ stated in the “Introduction” to Paper 1401, many tests have been made to 
etermine basin dimensions, but when groups of tests which combine the 
idgment and requirements of several experimenters are lumped together, 
1e conclusions are often misleading and sometimes contradictory. Since 
idgment cannot be eliminated from a hydraulic jump testing program, it is 
ecessary for the reader to decide from the description of the tests whether 
1e results are dependable. If they are, he should accept the findings entirely. 
ny attempt to mix the features of a stilling basin recommended by one group 
f experimenters with those recommended by another will lead to difficulty. 

The average tail-water depth 0.99 Dg, given in Table 2, was not the reason 
or recommending full tail-water depth, Do, for Basin III. It was coincidental 
jat the recommended value and the average value in Table 2 were the same. 
‘ests on Basin III showed that performance was far better at conjugate depth 
1an at the minimum tail-water depth, 0.83 Dg, shown in Figure 11. The rea- 
ons for using full conjugate depths in Basin III are fully explained in Paper 
403 under the heading “Tail-water Depths.” 

Mr. Rand asks the question “Could not more flexibility (in sill height) be 
rovided?” The authors found in these tests and in tests of other stilling 
asins not reported that higher end sills, particularly those with vertical 
ices, tend to increase the wave problem and do not help the scour problem. 
ills which tend to increase the ground roller action circulate more of the 
roded material near the sill. Material is continually being packed against 
1e sill with the excess being carried downstream to be recirculated. The 
uthors believe that the grinding action will damage the concrete by abrasion. 
ills which are too low produce greater scour depth close to the end of the 
pron. The sill height therefore was chosen to produce sufficient lift to re- 
ace scour but not enough to produce grinding or waves in the downstream 
hannel. 

Mr. Blaisdell shows that the authors’ energy loss curves based on their 
sst data fit the theoretical energy loss curves well. Since the basin lengths 
yr Basins II, III, and IV were chosen so that practically all of the jump action 
scurs within the basin, the loss curves of Figures 8 and 8A may be applied 
) Basins I through IV. In Mr. Blaisdell’s SAF basin, Figure 8, much of the 
stion occurs downstream from the end sill, leaving a considerable portion of 
1e hydraulic losses to occur in the downstream channel. 

During the testing program to develop Basins I through IV, the authors 
Instructed and tested several SAF basins. In their judgment and that of other 
2sign engineers in the Bureau of Reclamation, the overall factor of safety in- 
rent in the SAF design was too small for the relatively large structures de- 
igned in the Denver office of the USBR. In fact, Mr. Blaisdell states in his 
aper on the SAF basin “the safety factor incorporated in the design of the 
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SAF stilling basin is low, but it is felt to be adequate.” For the relatively 
small structures for which Mr. Blaisdell designed the SAF basin this is true. 

Experiences in the field with full-sized structures have demonstrated to 
the authors that tail-water elevation is often the most unpredictable feature 
in a basin. Experience has shown that the tail water may be several feet 
lower than anticipated because of riverbed retrogression or for other reaso 
or that the time required to establish the tail-water elevation may be longer 
than anticipated. Considering these factors, it is necessary that a stilling 
basin contain a built-in factor of safety against sweepout which is greater than 
the probable error in determining tail-water elevation. In remote areas, it 
is sometimes necessary to compute the tail-water elevation from incomplete 
records. In these instances, it is essential that a greater factor of safety 
than provided in the SAF basin be used. The choice between the two basins, 
therefore, may be a matter of judgment despite the fact that extensive data 
are presented to prove the value of both the SAF and Basins I-VI. 

The number of baffle piers to be installed in Basin III is not as critical as” 
Mr. Blaisdell indicates in his example. Changing the percent of space occu- 
pied by the baffle piers in Basin III by several points would not adversely af- 
fect performance and this leeway gives the designer the opportunity to obtain 
the desired number and arrangement of piers. 

The authors would like to thank the discussers of this series of papers for 
their interest and helpful suggestions. Many of the comments made will be 
value in future investigations of this type and in the continuation of this 
development program. 


Corrections 
Paper Page Correction 
1404 7 Near the end of the second paragraph, the word reduc- 
tions should be changed to reduction. 
1405 20. The name D. D. Handlanb should be B. D. Hindland. 
1406 6 Bibliography #5 the name Kinsvater should be 


Kindsvater. 
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HYDRAULIC DESIGN OF STILLING BASINS: HYDRAULIC JUMP ON A 
HORIZONTAL APRON AND HYDRAULIC JUMP ON A SLOPING APRON2 


Closure by J. N. Bradley and A. J. Peterka 


J. N. BRADLEY! and A. J. PETERKA.2—Mr. Vega has used the authors’ 
data to attempt to substantiate equations derived by Mr. F. J. Dominquez and 
to present an alternate method of expressing the length and depth of hydraulic 
jumps on sloping aprons. Although the general mathematical theory appears 
to be in order, it was necessary for Mr. Vega to make certain assumptions to 
convert the basic equations into usable forms. Apparently, the assumptions 
were not entirely valid, since the results actually measured by the authors 
cannot be reproduced from the final equations and charts developed by Mr. 
Vega. According to Mr. Vega, the percentage of error of his equation (8) 
from the actual values measured by the authors was 5.07 percent; however, 
individual values varied as much as 17.7 percent as shown in his Table 1. An 
inspection of his Figures 4 and 5 shows that the curves only approximate the 
points, particularly the curve of Figure 5. In combining the two curves, Mr. 
Vega has, in effect, averaged the jump lengths on all sloping aprons and has 
eliminated the differences in jump length found by the authors for various 
apron slopes and shown in Figure 32. Mr. Vega’s statement that “* * * the 
length of the jump is independent of the slope * * * ” is therefore incorrect 
and the curve of his Figure 6 gives only an approximate answer. Similarly, 
the curves of Figure 7, which give the initial and final depths before and after 
the jump, are also approximate. Over the midrange of Xo values, the curves 
are probably sufficiently accurate for ordinary problems but at low values of 
Xo (which correspond to high values of the Froude number) the curves indicate 
2 deficiency in tail-water depth. 

The following example illustrates the typical differences in jump dimen- 
sions to be expected from the authors’ and Mr. Vega’s methods. Assume a 
spillway with an ogee overflow section to have a discharge coefficient of 3.8 
und a unit discharge of 300 second-feet. The head on the crest would be 18.4 
eet. Give the jump dimensions for fall heights of 50, 100, and 250 feet for 
ypron slopes tan @ of 0, 0.05, 0.10, 0.15, 0.20, and 0.25. The Froude number 


F for the three heights of fall is 4.15, 6.55, and 11.38, respectively. Using the - 


zuthors’ methods along with Figures 15, 31, and 32, the values in the table be- 
ow may be computed for the three Froude numbers given above. The values 
n the table are plotted in Figure 47. The jump length varies considerably, 


1. Proc. Papers 1401 and 1405, October, 1957, by J. N. Bradley and A. J. 
Peterka. 

|. Hydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly 
of U. S. Bureau of Reclamation, Denver, Colo. 

. Hydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 
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Jecoming longer as the apron slope increases. This is shown by the lower 
yroup of curved solid lines which are loci of jump lengths for the range of 
slopes shown. The necessary tail-water depth also increases as the apron 
slope increases. The required depths are shown as the upper group of solid 
ine curves, plotted above the apron to which the depth applies. 

Using Mr. Vega’s methods, the results are somewhat different. For iden- 
ical flow conditions, Xo is 0.387, 0.288, and 0.197. Figure 6 shows the jump 
ength to be the same for all apron slopes. The jump lengths for the three Xo 
values have been plotted on Figure 47 (as the three vertical dotted lines) for 
-omparison with the authors’ findings. For each Fy or Xo value, the jump 
.ength determined from the equation derived by Mr. Vega is too long for flat 
slopes and too short for the steep slopes. In the extreme case, the jump 
‘ength shown by Mr. Vega, is 282 feet compared to 314 feet determined from 
the authors’ methods. In addition, the tailwater depths, indicated in Figure 7 
and plotted above the apron end as found by the authors’ methods, are deficient 
xy the amount shown in the table and in Figure 47. For the extreme case 
illustrated, Fy = 11.4 or Xo = 0.197, for an apron slope of 0.25, the tail-water 
jepth found by Mr. Vega is 7 feet less than the 104.5 feet determined by the 
authors to be the minimum. 

Some of the comments made by Mr. Thomas concerning sloping aprons are 
juestioned by the authors. Mr. Thomas starts his discussion with a statement 
hat the analysis of a sloping apron requires an assumption of jump length, a 
jimension which is not precisely known. Apparently, Mr. Thomas does not 
accept the fact that the entire testing program on sloping aprons was essential- 
y devoted to measuring the jump length on various slopes. The results are 
shown in Figures 32 and 33. 

The authors’ test apparatus was sufficiently large that some air was pres- 
sant in all the jumps tested. As might be expected, the larger flumes having 
1igher velocities had the most entrained air. No actual air quantities were 
measured but it is believed that the air quantities were not sufficient to affect 
the water surface profiles. Further, the authors believe from their experi- 
ances with models and prototypes that in the usual prototype structure, the 
water surface profile will not be changed sufficiently to be of concern. The 
isual freeboard allowances made on prototype training walls will include the 
sffect of increased air entrainment. 

Mr. Thomas states that a gently sloping apron would serve as a “visual 
warning” to prevent moving the jump off the apron during gate operation and 
hat the gently sloping apron would be superior in this respect to a horizontal 
ypron. The authors are not sure what Mr. Thomas means by “visual warning” 
jut in neither case does the jump suddenly move off the apron. 

In both cases, there is a gradual movement of the jump toward the end sill 
is the sweepout condition is approached. The taller end sill or the baffle 
iers on the horizontal apron might make the horizontal apron less susceptible 
© sudden sweepout than a sloping apron. As stated by the authors, the chief 
idvantage of a sloping apron is that in some installations the slope may reduce 
xxcavation costs, and not, as Mr. Thomas states, that the sloping apron would 
serve as a “visual warning” during gate operation. 

The authors did not:use baffle piers or even chute blocks on the sloping 
ypron because experience has shown that sloping aprons are usually associ- 
ited with very high dams. The sloping apron without appurtenances of any 
cind except a small end sill is therefore not limited by entrance velocities — 
vyhich might produce cavitation damage on or adjacent to the chute blocks or 
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baffle piers. No doubt a sloping apron using appurtenances could be de- 
veloped. However, there has been no need for this type of structure in the 
authors’ work. The authors regret that Figure A was not included with 
Mr. Thomas’ discussion. 

The information presented by Mr. Beck on the Santee Cooper, Petenwell 
and Castle Rock Dam stilling basins (his Table 1 and Figure 2) is interesting. 
The results obtained from individual model studies agree quite well with the 
general curves presented on sloping apron basins. It is not possible to com- — 
prehend the work involved in proportioning a sloping apron by trial, even with 
the aid of an individual model, without actually experiencing it. It is still ad- 
visable to verify sloping apron designs for large structures with an individual 
model, but the generalized curves presented by the authors should eliminate 
all or the greater part of the trial process. 
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Notes for Table 12 


On ogee crest, discharge coefficient is 3.8, discharge per foot 
width 300 cfs, head on crest is 18.4 feet. 


For 50 ft of fall Vy = 55 ft/sec P = 5645 ft 
IT) = 4.15 Xo = 0. 387 
do = 14.08 

For 100 ft of fall Vy = 74.5 ft/sec Dy = 4.02 ft 
Fy = 6.55 Xo = 0.288 
do = 14.08 

For 250 ft of fall V, = 108 ft/sec D, = 2.78 ft 
Fy = 11.38 XO = 0.197 
ad, = 14.08 
Table 12 
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HE HYDRAULIC DESIGN OF STILLING BASINS: SMALL BASINS FOR PIPE 
OR OPEN CHANNEL OUTLETS—NO TAIL WATER REQUIRED (BASIN VI)@ 


Closure by J. N. Bradley and A. J. Peterka 


J. N. BRADLEY! and A. J. PETERKA.2—The impact-type stilling Basin VI 
las been subjected to a series of tests by Mr. Argue and its performance 
ompared with a more standard type of basin having flared side walls, a de- 
ressed bottom, and an end sill. Despite the clearly stated dimensions and 
apacities given for the impact basin in Table II, Mr. Argue has based most 
ff his comparison tests on an overloaded structure. The maximum capacity 
ff Basin VI, 13 feet wide, is given in Table Il as 191 cfs. The test discharge 
ff 330 cfs used by Mr. Argue represents an overload of 73 percent above the 
‘ecommended capacity. It is not clear to the writers how it would be possible 
O pass this extra discharge through the basin without water passing over the 
op of the baffle. Figure 11A shows little if any flow over the baffle. Also, it 
jppears in Figure 11 that the structure was set in a very wide channel with 
1O wing walls connecting the basin to the banks. Induced eddies could then in- 
luence the flow pattern leaving the basin by concentrating the outflowing jet 
nto a narrow stream. Mr. Argue has mentioned this as one of the character- 
stics of Basin VI which he did not like. When side eddies are prevented, as 
n the recommended installation, Figure 46, this type of flow does not occur. 

Mr. Argue states that the impact basin has a ‘rigid relationship between 
naximum discharge and basin dimensions and no variation in these dimen- 
sions is given with different erodible materials downstream from the struc- 
ure.” Apparently, Mr. Argue did not realize the significance of the upper and 
ower limit curves of Figure 42 which are illustrated by the photographs of 
‘igure 43. The structure width may be varied between the upper and lower 
imit values to produce varying outflow velocities. For example, for a dis- 
hharge of 191 cfs, the width may be varied from 11 to 13 feet, making the dis- 
harge per foot of width of basin vary from 17.3 to 14.7 cfs. With the alternate 
nd sill, the unit discharge may be reduced further. The relative outlet 
‘elocities may be computed (assuming uniform velocity distribution) to deter- 
nine whether the higher or lower velocities are more compatible with the bed 
naterial encountered at a particular site. 

The writers believe that Mr. Argue has placed undue emphasis on full pipe 
erformance. Under ordinary conditions on a structure of this type, the pipe 
fill run full only at the maximum discharge. Even then, it is difficult to fill 
he pipe unless a carefully designed and constructed entrance is provided. 
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(It is possible in a 4-inch-diameter model to have full pipe flow which could 
not occur in a 54-inch prototype. Field experience indicates that the head on 
the entrance must exceed four times the pipe diameter and the pipe slope 7 
must be less than 1 percent to have a full pipe). The writers’ Basin VI is 
therefore designed to give best performance with open channel flow. 

The writers have several comments regarding the basin comparison tests. 
Mr. Argue states that the most obvious differences between the two struc- 
tures was the variation in discharge; the impact basin gave 10 percent less 
maximum capacity. It should be understood that the purpose of the impact 
basin is not to produce a maximum discharge but rather to dissipate and 
distribute energy effectively fora particular discharge or range of dis- 
charges. The fact that the impact basin reduced the maximum discharge—an 
overload of 73 percent—by 10 percent could hardly be called a defect. 

The statement by Mr. Argue, “The hanging baffle had no effect until the 
flow reached 220 cfs,” is probably an inadvertently worded sentence. Itis 
difficult to visualize a discharge of 220 cfs passing through a 54-inch-diamet 
pipe at a velocity of about 20 feet per second producing no effect when it ’ 
strikes a vertical wall. In the writers’ tests, the flow striking the baffle pro- 
duced violent eddy action upstream from the baffle. Flow was forced out fro 
beneath the baffle and the end sill produced sufficient resistance to the flow t 
produce fairly uniform distribution of flow into the downstream channel. 
Figure 43 illustrates the action in the model and Figure 46 shows perform- 
ance typical of several prototype structures which have been observed in the 
field. 

It is possible that the overloading of Basin VI prevented proper action of 
the water in the limited space upstream from the vertical baffle, making it © 
appear that the baffle had little effect on the flow. The increased head above 
the baffle would also increase the velocity beneath the baffle. These effects 
could account for the flow not filling the 45 degree expansion at the end sill i 
Mr. Argue’s tests. In the writers’ tests, when the quantity of water recom- 
mended for the structure was discharged at the recommended velocity, the | 
45 degree alternate sill was a distinct improvement to the basin. . 

In regard to the comparison of scour patterns, the statement of a few 
basic facts will probably clarify the detailed comparisons Mr. Argue has 
made. There is no magic quality of an energy dissipator which can make i 
clearly superior to all other energy dissipators for all operating conditions. 
Hydraulic structures which will handle an overload usually do not work as 
well for small discharges and vice versa. A design which has been develope 
for use with riprap protection will probably produce scour if the riprap is n¢ 
used. In short, unless a structure is tested and used under the conditions fo: 
which it was developed, it may not show to best advantage. The writers be- 
lieve that Basin VI constructed and used according to the recommendations 
will outperform a basin of the same general size which does not containa 
baffle or impact surface. In fact, Figures 2A, 2B, 3A, 3B, 4A, 4B, 5A and 5E 
show the impact-type basin to be superior or emia in tperideaaeee with he 
P. W. D. structure for discharges up to or slightly above the design discha1 
given for Basin VI. 4 

Mr. Argue’s remarks concerning the need for a cutoff wall at the end of 
the basin are correct. Because of undersized riprap, one prototype structu1 
subjected to near maximum flows was slightly undermined. As a protectior 
against scour of this type in the future, the end wall of the basin has been 
tended downward (several feet, depending on the type of material) to provid 
cutoff wall. However, on structines where adequate riprap protection can b 
provided, there is no need for this additional protection. k 
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